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, this amplifier and the associated logic cir¬ 
cuits have been presented by the writer 
in another paper . 3 

The reasons for choosing serial instead 
v© Y°f parallel operation are quite simple. In 
) |a parallel system, each digit of a number 
/*! J appears in a different circuit. In a paral- 
/ ^1 a< 3der a separate adder stage is used for 
each digit, whereas in a serial adder only 
one adder stage is used, and this one 


Synopsis: The first electric digital com¬ 
puters were built around the properties of 
relays. The superior speed capabilities of 
vacuum tubes has led in recent years to 
their use in new computer designs to replace 
relays. Because of the small size, low 
power consumption, and expected long life 
of transistors, it now appears that the 
transistor will replace the vacuum tube as 
a computer element. This paper presents 
a study of binary computer functions with 
recommended mechanizations that were 
selected because they appeared to be readily 
attainable with transistors now under de¬ 
velopment. Block diagrams are presented 
of switches, memory units, arithmetic units, 
and other basic components. Estimates are 
given for the number of parts required in 
the units. It is concluded that a high- 
performance all-semiconductor computer 
can be built with germanium diodes and 
transistors. 


D IGITAL computer operations 1 can be 
divided into two classes, memory 
and logic. Memory includes the thr ee 
operations write, read, and erase, where 
the write operation, for example, can be 
defined as the representation in space of a 
function of time. In computers, the logic 
operations consist of the recognition of 
spatial distributions of voltages and cur¬ 
rents. All the logic operations that can 
be described in words can be mechanized 
with suitable combinations of three basic 
functions. The basic functions are “or,” 
“and,” and “inhibition” or negation. An 
^-terminal or-circuit will have an output 
when any one of its n input leads is ener¬ 
gized. An ^-terminal and-circuit will 
have an output only if all of its n leads are 
energized. Inhibition leads can be added 
to either of the two circuits with the result 
that the circuit cannot have an output as 
long as there is a signal on any of the in¬ 
hibition leads. These logic circuits have 
the advantage that they can be built out 
of passive nonlinearities. A nonlinear de¬ 
vice has two states, a high-impedance and 
a low-impedance state. The more pro¬ 
sper 52-162, recommended by the AIEE Com¬ 
puting Devices Committee and approved by the 
AIEE Technical Program Committee for presenta¬ 
tion at the AIEE Summer General Meeting, 
Minneapolis, Minn., June 23-27, 1952. Manu¬ 
script submitted October 31,1951; made available 
for printing April 17, 1952. 
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nounced the nonlinearity the more satis¬ 
factory the unit is as a 2-state element. 
Since nonlinearity can be obtained from 
purely passive elements, it is possible to 
obtain all of the logic functions of a com¬ 
puter without the use of vacuum tubes or 
other active elements. In practice, it is 
necessary to use some active elements be¬ 
cause the nonlinear devices available are 
not ideal and have losses. It is necessary, 
therefore, to use amplifiers to make up for 
losses in the logic circuits. It is impor¬ 
tant to note that as the efficiency of the 
logic circuits is increased, the number of 
amplifiers required is decreased. 

This computer design philosophy was 
followed in the design of the National 
Bureau of Standards Computer SBAC. 2 
It is believed that the approach that will 
result in a vacuum-tubeless computer at 
the earliest date is to follow the SEAC ex¬ 
ample, in so far as the use of germanium 
diode logic circuits is concerned, but re¬ 
placing the vacuum-tube amplifiers with 
transistor amplifiers. Since the transistor 
itself has voltage and current relation¬ 
ships quite similar to a germanium diode, 
it is expected that the diodes in a transis¬ 
tor computer will operate in a more 
natural environment than in a vacuum- 
tube computer and will respond to such 
favorable conditions by exhibiting longer 
life and more reliable operation. It is for 
such an all-semiconductor computer that 
the block diagrams described herein have 
been planned. 

Serial operation has been used in all of 
these designs. In serial operation all the 
digits of a word or number are transmitted 
in series along a single wire with the least 
significant digit first. The interval be¬ 
tween the rise of successive digit pulses is 
referred to as one digit time. The interval 
between successive numbers or words is 
commonly referred to as the word time. 
Digit times of 1 microsecond have proved 
to be practicable, and words as long as 50 
binary digits have been employed (equiv¬ 
alent to 15 decimal digits). 

A number of the block diagrams pre¬ 
sented herein have been mechanized with 
a transistor amplifier that regenerates 
pulses at a megacycle rate. This amplifier 
operates on a total power drain of ap¬ 
proximately 50 milliwatts. The details of 


stage handles all the digits. Not only are 
the arithmetic operations performed with 
fewer components in the serial machine 
but switching is also simplified. To 
switch an w-digit parallel number req uir es 
n switches to switch the separate wires, 
whereas in the serial machine only one 
switch is required. 

Basic Building Blocks 

The proposed computer elements are 
based on an assembly of a relatively small 
number of different kinds of basic uni ts 
The basic units have been limited to an 
or-circuit, an and-circuit, an inhibitor-cir¬ 
cuit, an amplifier, and a storage cell based 
on a delay line. The designs suggested are 
similar to those in the ED VAC designed 
at the University of Pennsylvania and the 
Bureau of Standards SEAC. Some of the 
basic components have been described in 
in the literature. 4 ' 5 They are described 
again in the following paragraphs to pro¬ 
vide a specific basis for the estimates of 
numbers of parts that are a goal of this 
study. 

An M-terminal or-circuit (see Figure 1) 
develops an output when any one of its 
input terminals is energized. Crystal 
diodes are put in series with each input to 
prevent a pulse at one input from feeding 
back to any of the other inputs. An n- 
terminal circuit, therefore, requires n 
crystal diodes. 

An and-circuit having n input ter- 



Figurel. Or-circuit 
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Figure 2. And-circuit 


minals develops an output only when all 
n of the input terminals axe energized. 
In Figure 2 each of the inputs is returned 
to a negative voltage and the output is 
clamped slightly below ground by Xn+i. 
Only when all of the inputs rise above 
ground will X n +1 be cut off, permitting the 
output to rise. Thus the output consists 
of the overlapping part of the inputs. An 
w-temiinal circuit is seen to require n+l 
crystals. 

Two terminal and-circuits are used ex¬ 
tensively in serial computers for retiming 
signals. One terminal of the circuit is fed 
by the signal to be retimed while the other 
is fed by digit pulses from a master clock. 
There will be an output from the circuit 
only for the overlap of the master digit 
pulse and input signal. The circuit 
whereby the output is made a replica of 
the pulse from the clock is discussed in 
connection with the basic amplifier. 

An inhibitor terminal can be added to 
any and-circuit or or-circuit. Such a cir¬ 
cuit operates as though there were no in¬ 
hibitor terminal when the inhibiting pulse 
is absent. When the inhibiting pulse is 
present, however, the circuit prevents any 
output from being developed. The in¬ 
hibiting circuits used in the designs 


proposed herein are of the simple variety 
shown in Figure 3 where positive in¬ 
puts synchronized in time are required. 
Note that the signal to be inhibited is 
passed through an eighth-digit delay 
line, while the inhibiting pulse is passed 
both through and around a quarter¬ 
digit delay line. This insures that the 
inhibitor pulse will, in effect, arrive 
earlier than the signal pulse and last 
longer. In the absence of input pulses, 
crystal diode X A will clamp the output at 
ground because input B is returned to a 
negative potential. Note that Xi and X 2 
are returned through the transformer to a 
positive potential. If input B goes posi¬ 
tive (without an inhibiting pulse appear¬ 
ing on input A) Xi will be cut off and the 
output voltage will rise until it is clamped 
at the positive potential to which X\ and 
X 2 are returned. If there is an inhibiting 
pulse (positive) it is inverted by the trans¬ 
former and will carry X\ and X 2 negative, 
which will keep Xi conducting, no matter 
what happens at B. Thus, if pulses A and 
B were written as a 2-digit binary number 
AB, the circuit translates 01 into a 1 at 
the output. It translates 00, 10, and 11 
into zero at the output. 

Active elements will be used not as 
flip-flops or switches but as repeating 
amplifiers to make up for attenuation in 
crystal diode circuits and delay lines. 
The standard use will include a retiming 
feature as well as amplification. Wherever 
a pulse is likely to suffer intolerable at¬ 
tenuation, deformation, or a variable de¬ 
lay, a circuit like that of Figure 4 is in¬ 
serted in the machine. 

The assembly shown has two inputs, A 
and B. Input A is the pulse to be re¬ 
timed and amplified. Input B comes from 
the master clock. This component sup¬ 
plies reference pulses (known as digit 
pulses) every digit time. These pulses are 
available in various phases, that is, with 
various but accurately controlled delays 
of a fraction of a digit time. The pulse 
fed to B is selected to rise sometime be- 



A = PULSE TO BE RETIMED 

B=DIGIT PULSES-REFERENCE 

PULSES FROM MASTER CLOCK 


Figure 4. Amplifier with pulse retiming 

If there is no input on A, there will be 
no output from the amplifier because of 
the and-circuit. If there is an input on A, 
when the digit pulse arrives the amplifier 
output will rise with a rise time deter¬ 
mined by the digit pulse (assuming that 
the amplifier passband does not limit it). 
Part of the amplifier output is fed back 
through an or-circuit to the and-circuit. 
This insures that the output pulse will not 
fall until the reference digit pulse does, 
even though pulse A may have ended 
after B rose. 

This reshaping with logic circuits and 
an amplifier is the way in which every 
pulse is maintained with the desired time 
synchronization. Pulse A may vary 
somewhat in the delay it has suffered but 
the output pulse still will leave the am¬ 
plifier at a time determined only by the 
reference pulse from the master clock. 
Thus, the pulses in the computer are made 
to have fixed durations and to occur at 
designated times. 

In the tabulations of components in the 
following paragraphs an active element is 
intended to represent as one unit a 
vacuum tube with transformer or a tran¬ 
sistor. The term will include neither the 
and-circuit nor the or-circuit shown in 
Figure 4. 

The basic storage cell proposed is not a 
static device like a flip-flop but is an 
electric delay line plus an amplifier. When 
vacuum tubes are used, this type of stor¬ 
age saves one active element in a 1-digit 
storage cell and is believed to be a more 
reliable use of active elements. In larger 
storage units more elements will be 
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Tabic I. Basic Building Blocks 


Unit 

Crystal 

Diodes 

Digits of 
Delay ! 

Active 

Elements 

N-terminal or -circuit. . 

.. ,N_ 

.... 0 . . . 

. . ,0 

N-terminal and -circuit. 

.N + l.. 

.. . 0 . . . 

. . .0 

Inhibitor-circuit. 

. . . 4 . . 

. . . 8 4, . 

0 

N-digit storage cell. 

. . .10*.. 

....N... 

.. .1* 


* One active element will be used for IV <8; for 
jV S:8 the number of active elements planned is the 
smallest integer 3: A r /8. Six additional crystal 
diodes will be required with each active element. 


delay line before it is amplified and recir¬ 
culated. The erase signal is received 
whenever new data are to be stored. It 
serves to erase the data in storage, block¬ 
ing the delay line output from its input 
until the new data have been inserted. 

The delay line may be long enough to 
store one word or just one digit of data. 
It is believed that up to 15-digit delay 
lines with lumped impedances can be 
built that will hold the delay constant to 
within a small fraction of one digit time. 
Depending on the length of a word, it may 
be necessary to break 1-word lines into 
sections and insert an amplifier between 
sections to retime and regenerate the 
pulses stored. To insure conservative es¬ 
timates, the estimates made in following 
sections are on the basis of regeneration 
after every eight digits of delay. 

The circuits which have been discussed 
use the components listed in Table I. 
This table will be used frequently in es¬ 
timating the parts required in the larger 
assemblies. 

Switches 

The switches are planned to combine a 
switching and a storage function, When 
a switch is given instructions to go to the 
Jfe’th position, it goes there and it re¬ 
members that it is to remain there (self¬ 
locking operation). All the elements of 
the switches have been discussed in pre¬ 
ceding paragraphs. How these elements 
are combined to make switches is de¬ 
scribed in the following: 

A single-pole double-throw switch, as 
shown in Figure 6, consists of a storage 
cell and a switch unit. When a 1 is stored 
in the storage unit, as the result of a pulse 
on the switching instruction lead, the left- 
hand and-circuit of the switch unit will 
pass signal a while the inhibitor blocks 
signal b. When a zero is stored, the and- 
circuit will block signal a and the in¬ 
hibitor-circuit in the switch unit will pass 
signal b. 

Whenever the switch is to be reset, an 
erase signal is fed to the storage unit, 
which then, drops its old instruction and 


goes to position b unless the new instruc¬ 
tion sets it to a. 

The switch unit itself uses 13 crystal 
diodes, 3/8 digit of delay, and one active 
element, as can be ascertained from Figure 
6 and Table I. A complete single-pole 
double-throw switch has 23 diodes, l 3 / 4 
digits of delay, and two active elements. 
A 2-pole double-throw switch would have 
two switch units and one storage unit. 
A 3- or 4-pole switch would have three or 
four switch units, one storage unit, and 
perhaps an extra amplifier to prevent the 
switch units from loading down the stor¬ 
age unit excessively. 

Switches with more than two positions 
can be assembled in a slightly different 
manner from the double-throw switches. 
Suppose, for discussion, that an 8-position 
switch is needed. A 3-digit code must be 
sent to the switch to specify the position 
it is to select. A convenient way to oper¬ 
ate such a switch is to translate the 3- 
digit code into a 6-wire code so that each 
digit is represented by signals of opposite 
phase on a pair of wires. In this system, 
a 1 is represented by a positive pulse on 
the positive bus and a negative pulse on 
the negative bus. A zero is represented 
by a negative pulse on the positive bus 
and a positive pulse on the negative bus. 
A 3-terminal and-circuit is provided for 
each of the eight lines that may be se¬ 
lected. The and-circuits are operated by 
positive pulses, and all are connected to 
either one or the other of each pair of wires 
representing the three digits of the posi¬ 
tion code. The and-circuit on the fifth 
(101) wire, for example, is connected to 
the positive bus of the pair representing 
the coefficient of 2 2 , to the negative bus of 
2 1 , and the positive bus of 2°. The three 
leads of the and-circuit will be energized 
by positive pulses only when 101 is sent to 
the switch as instructions. 

A single-pole 8-position switch is shown 
in Figure 7. The eight 3-terminal and- 
circuits (one for each horizontal lead), plus 
the 9-terminal or-circuit make an easily 
studied crystal matrix which feeds an¬ 
other and-circuit for retiming of pulses and 
an amplifier on the output. The and-cir¬ 
cuits are fed by switches So, Si, and Si, 
which are driven by the switching in¬ 
structions. These three switches are 
similar to the one shown in Figure 6. Ao is 
operated by the coefficient of 2°, Si by the 
coefficient of 2 X , and Si by the coefficient 
of 2 2 . As the figure is drawn, only the 
and-circuit on line 7 (111) is energized. A 
pulse will be received at the output when¬ 
ever a pulse is put on line 7. Pulses can be 
put on any of the other lines without get¬ 
ting to the output, because every other 
line is held negative by at least one lead 


Tabic II. Switches 


Type of 

Switch 

Crystal 

Diodes 

Digits of 
Delay 

Active 

Elements 

1P.2T . 

. . 23. 

. .i«A. 

... 2 

2P.2T . 

. . 36. 

. .ltd. 

... 3 

4P.2T . 

. . 63. 

. .6'/!. 

... 6 

1P.4T . 

. . 61. 

. .23/ 4 . 

... 4 

1P.8T . 

..105. 

. .6 . 

... 7 

1P.16T . 

..160. 

. .7 . 

... 9 

1P.32T . 

..311. 

■ .83/4. 

... 11 

Where 2P.2T, 

for example 

, represents 

a 2-pole 


2-throw switch. 

from S 0 , Si, or Si. If an erase signal were 
sent to the three storage cells of S 0 , Si, and 
S 2 , the number 111 would be erased, leav¬ 
ing 000, and the switches would be free to 
move to new positions and select another 
one of the eight inputs. 

Based on this type of mechanization, an 
r-position switch requires n single-pole 
double-throw switches, where n is the 
smallest integer equal to or greater than 
log 2 r. The switch requires r (n-\- 1) 
crystal diodes, plus four diodes for retim¬ 
ing the output, and an output ampli¬ 
fier. 

Table II gives estimates of the crystal 
diodes, digits of delay, and active ele¬ 
ments required in switches of different 
degrees of complexity. 

Handling of Negative Numbers 

The most significant digit place of 
every number is reserved to indicate the 
sign of the number. Positive numbers 
have a zero in the last place. A negative 
number is obtained by taking the 2’s com¬ 
plement of the positive number. This 
results in every negative number having 
a 1 in its last place. The system is equiv¬ 
alent to the 10’s complement method 
used in decimal calculators. In a decimal 
calculator operating with three signifi¬ 
cant figures, a fourth place might be pro¬ 
vided for the sign. The number —187 
might be represented by its 10’s comple¬ 
ment 9,813. Then, for example, if —187 
were required to be added to 500 the 
operation would be to add 9,813 to 500, 
which gives 10,313, and thus is recognized 
as 0,313 since the machine was assumed to 
have only four digit places. 


Tabic III. Binary Addition 


Addend 

Inputs 

Augend 

Carry 

Outputs 

Sum New Carry 

0. . . . 

_0_ 

. .. .0. 

....0. 

.0 

0. .. . 

....0_ 

...1) 



0. . . . 

_1_ 

...of.... 

.1. 

.0 

1_ 

....0_ 

...1) 



0. ... 

_1_ 

.. .1) 



1. .. . 

_1_ 

...of... 

.0. 

.1 

1. .. . 

....0. .. . 

...11 



1. ... 

_1_ 

_1. 

....1. 

.1 
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SWITCH UNIT 



Figure 6 (left). 
A simple switch 


O SETS SWITCH TO 'V' 
1 SETS SWITCH TO "a" 


SYMBOL 

- 

-°b 

MECHANIZATION 


Figure 8 (right). 
Handling of neg¬ 
ative numbers 


fal FORMATION OF 

NEGATIVE NUMBERS 

.+ 8 = 

0001000 

ONES 

COMPLIMENT = 

1 I10 1M 

ADD ONE 

1 

.•.-8 = 

1111000 

CHECK -8 + 8 = 0 

-+8 = 

0001000 

-8 

1111000 

SUM = 1 

0 0 0 0 000 = ZERO 

TO THE MACHINE 


(b) 

ADDITION 



-8 +5 = -3 
-8 = 

+5 = 

1 i 1 1 0 0 0 

0 0 0 0 VO 1 


SUM = 

CHECK +3 = 

1 1 1 1 1 0 1 = -3 
0000011 


SUM = 1 

0000000 = ZERO 


-8 + 15=7 
-8 
+ 15 

1111000 
0001111 


SUM = 1 

0000111=7 


NEGATIVE NUMBERS IN A 6 
SIGNIFICANT DIGIT SYSTEM 
WITH ONE IN SEVENTH PLACE 
SIGNIFYING THAT THE NUMBER 
IS NEGATIVE 


(C) MULTIPLICATION 


-8x7 =-56 
-8= I 1 I t 000 
7 = 0 0 0 0 111 



1 1 1 10 0 0 


0 0 0 0 1 1 1 


1111000 

1 

1 1 10 0 0 

1 1 

110 0 0 

000 

0000 

0000 

0 0 0 

00000 

0 0 

0 00000 

0 


10 0 10 0 0 


f SHOWS ANSWER 
IS NEGATIVE 

CHECK 
-56+56 =0 
-56 100 1 00 0 

+ 56 1 1 10 0 0 

ito 0 0 0 0 0 0 


(d) ANOTHER MULTIPLICATION 


-8X-7 

= + 56 

-8 = 

11 1 1 1 00 0 

-7 = 

1111001 


1111000 

0 

000000 

0 0 

0 0 000 

1 1 1 

10 00 

Mil 

00 0 

1110 

0 0 

11100 

0 


0111000 


INDICATES 
- POSITIVE 


A negative number (2’s complement) 
can be obtained in the binary computer 
by first, forming the l’s complement 
(changing all zeroes to l’s and vice versa 
by means of an inhibitor circuit), and 
then adding 1. Figure 8 shows several ex¬ 
amples of binary arithmetic performed 
with negative numbers. 

A point of interest is that even though 
either or both the multiplicand and multi¬ 
plier are negative the correct sign will be 
obtained for a product provided the nega¬ 
tive numbers are increased from their 
normal length W to 2W-1 by filling in l’s 
before multiplying. This is necessary be- 

SYMBOL 


cause the product of two IF-digit num¬ 
bers, where the last digit specifies sign, is 
a number 2IF-1 digits long. Unless the 
multiplicand and multiplier are increased 
to length 2W-1 (when negative) the 2IF-1 
place may be incorrect. This lengthening 
is not required for positive numbers, be¬ 
cause the digits from W to 2W-1 would be 
zeroes if they were filled in and would 
contribute nothing to the product. A 
consequence is that the product XY can 
be obtained to 2IF-1 places with W{W-l) 


NUMBER 


elementary additions provided X and Y 
are positive while W(2W-1) additions will 
be required if both are negative. Where 
fast multiplication is desired it may be ad¬ 
visable to convert negative numbers to 
positive ones, multiply, and then adjust 
the sign of the product. 

Adder 

The adder can be considered as a trans¬ 
lator with three inputs: addend, augend, 
and carry. It is a simple translator in that 
its output is a function only of the number 
of ones among its three inputs, as can be 


■O 



seen from Table III. 

The combination 0 0 0 is automatically 
taken care of in the adder shown in Figure 
9. The three dashed circuits at the left of 
the block diagram recognize three situa¬ 
tions among the three inputs. The situa¬ 
tions are: at least one 1, at least two 
l’s and three l’s among the inputs. 11 
there is only one 1, it will go through the 
bottom or-circuit, the following inhibitor 
circuit, and then another or-circuit. After 
being recloclced and amplified it will pro? 
vide a 1 as the sum. In this case none oj 
the and-circuits on the A, B, and earn 
leads will have operated. If there are at 
least two 1’s on the A, B, and carry leads 
at least one of the three 2-terminal and- 
circuits in the dashed box will operate 
with two results. The output of the 3; 
terminal or-circuit at the bottom left o 
the diagram will be inhibited so that if 
makes no contribution to the sum. Ii 
addition, a carry signal will be developer 
which is delayed one digit, reclocked, an< 
amplified to serve as the carry for th 
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Figure 9 (left). 
Adder 


Figure 11 (right). 
Multiplier 


WORD 

PULSE 


DIGIT 

PULSES 


next augend and addend. If there are 
three l’s, the 3-terminal and-circuit at the 
top and left-hand side of the diagram will 
operate and develop a sum of one. The 
three 2-terminal and-circuits on the left 
of the diagram also will have operated and 
provided the carry. Thus, the adder 
table is mechanized with two active ele¬ 
ments (amplifiers). 

The inhibitor circuit in series with the 
carry lead should be noted. This circuit 
is fed by a word pulse as well as the carry 
digits. The word pulse is received in 
synchronism with the first digit of every 
number. The word pulse will inhibit the 
carry pulse if one is present and will pre¬ 
vent a carry developed in one problem 
from being used in the next. This feature 
is required in the addition of negative 
numbers. 

The carry lead is brought outside the 
adder to facilitate subtraction. Suppose x 
is to be subtracted from y. The number y 
might be fed to the addend terminal and 
the number x fed through an inhibitor 
circuit to the augend terminal. The in¬ 
hibitor would also be fed by digit pulses 
from the master clock, and the augend 
would therefore be the l’s complement of 
x. A 1 would be inserted into the carry 
terminal in synchronism with the first 


DELAY LINE 



digit of y and of the l’s complement of x. 
The sum out of the adder would then be 
y — x. The adder produces a sum within a 
fraction of a digit time after it receives an 
input. Thus, there is only a small delay in 
obtaining the sum of two numbers. 

Table I can be consulted to show that 
the adder of Figure 9 requires 38 crystal 
diodes, l 3 / 8 digits of delay, and two active 
elements. The operation of so many 
crystal diode circuits in series without am¬ 
plifiers may be questioned, but it is be¬ 
lieved that it will become feasible as 
diodes and transistors are improved. 

Accumulator 

A block diagram of an accumulator is 
shown in Figure 10. The output of the 
adder is fed back to its input through a w- 
digit delay line. The output of the delay 
line is reclocked continuously in the and- 
circuit in accord with digit pulses from 
the master clock. The timed pulses then 
are amplified in a 1-stage amplifier to 
make up for attenuation in the delay line. 
■Whenever a new accumulation is to be 
started, an erase signal is fed to the in¬ 
hibiting circuit. This signal is w-digits 
long and blocks the output of the delay 
line from the adder and insures that the 
new accumulation will start from zero. An 
accumulator to accumulate sums 48 digits 


Figure 10 (left). 
Accumulator 



long will require 76 crystals, 48 digits of 
delay, and eight active elements. 

Multiplier 

The multiplier that is discussed in the 
following paragraphs is designed to mul- 
. tiply two positive members, each having 
W digits, in W(2W+l) digit times. Oper¬ 
ating at a megacycle rate, the product of 
two 48-binary digit numbers would be ob¬ 
tained in 4,656 microseconds. It is be¬ 
lieved most efficient to convert all nega¬ 
tive members to positive ones before mul¬ 
tiplying and to adjust the sign of the 
product according to the rules of alge¬ 
braic multiplication. The components 
to be added to the multiplier to make the 
sign correction have not been included in 
this study. The multiplicand x is multi¬ 
plied by the first (least significant) digit 
of the multiplier y in the first (2FF+1) 
period and by the W’th digit of y during 
the W’th (2PF+1) period. The W partial 
products are added together, with each 
partial product moved over one place with 
respect to the preceding one as it is ac¬ 
cumulated. 

The multiplication table in binary 
arithmetic is very simple. If it is desired 
to multiply the binary number x by a 
single binary digit m, x is connected to the 
input of a switch that is held open if 
m = 0, and closed if m — 1. The output of 
the switch will be mx. Thus, in multiply¬ 
ing x by the successive digits of the num¬ 
ber y, it is only required to pass x through 
a switch that is opened and closed at the 
appropriate times according to whether 
the successive digits of y are zeroes or 
l’s. 
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DIVIDE 15 (lltl) BY 24 (llOOO) 
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| 1 1 1 1.000 
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11 000 
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Figure 1 3 (above). Binary division 
Figure 15 (right). Divider 

The block diagram of a multiplier is 
shown in Figure 11. The multiplicand x 
is stored in delay line 1, and the multiplier 
y is stored in delay line 2. After the first 
word period, 51— A, Sl—B, and 51 — C 
are thrown to position b and the multi¬ 
plier, multiplicand, and start pulse re¬ 
circulate in their delay lines until the 
multiplication is finished. DLl (delay 
line 1) storing the multiplicand has one 
more digit of delay than line DL2 in 
which the multiplier is stored. Therefore, 
after the first circulation of x through 
DLl, the least significant digit of x leaves 
DLl just as the second least significant 
digit of y leaves DL2. After the W— l’th 
circulation of x, its least significant digit 
leaves DLl just as the most significant 


DIVISOR 


1 HERE SETS 
S2 TO"b" 
i HERE SETS 
S2 TO "a" 


START 

PULSE 



QUOTIENT 


digits of y, choosing a new digit every 
2w-\-l digit times. The examining pulse 
is also sent to 52, where it is the erase 
signal of the storage cell associated with 
52. The output of the and-circuit goes to 
52 as switching instruction and is stored 
there. Thus, 52 stores either a 1 or a 
zero, 'according to whether the digit of y 
examined is a 1 or a zero. The switch is 


a one is written in the quotient and 10.00 
is subtracted from 11.00. The remainder 
is 1.00. The binary point is again shifted 
in the divisor, giving 1.000 which goes ex¬ 
actly once, and another one is written in 
the quotient giving an answer of 0.011 or 
0 X V 2 h 1 X V 4 h 1 X Vs equal to 3/S. 

It should be noted that at each step the 
divisor is either subtracted from the re¬ 


digit of y leaves DL2. 

The output of DL2 is fed continuously 
to an and-circuit. Every 2 W- f-1 digits, an 
examining pulse is fed to the and-circuit 
from DL3, and if the digit of y coming out 
of DL2 is a 1 at that time, the and-circuit 
develops a 1 in its output. If, when the 
examining pulse comes along, the coef¬ 
ficient of y leaving DL2 is a zero, the and- 
circuit will have no output. Because 
DL2 has a delay of 2w digits, the and- 
circuit in effect filters out the successive 


DIVIDE 3 BY 8 

WORD LENGTH =5 PLACES. 5TH PLACE INDICATES SIGN. 
3 = 0 0011 

8 = 0 1000 

ONES COMPLEMENT = 10 111 

ADD ONE _1_ 

-8 = 11000 

CHANGE NUMBERS TO 2W-1 OR 9 DIGIT LENGTH. 

3 = 0000000 1 1 
8= 000001000 
“8 = 1111110 0 0 

DIVISION RULE: IF (2W-l)TH ( 9 th) DIGIT OF REMAINDER 
IS ZERO ADD-DIVISOR; ADD + DIVISOR IF DIGIT IS ONE. 

_ 0.0 1 1 0 0 

0 00001000 | 0000000 11.00000 
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NOTE 9TH DIGIT—>111 1 I 1 I 0 1 I 
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+ B WITH SHIFT 10 0000 1000 

NOTE 9TH DIGIT->l|l 1111110 

AND ADD 

+ 8 WITH SHIFT _ 0 0 00 01 000 

NOTE 9TH DIGIT-volo 00 00100 

AND ADD 

-8 WITH SHIFT _ [i ll 11 10 0 0 

NOTE 9TH DIGIT-»-0|00 00 000 0 

AND ADD 

-8 WITH SHIFT _ 11 1 11110 00 

NOTE 9TH DIGIT - >l|l 11 1 10 00 

AND ADD 
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NOTE 9TH DIGIT-i-tli 1 1 11 00 0 


THE QUOTIENT IS THE ONES COMPLEMENT OF THE 
SUCCESSIVE 9TH DIGITS OF THE REMAINDERS. 


Figure 14. Simplified division process 


arranged so that when a 1 is stored, 52 is 
closed, and when a zero is stored, it is 
open. The examining pulse is developed 
from the start pulse which is received by 
the multiplier coincident with the first 
digits of x and y and recirculated in delay 
line 3. 

A new partial product is obtained at the 
output of 52 every 2wf-l digit times by 
the process just described. To accumu¬ 
late them, an adder is provided which ac¬ 
cepts as one input the output of 52, and as 
the other input its own sum output de¬ 
layed by 2w digits. The one digit place dif¬ 
ference in the arrival of the two inputs of 
the adder automatically provides the 
“shift” feature necessary in adding the 
partial products. Figure 12 shows a digit 
by digit account of the operation of a 3- 
digit multiplier when multiplying 110 by 
101. 

A multiplier based on Figure 12 can be 
mechanized for 48 digit words with 390 
crystals, approximately 390 digits of de¬ 
lay, and 55 active elements. 

Divider 

The divider operates with repeated 
subtractions and shift operations much as 
one would do in pencil and paper division. 
Figure 13 shows two division examples. 
In the first one, 3 is divided by 8. First, it 
is seen that 1,000 will not go into 11, so 
the first zero is written in the quotient. 
100.0 will not go into 11, so the second 
zero is written. 10.00 will go into 11.0, so 


mainder or nothing is subtracted and that 
a shift to the right is made each time be¬ 
fore the divisor is subtracted. In effect, 
the successive subtrahends are obtained 
by repeated divisions of the divisor by 2. 
For machine operation, the awkward part 
of this process is that it cannot be ascer¬ 
tained whether a subtraction should be 
made without making it and seeing 
whether or not the remainder is positive 
or negative. If the remainder is nega¬ 
tive, the subtrahend should be added 
back. This necessity for sometimes add¬ 
ing back the subtrahend may lie avoided. 
Let x be the subtrahend and R the minu¬ 
end. Suppose x is subtracted from R and 
the remainder found to be negative. Dur¬ 
ing the next step the machine should sub¬ 
tract x/2 from R, but the machine has 
lost R by subtracting x from it. R can be 
regained and x/2 subtracted at the same 
time by adding the factor x/2 to the nega¬ 
tive remainder R — x which gives the de¬ 
sired result R—x/ 2. Thus, on the succes¬ 
sive steps the machine either subtracts 
half the previous subtrahend or adds half 
the previous subtrahend, depending on 
whether the previous remainder was posi¬ 
tive or negative. The process is simpli¬ 
fied further by avoiding subtraction 
through adding the 2’s complement of the 
subtrahend rather than subtracting the 
subtrahend. 

An example of the simplified division 
process is shown in Figure 14. The ma¬ 
chine is assumed to have five digit places. 
The first step is to fill out the divisor and 
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Figure 1 8 (right). 
Shift register 


REGISTER FOR PARALLEL READ OUT 


dividend to 2w—l places, and then to 
form the complement of the divisor. The 
divisor complement is first added to the 
dividend. The 2w— l’th digit of the sum 
is noted. In the example, the digit is a 1, 
indicating a negative result. The ma¬ 
chine then writes a zero in the quotient, 
drops out the 2w— l’th digit of the previ¬ 
ous sum, and adds half the divisor (divisor 
shifted one place to the right). The 
2w— l’th digit is examined again and is 
found to be a 1. Another zero is written 
in the quotient, the 2w— 1’th digit is 
dropped out, and the divisor is again 
shifted and added. This time the 2w— 
l’th digit is zero, indicating a positive 
sum. A 1 is written in the quotient, and 
the complement of the divisor is shifted 
and added. This process is continued 
until the 5-digit answer 0.0110 is ob¬ 
tained. Note that the first operation was 
made with the binary points of the divisor 
and dividend lined up. In general, it 
would be necessary to start with the least 
significant digit of the divisor written 
under the most significant digit of the 
dividend. 2iv(2w— 1) digit times are re¬ 
quired for division in a general-purpose 
machine. 

The divider block diagram is shown in 
Figure 15. The divisor is stored in delay 
line 1 and recirculates there, becoming 
available every 2w— 1 digit times. The 
1’s complement of the divisor is taken by 


the inhibitor circuit fed by DLl (delay 
line 1). Eveiy 2w—\ digit times switch 
S2 selects either the l’s complement of 
the divisor or the divisor itself to go into 
the adder. The switch selects the l’s 
complement when it is desired to sub¬ 
tract. However, to obtain subtraction it 
is necessary to add the 2’s complement. 
This is done by bringing out the carry ter¬ 
minal of the adder, and whenever S2 is 
thrown to select the l’s complement a 1 is 
sent into the carry terminal and, in ef¬ 
fect, makes the augend the 2’s comple¬ 
ment. 

The storage cell of switch S2 receives an 
erase signal every 2w~l digit times. At 
the same time, the switch receives a 1 or a 
zero from the inhibitor circuit on the 
adder output. If this inhibitor sends a 1 to 
the switch, the 2iv~ 1’th digit of the sum 
just completed was a zero indicating a 
positive number and switch 52 goes to 
position b to select the 1 ’s complement of 
the divisor for the next addition. The 
output of the inhibitor circuit also is sent 
to the carry terminal of the adder, so that 
when 52 is at b a 1 will be injected to con¬ 
vert the l’s complement to a 2’s comple¬ 
ment. Each output of the inhibitor is 
also a digit of the quotient, for if the 
2w— 1’th digit of the sum is a zero, a 1 
should be written in the quotient and the 
inhibitor supplies that 1. It is written in 
delay line 3, which has a circulation time 



Figure 17. Binary counter 


of 2w. Since a new quotient digit is re¬ 
ceived every 2w— l’th digit time, delay 
line 3 has as its output the digits of the 
quotient in reverse order to that in which 
they are obtained, which puts the least 
significant digits first as they should 
be. 

The foregoing discussion shows how the 
augend is obtained for the adder. The ad¬ 
dend is obtained from SIB. This switch 
selects the dividend during the first 2w— 1 
period, and thereafter selects the previous 
sum from delay line 2. Note that an in¬ 
hibitor circuit between the adder and 
DL2 blocks out the 2w— l’th digit each 
time to keep the remainder from growing 
longer than 2w— 1. Since delay line 2 has 
2 w delay and delay line 1 has 2w — 1 de¬ 
lay, the augend arrives one digit early 
each 2w— 1 period with respect to the ad¬ 
dend. This feature supplies the succes¬ 
sive divisions by two that are required. 

Delay line 4 exists merely to supply the 
2w— 1 pulse that is sent to S2 and the in¬ 
hibitor circuit on the adder output. 
Switch 51 is driven to a at the end of the 
computation period to empty the delay 
lines of the dividend, divisor, and 2w~\. 
pulse preparatory to solving a new prob¬ 
lem. After the first 2w—l digit times of a 
new computation, all sections of 51. are 
thrown to b and recirculate their data. 

An estimate based on Tables I and II 
of the parts required by the divider of 
Figure 16 shows that to handle 48 digit 
numbers 371 crystal diodes, 390 digits of 
delay, and 51 active elements would be re¬ 
quired. This assumes that the delay lines 
have their data regenerated every eight 
digits. If the data are regenerated only 
every 16 digits, only 36 active elements 
are required. The divider is seen to be of 
about the same complexity of the mul¬ 
tiplier. 


July 1952 


Felker—Typical Block Diagrams for a Transistor Digital Computer 























Translator 

This unit translates one set of numbers 
or words into another set. Storage may 
be regarded as translation in which ad¬ 
dresses are translated into the words 
stored at the addresses. 

Translators are useful in developing 
switching orders. At the completion of 
some process, it may be necessary to de¬ 
velop a new set of orders. When the 
orders must be developed in a certain se¬ 
quence, a binary counter can be used to 
develop the input for the translator. 
Whenever the next set of orders are re¬ 
quired, a pulse is fed to the counter. The 
binary count number is used as the 
switching instruction for a switch. Every 
time the counter receives an input pulse 
the switch will be advanced one position 
and will put a pulse on the selected lead. 
The various output leads of the switch 
will be connected to a group of or-circuits. 
One or-circuit will be required for each 
digit of the translator output. Figure 16 
shows a scheme for generating eight sets 
of 3-digit words. To generate n digit 
words n or-circuits would be required. 
No more than eight leads would be re¬ 
quired for any or-circuit. The translator 
shown develops the output word in paral¬ 
lel. If the or-circuits are fed to taps in a 
delay line, the output word will be re¬ 
ceived serially from the output of the de¬ 
lay line. 

Binary Counter 

Figure 17 shows how each stage of a 
binary counter can be made from a 1-digit 
memory cell plus an inhibitor and an and- 
circuit. In Figure 17, if the count is at 01 
and a 1 is sent into the counter, the 1 in 
the first stage will be admitted to the 
second stage. The 1 that was in the first 
stage inhibits the 1 input for the first 


stage and the stage therefore goes to 0 and 
the count is 10. If another 1 is sent into 
the counter, it will enter the first stage and 
the count will advance to 11. If another 
1 is fed in, the first two stages will set to 
zero and the third will receive a 1 making 
the count 100. This type of counter can 
be used only to count synchronized input 
pulses but that is all that a counter is ex¬ 
pected to do in a serial computer. 

In a serial machine binary counters 
often can be replaced by simpler de¬ 
vices. Suppose that it is desired to time 
some operation and know when 64 digit 
times have elapsed. A 6-stage counter 
could be used or a digit pulse could be cir¬ 
culated in a 7- and a 9-digit delay line. 
If an and-circuit were fed by the two de- 
lay lines, 63 digit times would elapse be¬ 
fore the and-circuit were operated. If 
the output were fed through a 1-digit de¬ 
lay line the desired 64 digit delay would 
be obtained. 

Shift Registers 

Shift registers are sometimes needed to 
take data at a low rate, from a magnetic 
drum, for example, and retransmit it at a 
rapid rate. A shift register that is loaded 
serially may be read out in parallel on IF 
leads or in series on one lead. Each stage 
of such registers can be built around a 1- 
digit memory cell. In Figure 18 are 
illustrated two stages of each kind of 
register. The numbers in each cell are 
erased and advanced to the next cell 
whenever an advance pulse is received. 

Advance pulses will be received every 
time a new digit is to be loaded. When a 
number is in the register, the advance 
pulses can be stopped and the number 
can be stored until it is required. In the 
case of the parallel-read-out registers, ex¬ 
tra and-circuits are used at each stage. 
These are fed by the read-out signal. 


Besides their use in changing the rate of 
transmission of data, the shift registers 
described can store words for integral 
multiples of a digit time. Where the 
storage period is an integral multiple of a 
word time it is more efficient to use delay 
lines as discussed previously. It may be 
mentioned that in a synchronous com¬ 
puter tapped delay lines can replace shift 
registers for the conversion of serial to 
parallel data and vice versa. 


Conclusions 

A fairly extensive catalogue of digital 
computer components has been given in 
block diagram form with estimates of the 
number of parts they require. The de¬ 
signs have the common feature of using 
active elements as amplifiers only. The 
principal conclusion to be drawn is that an 
all-semiconductor computer can be built 
with diodes and transistors. This paper 
has been written as a contribution to¬ 
wards the development of such a com¬ 
puter. 
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Computing Circuits and Devices for 
Industrial Process Functions 

A. J. HORNFECK 

MEMBER AIEE 


P OTENTIOMETER and bridge-bal¬ 
ance systems using galvanometer de¬ 
tectors and electromechanical amplifiers 
have been used successfully in the re¬ 
cording of industrial process variables for 
many years. Electronic detection and 
amplification has greatly extended the 
field of application of these basic instru¬ 
ments by increasing'response speed, sen¬ 
sitivity, and range of application. 1 - 2 The 
versatility and standardization of design 
made possible by the electronic servo am¬ 
plifier has permitted the utilization of the 
accurate null-type metering system to the 
continuous measurement and control of 
practically all process variables which can 
be converted to electrical functions. The 
electronically operated null-balance re¬ 
corder is admirably suited to perform as 
an analogue calculator, and circuits and 
devices are available which permit the 
computation of derived process variables 
including sums, differences, products, 
ratios, powers, roots, and trigonometric 
functions. The electric computing sys¬ 
tems offer numerous advantages over 
equivalent mechanical methods. These 
include: 

1. Transmitters (also called transducers or 
pickups) whose outputs are to be included 
in the calculation can be located remotely 
from one another in their most favorable 
location relative to process. In this way 
the electric system offers the additional ad¬ 
vantage of shprt-range telemetering for 
centralized process control. 

2. Complete electronically operated com¬ 
puting systems are built up using the basic 
standard components of the simple measur¬ 
ing systems. Mechanical computers gener¬ 
ally require special design of linkages, 
cams, and gears. 

3. Functions can be computed accurately 
over wide ranges and with relatively 
straightforward calibration procedure. 
Mechanical computers involving complex 
linkages and cams are generally limited and 
require special calibrations for each applica¬ 
tion. 

Analogue-type computers using the 
servo operated systems are very well 
suited to the continuous recording of de¬ 
rived process functions. Digital comput¬ 
ers and also analogue systems which are 
entirely electronic in operation do not ap¬ 
pear, at least at the present, to be suitable 
for industrial process control. For this 


service the requirements are sturdiness of 
design, the necessity for remote indication 
and recording of results, long life, relia¬ 
bility, and ease of service. The funda¬ 
mental mathematical accuracy of the 
more complex computing machines usu¬ 
ally is not required since the accuracy al¬ 
most always is limited by the primary 
measuring element or the transducer 
which converts the physical quantity to 
an electrical function. 

Some of the applications of servo 
actuated computers of the automatic 
null-balance type are: 

1. Summation of fluid flows in chemical 
and other type industrial processes. 

2. Computation of ratios of related factors 
as, for example, fuel-air ratio in combustion 
control and control of the ratio of chemicals 
to obtain a desired solution in a process. 

3. Calculation of heat released or absorbed 
in heat exchangers and refrigeration sys¬ 
tems. 

4. Calculation of over-all plant efficiencies 
as, for example, that of a steam-electric 
generating station. 

5. Accurate compensation of measure¬ 
ments as, for example, automatic correction 
of gas flow meters for wide range variations 
in the density of the gas. 

6. Automatic and continuous solution of 
theoretical and empirical equations relating 
a number of process variables such as pres¬ 
sure, temperature, pH, conductivity, and 
flow to a final product output. 

The Elementary System 

Figure 1 is a diagram of the automatic 
null-balance system of measurement and 
shows in block form the basic components 
of the system. The wide application of 
the null-balance system to the measure¬ 
ment and control of process variables is 
illustrated. In some applications it is pos¬ 
sible to perform simple computations such 
as addition, subtraction, and division 
directly in the basic measuring circuit. 
For example, temperature difference can 
be computed by simply incorporating two 
temperature detectors in a bridge or po¬ 
tentiometer circuit in such a way that the 
balancing position of the follow-up is 
made directly proportional to the dif¬ 
ference. In other applications, particu¬ 
larly where the energy output of the 


measuring element is very small or where 
more complex calculations are required, 
retransmission from a primary instrument 
is employed. 

If the pickup is a primary measuring 
device such as a thermocouple, resistance 
thermometer, or photocell, the physical 
variable is converted directly into an 
electrical output such as current, voltage, 
or resistance. The miscellaneous or 
“quality” variables such as pH, vacuum, 
conductivity, and gas analysis are best 
measured by electric or electronic primary 
elements. Consequently, the output of 
these elements can sometimes be used 
directly in computers without the inter¬ 
mediate mechanical link. On the other 
hand, quantities such as pressure, fluid 
flow, and liquid level are converted most 
frequently into mechanical factors such as 
force or motion by means of a mechanical 
metering element or link which may be a 
float, diaphragm, bellows, liquid column, 
or bourdon tube. This metering element 
in turn actuates a position transmitter or 
pickup device whose electrical output is 
proportional to the physical variable. For 
man}' straightforward algebraic func¬ 
tions this output can be combined with 
the output of similar elements to perform 
the desired computation. However, in 
some applications the motion of the 
metering element bears a nonlinear rela¬ 
tion to the process variable, as for ex¬ 
ample: fluid flow is proportional to the 
square root of differential pressure drop 
across an orifice or flow nozzle. If such 
outputs are to be combined for totalizing 
and similar functions, the nonlinear 
quantity must be converted to an elec¬ 
trical factor proportional to the desired 
quantity. In these applications the func¬ 
tion is first extracted in a primary null 
balance type of receiver and retransmitted 
as an electrical quantity proportional to 
the process variable. 

The principles discussed are illustrated 
by a British-thermal-unit calculating sys¬ 
tem shown diagrammatically in Figure 2. 
It is desired to compute the rate at which 
heat is absorbed by the fluid flowing 
through the heat exchanger. Two tem¬ 
perature-sensing elements 7\ and jT 2 are 
placed in the inlet and outlet side of the 
heat exchanger and connected into an 
automatically balanced differential bridge 
or potentiometer circuit so that the mo- 
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tion of the motor driven follow-up is pro¬ 
portional to temperature difference. The 
balancing motor also actuates a retrans¬ 
mitter which may be a resistance slide- 
wire. Flow is measured by means of a 
mechanical-type element responsive to 
the differential pressure drop across an 
orifice or other type of primary element in 
the pipe. This pickup device converts 
differential pressure into a telemetered 
electric variable and the motion or posi¬ 
tion of the follow-up in the flow receiver 
is responsive to this differential. Since 
the flow is proportional to the square root 
of the differential, the measured function 
is converted to flow in the receiver by 
means of a properly shaped cam operating 
a retransmitter whose output is made 
proportional to the flow. The output of 
the flow retransmitter is then combined 
with the output-of the temperature-dif¬ 
ference receiver in a multiplying circuit 
in which the receiver moves proportional 
to the product of the retransmitted vari¬ 
ables to produce an indication or record¬ 
ing of the rate of British-thermal-unit ab¬ 
sorption. A more detailed description of a 
typical computing system is presented in 
.-a later section of this paper. 

Pickups and Their Circuits 

It is not the purpose of this paper to de¬ 
scribe in detail the wide variety of pri¬ 
mary metering elements, transmitting, and 
pickup devices as well as the basic meas¬ 
uring circuits of the null type frequently 
used in potentiometer and bridge-balance 
instruments. These are either well known 
or are adequately described elsewhere in 
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the literature. Primary measuring ele¬ 
ments may be grouped into two general 
categories; namely, the self-generating 
type in which a small amount of energy 
from the process is converted into an out¬ 
put voltage, and the transducer type in 
which the measured variable influences 
either in direct proportion or by some 
known relation the electrical character¬ 
istics of the pickup device. Examples of 
the first class are thermocouples, photo¬ 
cells of the barrier layer type, and pH 
electrodes. Examples of the transducer 
type are resistance temperature detectors, 
resistance wire strain gages, electrolytic 
conductivity cells, and heated filaments 
for gas analysis by thermal conductivity. 

Of more specific interest in the com¬ 
puting systems are the position pickup de¬ 
vices which convert position or motion 
into an electrical output. The following 
types of position transmitters and func¬ 
tion generators are used in the analogue 
computing systems. 

1. The adjustable core transformer, some¬ 
times called the differential transformer or 
linear induction potentiometer. 

2. The resistance slide-wire used either as a 
variable resistance or as a rheostat poten¬ 
tiometer. 

3. Induction potentiometers for converting 
angle of rotation to a proportional output 
voltage. 

4. Sine-cosine electrical resolvers for con¬ 
verting angle of rotation to electric voltages 
proportional to sine and cosine functions. 

The first two types of simple position 
pickups are used most extensively in the 
type of process function computers 
covered by this paper and will be de¬ 
scribed in detail. Both of these devices 
may be used either as transmitters 
directly coupled to mechanically operated 


metering devices or as retransmitters foi 
generating voltages proportional to e 
function of the original variable. The} 
are also used as servo-motor operated fol¬ 
low-ups in intermediate and final record¬ 
ing receivers. 

Adjustable Core Transformer 

Figure 3 illustrates the principles anci 
application of the adjustable core trans¬ 
former which is used extensively in the 
measurement of factors which are con¬ 
vertible to position or motion. The mag¬ 
netic core of the transmitter is moved by 
the mechanical metering element which 
in some cases actually may be a complete 
meter which indicates or records the vari¬ 
able and in others only a basic link be¬ 
tween the primary detecting element and 
the transmitting unit. Figures 3(A) and 
3(B) show the structure and principle ol 
operation. The magnetic flux linkages 
between the primary winding and each 
of the two secondaries is determined by 
the magnetic core position. Consequently 
the voltages induced in the secondary 
windings vary linearly and differentially 
in opposite directions in accordance with 
the following equations. 

Ei — Eo+kx 

Et-Eo-kx ( 1 ) 

where 

En — voltage when the core is in the exact 
center position 

v = per-unit motion of core = actual motion 
divided by total motion 
k = total change in voltage per winding at 
the extreme of core position 

In the ratio-balance circuit of Figure 
3(C) the secondaries are connected so 
that their voltages aid. At balance the 


Figure 1. Block diagram of the automatic, null-balance system for measurement, 

computation,"and control 


->- 

1 

d-c 

unbalance 

d-c to a-c 

C0IVEIT0I 



a-c 



VOLTACE AMPLIFIES 
A M0T0I C0IT10L 



a-c 

unbalance 


t 


Moilor Functions 



SLUE WIIE - 

SLUE WIIE 

ADJISTASLE CORE 
TIAISFOIMEI 

IEC0IIIIB PER 
IIIICATOI 
IITEIIATOI 

PILOT VALVE — 
CONTACTS - 


IETIARSMITTEI 


SLUE WIIE 

AIJISTAILE 00IE 
TIAISFOIMEI 


d-c Measurements 
a-c Measurements 


Pneumatic Control 
Electric Control 

Electric Control 
or Computing 


184 


Hornfeck—Computing Circuits and Devices for Industrial Process Functions 


July 1952 










Figure 2. Diagram of a 
typical process comput¬ 
ing system—British- 

thermal-unit calculator 


'TEMPERATURE SENSING ELEMENTS 


DIFFERENTIAL PRESSURE h = (FLOW)‘ 



ELECTRONIC 

AMPLIFIER 


MEASURES TEMPERATURE AND 
TRANSMITS PROPORTIONAL OUTPUT 


fundamental voltage input to the ampli¬ 
fier is 0 or a minimum, and the angular 
position of the follow-up slide-wire indi¬ 
cates the position of the transmitter core 
and consequently the measured variable. 
For this balance condition the following 
equation applies. 


Eo+kx B + MS 


Ev-kx A+(l -M)S 


( 2 ) 


If the values of A, B and S are properly 
chosen so that B = (E a S)/k and A = S 
(Eo/k — 1), the slide-wire rotation M will 
equal the core position x. In this analysis 
it is assumed that the magnetic core is in 
the center of the coil when the transmitted 
variable is zero. In many applications 


CONVERTS HEAD TO FLOW 
AND RETRANSMITS 
PROPORTIONAL OUTPUT 


BTU RECORDER 


the core is centered at 50 per cent in 
which case resistances A and B are equal 
or A=B = S[E 0 /{k-1/2)]. 

In the potentiometric circuit shown in 
Figure 3(D) the secondary windings of the 
transmitter are connected in phase op¬ 
position so that the total voltage is a pro- 



AMPLIFIER 


BTU PER HR * FLOW XAT x CONSTANT 

portional function of the core position. 
The identical transformer unit in the re¬ 
ceiver is positioned by the motor and the 
system is balanced when the transmitter 
and the receiver core positions are ex- 
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Figure 3. Basic meas¬ 
uring circuits using ad¬ 
justable core trans¬ 
former-type transmitters 
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Figure 4. Pickups of the adjustable core transformer type for specific applications 


actly the same. This system has several 
advantages over the one shown in Figure 
3(C), the most important being the fact 
that no current flows in the leads of the 
measuring circuit. Consequently, it can 
be used over greater distances without 
being affected by line resistance. 

In the slide-wire balance type of poten¬ 
tiometer system shown in Figure 3(E), 
two pairs of secondary windings are used 
on the transmitter. One pair with the 
secondary windings connected opposing 
sets up a voltage proportional to the 
transmitter core position. The other pair 
of windings connected with the voltages 
aiding sets up a reference potential for the 
balancing slide-wire. The system is cali¬ 
brated by adjusting the voltage across this 
slide-wire by means of the series resistance 
so that Eq is equal to the maximum output 
of the transmitter. 

A number of typical transmitters using 
the transformer pickup device are illus¬ 
trated by Figures 4(A) to 4(E) inclusive. 
Figures 4(A), 4(B), and 4(C) show de¬ 
vices of the integral construction type in 
which the transformer unit is designed as 
a basic part of the mechanical measuring 
device. This type of pickup is particu- 
larly adaptable to the telemetering of 
flow and differential pressures where the 
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moving element of the basic meter device 
including the transformer core can be en¬ 
closed entirely in a high-pressure 
chamber. 

In Figures 4(D) and 4(E), the trans¬ 
former is incorporated as a retransmitting 
element in recording meters. Figure 
4(D) shows the transformer transmitter 
operated through suitable linkage from 
the spindle of a mechanical flow meter 
and it is used in preference to the integral 
devices previously described where the 
core stroke is to be adjustable for calibra¬ 
tion purposes. Figure 4(E) shows the 
core of the transformer actuated by the 
motor-driven follow-up of an electrically 
balanced recorder. This system is par¬ 
ticularly advantageous where the out¬ 
put of the retransmitter must be a non¬ 
linear function of the measured variable. 
In a complex process computing system, 
many of these retransmitter converters 
might be used to order to translate the 
various metered functions such as pres¬ 
sure, temperature, and flow into elec¬ 
trical quantities. 

Resistance Slide-wire Circuits 
Resistance slide-wires are used either 


devices or as retransmitters in electronic 
and mechanical type meters, similarly to 
the adjustable core transformer already 
described. Since the transmitter may be 
situated in an unattended location, the 
slide-wire lacks the inherent sturdiness of 
the transformer device. Nevertheless, it 
is cjuite useful particularly where the 
computing system requires a resistance 
function. The friction of the sliding con¬ 
tact as compared with the essentially 
frictionless “floating action” of the mov¬ 
able core unit is also a disadvantage, par¬ 
ticularly if the slide-wire is actuated by a 
low-power device such as a Bourdon tube 
pressure element. 

Figure 5 illustrates the most commonly 
used null metering circuits of the bridge 
and potentiometric type which are used 
with the resistance slide-wire transmitter. 
These circuits are basically very simple 
and well known so that no detailed analy¬ 
sis of their characteristics is required. 
For algebraic calculations of process 
variable combinations the circuit shown in 
Figure 5(D) is particularly useful. In 
systems using this basic a-c potentiometer 
or ratio principle, the slide-wires function 
as sources of voltage proportional to the 
transmitted variables. The advantages 
of the potential balance circuits over 
Wheatstone bridge circuits perhaps are 
fairly obvious, but the most important 
ones are: 

1. Elimination of contact resistance errors. 

2. Relative ease of calibration and adjust¬ 
ment to take care of wide ranges in mag¬ 
nitudes. 

Figure 6(A) shows the low-torque slide- 
wire transmitter for operation directly 
from a mechanical metering element. 
Because of its low friction, it can be 
geared or linked directly to the output of a 
sensitive metering element for converting 
the physical variable to an electrical out¬ 
put. Figure 6(B) shows the internal 
mounting of this element in a recording 
mechanical meter. If no record is re¬ 
quired at the transmitting end, a pickup 
of the type shown in Figure 6(C) can be 
used in which the low-torque slide-wire is 
incorporated as an integral part of the 
mechanical device. Figure 6(D) shows 
the motor-driven slide-wire unit in the 
electronically balanced recorder. One 
slide-wire of this assembly is used to 
balance the primary metering circuit and 
the second slide-wire is used as a retrans¬ 
mitter for computing purposes. In some 
applications the adjustable core trans¬ 
former is also used as a part of this as¬ 
sembly either as a follow-up in a basic 
metering circuit or as a retransmitter as 
shown in Figure 5(E). 
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Basic Computing Circuits 

Transformer Type 


Figure 7 shows the elementary circuits 
using the transformer transmitter for the 
simple arithmetic functions of addition, 
subtraction, and division of transmitted 
variables. 

In the circuit of Figure 7(A) the per- 
unit rotation of the motor-driven balanc¬ 
ing slide-wire is equal to the total or 
average motion of all the transmitter 
cores. The lower or zero end windings 
all are connected in series aiding so that 
the total voltage of one leg of the bridge 
is given by 

£2 = Aiu+-Ei'2+-Ez2+ • • • Em (3) 

Similarly, all the top windings are con¬ 
nected in series aiding so that 

Ei — E x i + Eyi+Ezi + . .. E n i 

at balance 

Ei Ex i T . •. Em B -f- MS 

E« E X i~\-.. . E n » A-\-S—MS 

where 


(4) 

(5) 


Em — voltage output of the top winding of 
the iVth transmitter 

E,vi — voltage output of the bottom winding 
of the .ZVth transmitter 


If the cores are centered at 0.5 of maxi¬ 
mum stroke 

E xx =[E»-{k/2))+kx 
Em=[E,-(h/2)]+kn 
Ex 2 = [Fo+(A/2 )]—kx 

E n ^[E 0 +(k/2)]-kn (6) 


Substituting these values of voltages in 
the balance equation 5 yields 


N[Eo-(k/2)]+k(x+. . . n ) 
N[E 0 +(k/2)]-k(x+. .-n)~ 

B + MS 
A + S-MS 


(7) 


If the fixed resistors A and B are adjusted 
equal so that 

A~B~S[E*/k-('/*)] (8) 

The slide-wire rotation m is given by 


M= 


x -\-.. .n 


n — motion of the IVth transmitter 


(9) 


where A r =the total number of trans¬ 
mitters. 

For example, in Figure 7 (A) there are 
three variables X, Y, and Z. Conse¬ 
quently 


uA+I±^ ( 10 ) 

3 

If the ranges of the transmitters are not 
equal, the strokes of the cores are ad- 
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r —® 



Figure 5. Basic circuits of the resistance slide-wire 


justed proportionally. 

Subtraction is accomplished in a similar 
manner as shown in Figure 7(D) by con¬ 
nected opposite secondary windings in 
series. The difference between the totals 
of two groups of variables can be accom¬ 
plished by properly interconnecting the 
secondary windings. 


In the totalizing and subtracting cir¬ 
cuits shown by Figures 7(B) and 7(D) 
secondary windings are connected in 
phase opposition to produce voltages pro¬ 
portional to each transmitted variable. 
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Figure 6. Pickups of the slide-wire type for specific applications 
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AT BALANCE, eb = 0 
E 0 W » E x +E y + E z 

7n - x + y + z 


(B) totalizing -transformer balance 




AT BALANCE eb=0 
EoW = Ex “Ey 
m = x-y 



AT BALANCE eb= 0 
E| s E X i ^E Y , e B+77S 
E z E y | + Exe A+S~77JS 

7n -- X-Y 


C> RATIO CIRCUIT (01 SUBTRACTING (B SUBTRACTING 

SLIDEWIRE BALANCE TRANSFORMER BALANCE SLIDEWIRE BALANCE 


Figure 7. Basic computing circuits using adjustable core transformers 


The balancing follow-up is a similar trans¬ 
former unit. The proportion of the volt¬ 
age of each transmitter is adjusted by 
means of potential dividing resistors to 
the full-scale range. In Figure 7(B) it is 
assumed that all ranges are equal. 

In the ratio computing circuit shown by 
Figure 7(C) a fixed proportion of the 
linear voltage output of transmitter jc is 
compared against an adjustable propor¬ 
tion of the output of transmitter F. The 
proportion of the F output required to 
balance depends on the position M of the 
balancing slide-wire. A and B are cali¬ 
brating resistors determined by the range 
in ratios to be measured. 

Slide-wire Type 

The elementary computing circuits 
most frequently used with slide-wire 


transmitter and receiver elements are 
shown in Figure 8. These circuits are all 
of the a-c potentiometer* type and are 
generally preferable to Wheatstone bridge 
circuits because of: elimination of ef¬ 
fect of slide-wire brush contact resistance, 
and wide flexibility and ease of calibration 
using a minimum number of fixed ele¬ 
ments. Multiwinding transformers of 
shielded construction having negligible 
phase shift are used to supply proportional 
voltages to each computing slide-wire. 

Totalization and subtraction are ac¬ 
complished simply by algebraic addition 
of the voltage outputs of each transmit¬ 
ting slide-wire and comparing this total 
output against that developed by the 

*Note that these circuits are only pseudopotenti- 
ometric since they are basically ratio types and do 
not measure a-c voltages except in terms of a per¬ 
centage or portion of the supply voltage. 


balancing slide-wire in the receiver. 
Multiplication and division of transmitted 
functions are accomplished as shown in 
Figures 8(A) and 8(B) by connecting the 
computing slide-wires in cascade. To 
avoid the error produced by the loading 
effect of the second slide-wire, it is neces¬ 
sary that its resistance be high relative to 
the first. The deviation from linearity of 
the voltage output of the second slide is 
determined by the following analysis. 

The voltage across the second slide- 
wire is given by 

E a AM- 

C ~A + M-M 2 
E a = supply voltage 

ikf=per unit rotation of first slide-wire 
A = Si/Si where S% and Si are the resistances 
of the first and second slide-wire 
respectively. 


188 


Hornfeck—Computing Circuits and Devices for Industrial Process Functions 


July 1952 













_ S ^ _ 

> 


a 

| *> 

XI 

o 

\ 

\ 

\ 

I \”"■ ■■ 







S: fsl 



x 

o 

O UI x|v 

H II a 

XIV* N 
»N 

ui £ 

o 

z 



& 


o uj r 

HUH 

•gx >- 

* O 
UJ UJ 

o 

z 

< 

_J 

< 

CD 





> 

X 

o y-. 
o uj x 

M II » 

X> Nl N 
O o 
UJ 
UI 

o 

z 

< 



X 

o« 

O UJ X 
II ■ ■ 

>• 

•ui 




The theoretical linear voltage across 
the second slide-wire assuming its resist¬ 
ance to be infinitely large is 

et — EoM (12) 

The deviation from linearity is given by 


MH1-M) 

''A + M(l-M) 


For large values of A this error is a 
maximum when 

if=0.67 

so that 

Ae max==0. 15/A (15) 
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For example, if 
,4 = 100 

Ae max = 0.0015 = 0.15 per cent 

For a rather complete analysis of 
potentiometer loading errors refer to the 
paper by Nettleton and Dole. 3 

A voltage proportional to the square of 
the measured function is obtained simply 
by actuating two cascaded slide-wires 
from the transmitter as shown by Figure 
8(D). Consequently at balance the out¬ 
put and rotation of the motor-driven 
slide-wire is equal to the square of the 
transmitted variable. 

The square root of a function is ob¬ 
tained by using the circuit shown in 
Figure 8(E) in which the balancing motor 
operates simultaneously two cascaded 
slide-wires. Consequently, the voltage 
output of the receiver is proportional to 
the square of the motor shaft rotation. 
At balance this voltage must equal the 
linear voltage output of the transmitter so 
that the motion Y of the receiver is pro¬ 
portional to the square root of the meas¬ 
ured variable. 

Nonlinear Function Computer 

The system shown by Figure 9 has 
proved to be a powerful tool in metering 
single variables or combinations of proc- 


100 



ess function where: the transmitter mo¬ 
tion is a nonlinear function of the meas¬ 
ured variable, and the measured variable 
is to be converted into another function 
for computing purposes. The relation be¬ 
tween the balancing core motion Y and 
the servo motor shaft rotation M depends 
on the shape of the Team. Consequently 
M=F(Y). (F denotes the function speci¬ 
fied by the shape of the cam.) 

M=F(x ) (16) 

Since the motor shaft rotates in accord¬ 
ance with the desired function, a linear 
rise cam is used to actuate the linakge 
positioning the indicating, recording, or 
integrating elements of the receiver. Re¬ 
transmitting slide-wires actuated by the 
same motor reproduce the desired func¬ 
tion as linear voltages or resistances for 
control, computations, or further telem¬ 
etering. 

As demonstrated by the figure, this 
system is frequently used to convert a 
measurement of differential pressure, such 
as would be obtained by the transmitter 
of Figure 4(B), to fluid flow for linear re¬ 
cording. In other applications this sys¬ 
tem is used to convert a measured variable 
such as temperature to some nonlinear 
function of this variable such as the 
specific heat of a heat transfer fluid. 

If a slide-wire is used as the transmitter, 
the basic circuit of Figure 3(C) is inverted 
using a double cam, transformer balance 
follow-up as in Figure 9. 

The cam-operated transformer bal¬ 
ancer is generally superior to the tapered 
slide-wire follow-up for the conversion of 
functions because: cams can be cut to a 
greater degree of accuracy so that slide- 
wires can be wound to the desired func¬ 



RECEIVER 


PICK-UP 


Figure 9. Nonlinear function computer 


tions, and the range of taper of slide- 
wires is relatively limited whereas the 
cam rise rate can be made to cover a range 
from zero to almost infinity. 

A Typical Computer—Density 
Compensated Gas Flow Meter 

With the components and elementary 
circuits already described, systems of any 
desired complexity can be designed to 
perform a specific computation. An im¬ 
portant application of these to a system 
for recording pressure, temperature, and 
compensated gas flow is illustrated by 
Figure 10. The complete elementary cir¬ 
cuit for this metering system is shown in 
Figure 11. 

Equation 17 4 gives the weight rate of 
flow as a function of the differential pres¬ 
sure drop across the primary metering 
element. 

FF=358.9 CfD-iWhp (17) 

W —flow rate, pounds per hour 
D -2 = orifice diameter, inches 
h = differential pressure drop, inches of 
water at 68 degrees fahrenheit 
p = density of flowing fluid, pounds per cubic 
foot 

C and / are design constants 
For an ideal gas 
P Ts 

P = ps — — (18) 

PsT v 

ps = density at standard or design conditions 
Ps, Ts — design pressure and temperature 
P, T= operating pressure and temperature 

Substituting equation 13 in equation 12 
yields 

W=358.9CfDP^ h Ps^I? . (19) 

Combining constants in equation 14 


AT BALANCE, E 0 Y=E 0 X 

Y= X 

IF Y (RISE OF CAM , A')=777 rt 
(MOTOR SHAFT ROTATION) 

?n - y# =x^ 
Z-7H = X Ya.= CHART OR SCALE 
READING. 

EXAMPLE 

IF TRANSMITTER IS ACTUATED 
BY DIFFERENTIAL PRESSURE 
AS BY A BELLOWS TYPE 
METER, 

X s h = CONSTANT•( FLOW )* 
MAKE Y =77?*= X 

THEN 7U-- ~VT = VF* 

CONSTANT* FLOW. 
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Figure 10. Schematic 
diagram of density 
compensated gas flow 
rate meter 

Recorder records and 
transmits pressure and 
temperature, records 
compensated flow, and 
integrates compensated 
flow 

Computing receiver 
computes and transmits 
compensated flow 


TEMPERATURE SENSING 
ELEMENT 



W=K 3 — (20) 

The system is essentially an analogue 
computer which solves equation 20 by 
converting P, T, and h to electrical quan¬ 
tities and extracting the square root with 
the cam positioning nonlinear computer 
of Figure 9. 

The circuit of Figure 11 is balanced 
when the voltage developed across the 
pressure compensating slide-wire is equal 
to the voltage across the temperature 
slide-wire. 


If the cam which positions the core is 
shaped so that its rise is given by 

Y—M 2 

/- VhP 

M— V F= —— = W (22) 

Consequently the angular position M of 
the servo motor shaft is a linear function 
of compensated flow 

The specific arrangement of Figure 10 


illustrates how the basic elements can be 
arranged to suit the particular metering 
requirements. In this system correlating 
records of pressure, temperature, and 
compensated flow as well as an integra¬ 
tion of flow are required in a single meter. 
Space limitations require that the com¬ 
puting receiver employing the cam posi¬ 
tioned transformer follow-up be placed in 
a separate case. Compensated flow is re¬ 
transmitted from this computing unit to 
the recorder by means of a retransmitting 
slide-wire and corresponding receiving 
unit. If the integrator were not required. 


K YT - KXP — KhP (21) 
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Figure 12 (left). A 3-element receiver for recording pres¬ 
sure, temperature, and compensated flow rate with integration 

of flow 


Figure 13 (above). Nonrecording receiver for computing, 
controlling, and retransmitting 


the entire metering could be accomplished 
by a single, multifunction type of record¬ 
ing receiver. 

Figure 12 shows the 3-pen recording 
receiver presented diagrammatically in 
Figure 11. The motor-driven slide-wire 
receivers for temperature, pressure, and 
compensated flow drive color-coded re¬ 
cording pens to obtain continuous, identi¬ 
fiable records on the chart. The flow in¬ 
tegrator is of the mechanical type and 
permits determination of total flow over a 
specified period of time. The servoam- 
plifiers for the balancing circuits are 
mounted on the back of the instrument 
cases and are connected into their various 
circuits by means of plug and socket con¬ 
nectors for ease of servicing. 

Control 

In most applications the ultimate ob¬ 
jective of measurement is automatic con¬ 
trol of the process. There is increased ac¬ 
ceptance in the process industries of the 
basic electronically balanced null re¬ 
corder and its many variations as an 
analogue computer to obtain derived func¬ 
tions which are more significant to the 
process output than the primary vari¬ 
ables of pressure, temperature, level, 
flow, and so forth. The addition of auto¬ 


matic control to these telemeter com¬ 
puters is a further step in the complete 
mechanization of the process plant. 

A newly developed nonrecording type 
of computing receiver, incorporating in a 
relatively small instrument case retrans¬ 
mitting elements for either pneumatic or 
electric control, is shown by Figure 13. 
The use of the double cam balancing unit 
for the transformer follow-up and the ad¬ 
dition of two slide-wires for resistance 
functions makes the application of this 
unit relatively universal. The balancing 
amplifier has been made as small as pos¬ 
sible by the use of miniature tubes and 
functionally designed and selected com¬ 
ponents. It is mounted inside the case, 
but is easily removed for inspection and 
replacement. The motor-operated cam 
and slide-wire unit, the adjustable core 
transformer, and the pneumatic pilot are 
also subassemblies which can be removed 
readily and replaced. This type of re¬ 
ceiver will be used in applications where 
the available space is severely limited and 
where recording of the individual func¬ 
tions is not required. The addition’of 
control to the more conventional record¬ 
ing receivers, such as shown in Figure 12, 
is common. 

It should be obvious that time and 
space limitations prevent the description 


of the more elaborate and complete com¬ 
puting and control systems which would 
be used to meter and control a complex 
process. The number of combinations of 
metering and control which can be incor¬ 
porated even in one case, such as described 
in Figure 12, is quite large. Theoret¬ 
ically, there appears to be practically no 
limit to what can be done with these com¬ 
binations of relatively standard units and 
circuits. Unquestionably, there is an in¬ 
creased tendency for automatic supervi¬ 
sion by remote metering and control of 
various process functions, ranging all the 
way from the manufacture of cement to 
the generation of steam. The automatic 
computer permits further mechanization 
of these processes. For example, the 
generation of steam in the large central 
stations is accomplished almost com¬ 
pletely automatically to achieve high ef¬ 
ficiency. In the smaller power plants and 
in other processes, the desirability of com¬ 
plete automatic supervision depends on 
factors such as cost, reliability, safety, 
and the degree to which the special tech¬ 
niques and arts required in the process can 
be reduced to a science. 
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A Long-Distance Multipoint Telemetering 
System Using Teletype Transmission 


A. J. HORNFECK 

MEMBER AIEE 

Synopsis: This paper describes a point-to- 
point telemetering system using pulse width 
modulation and sequential scanning for 
transmission of several measured variables 
over a single circuit by teletype. The sys¬ 
tem, which is essentially electromechanical 
rather than electronic, operates on the basic 
principle whereby the magnitude of a 
measured variable is converted into a pulse 
of proportional time width. The system 
was developed specifically for remote moni¬ 
toring of pumping stations on a pipeline to 
permit automatic operation from a central 
dispatch office. In the design of the equip¬ 
ment emphasis was placed on requirements 
of reliability, sturdiness, simplicity, and 
ease of servicing rather than on factors such 
as high speed and excessive sensitivity. 

Electric Telemetering in Process 
Industries 

E LECTRIC telemetering has been 
used in the electrical industry for 
many years as a means of operating 
power systems, generating stations, and 
substations. 1 The application of electric 
telemetering to the field of industrial 
processes, while relatively more recent, is 
becoming increasingly important. Proc¬ 
esses which depend on telemetering for 
operation include the chemical industries, 
the making of iron and steel, the trans¬ 
mission and distribution of natural gas, 
the refining and transmission of oil and 
oil products, the operation of atomic re¬ 
actors, and many others. 

Short-range or local telemetering, fre¬ 
quently called remote metering, is com 
monly used inside large plants. Typical 
systems employ potentiometric and bridge 
balance circuits and electronically oper¬ 
ated instruments as receivers. Air- 
operated telemetering for these applica¬ 
tions offers strong competition to electric 


G. R. MARKOW 

ASSOCIATE MEMBER AIEE 

systems up to distances of several thous¬ 
and feet, particularly for the measure¬ 
ment and control of the nonelectrical 
variables such as pressure, temperature, 
level, and flow. The electric systems have 
the advantages of greater range and al¬ 
most instantaneous response, but since 
the accuracy and resolution is affected by 
the attenuation of the transmitted signal, 
the longest distances normally used are 
approximately 5 miles for a-c systems and 
up to 50 miles for d-c potentiometric sys¬ 
tems. In applications where the system 
involves distances of the order of 100 or 
even 1,000 miles as, for example, an oil 
pipeline system, measured quantities are 
converted into some factor of time such 
as the frequency of a periodic signal or 
the time duration of a pulse. 

The system described in this paper, 
while utilizing some of the devices of 
local telemetering, is of the long-distance, 
time-factor type. Its operation depends 
on a number of principles which individ¬ 
ually are not new in the present art of 
telemetering 2 but in combination result in 
a system of considerable flexibility which 
is capable of telemetering measured quant¬ 
ities over commercial telephone or tele¬ 
graph circuits or by radio link any desired 
distance. 

The design of this system is basically of 
an electromechanical nature utilizing a 
minimum number of electronic compo¬ 
nents. Electronic amplifiers and relays 
actually are used only in the local re¬ 
ceivers and in the ‘ 'sweep balance” de¬ 
tectors, since the requirements of the 
specific application as to speed and res¬ 
olution did not warrant an all-electronic 
system. However, it should be apparent 
that the general principle of operation of 


the system is not restricted to the specific 
design which was developed for the pipe¬ 
line application. 

This long-distance telemetering system 
was developed as part of a co-operative 
project to provide completely automatic 
operation for pumping stations on a 
products pipeline of the Shell Oil Com¬ 
pany running from Wood River, Ill., to 
Toledo and Columbus, Ohio. These 
pumping stations are controlled from the 
central dispatching office located in New 
York City. 

The terminology employed in this paper 
is taken both from the point-to-point and 
the mobile or radio telemetering fields 
and the choice of a term used in describ¬ 
ing a principal function is selected on the 
basis of its accuracy and conciseness. It 
is important to note that basic principles 
of operation are common to both groups 
of telemetering and the differences in 
equipment exist principally because of 
the special techniques that are utilized 
in the two general systems. These dif¬ 
ferences in techniques depend primarily 
on the differences in requirements which 
exist between these two general fields of 
application. From the technical stand 
point it is believed that mobile telemeter¬ 
ing to a considerable extent has become 
the leader principally because of its 
utilization of the high speed of response, 
compactness of design, and mobility of 
the radio link. Consequently, it is ex¬ 
pected that mobile telemetering will in¬ 
fluence the design of point-to-point or 
ground systems to a considerable extent 
in the future. Furthermore, it is prob¬ 
able that in systems developed for telem- 
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Figure 1. Block diagram 


etering at very high speeds a multitude 
of quantities, important to the flight test 
of a controlled missile, will be applied in 
modified forms to industrial telemetering 
systems. Obviously, differences of design 
will continue to exist since factors such as 
size, weight, and high-frequency re¬ 
sponse are not critical in the stationary 
systems, whereas long life, reliability, 
ease of service by nonelectronic person¬ 
nel, and other factors are important in the 
industrial field. 

The Basic System 

The telemetering system described in 
this paper is essentially a combination of 
local telemetering units of the servo- 
actuated position type 3 coupled to a long¬ 
distance system which operates on a 
general principle of what the mobile telem¬ 
etering engineers call a “time division, 
pulse coded multiplexing system.” Ac¬ 
tually, it is also what the point-to-point 
telemetering engineers designate as an 
impulse duration multiplex or multipoint 
system. In this scheme each measurand 
(measured quantity) is sampled in a se¬ 
quence by means of a scanner and then 
transmitted to the final receiving station 
ar stations as a series of coded characters. 
Each channel of information is identified 
by the code selected and each measure¬ 
ment of the various quantities is sampled 
in sequence and transmitted as a com¬ 
plete series of coded characters by con¬ 
verting the width of the pulse as deter¬ 
mined by the scanner into a proportional 
number of coded characters. Figure 1 is 
the basic diagram of the system. The 
measurands may be electrical but in proc¬ 
ess industries are more frequently non¬ 
electrical quantities such as pressure, flow, 
level, and so forth. Each of the quanti¬ 
ties is converted either directly or through 
several stages of local telemetering into 
an electrical quantity proportional to the 


of the telemetering system 


original variable. The outputs of all the 
transducing or converting elements in the 
system are of the same nature so that they 
can be compared against a common 
scanning function. The sequence switch, 
commonly called the commutator or 
sampler, continuously and in a predeter¬ 
mined sequence connects each of the 
transducers into the scanning circuit 
where a comparison is made between 
their outputs and the time varying output 
voltage of the scanner. If the output of 
the pickup is a voltage directly propor¬ 
tional to the measurand, the output of the 
scanner during the rise period is made a 
linear function of time. As will be ex¬ 
plained in detail later, the function of the 
scanner is to convert the magnitude of the 
output of each transducer into a propor¬ 
tional time pulse. 

The proportional time interval is rep¬ 
resented by the time of dwell of an elec¬ 
tric contact in this system but also may 
be represented by the period of conduc¬ 
tion of an electron tube. This pulse may 
be transmitted without further modifica¬ 
tion over a wire link. If a radio link is 
used, an amplitude- or frequency-modu¬ 
lated carrier is required. 

The transmitted quantity may be iden¬ 
tified at the receiver by its chronological 
position in a train of impulses such as 
might be indicated on a cathode-ray 
oscillograph. In the specific system, de¬ 
veloped for the pipeline telemetering, 
identification is accomplished by convert¬ 
ing each pulse into a group of coded tele¬ 
type characters. The coding is accom¬ 
plished by synchronizing the teletype 
equipment with the sequence switch. 

Description of the Function of the 
Elements in the System 

Pickups 

There are many primary detectors or 
pickups in use today whose purpose is to 
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convert the measurand to some electr 
quantity which can be measured convi 
ently. Figure 2 shows a number of 
transducers, which convert the pro< 
variables such as pressure and flow 
proportional alternating voltages 
means of differential transformers, wt 
output voltages are proportional to tl 
core position. 3 

The pickups shown by this figure 
employed where the basic measuring 
ments are mechanical devices such s 
bellows, Bourdon tubes, floats, and 
forth. This type of mechanical link \ 
duces a position or motion proportic 
to the measurand which then is convei 
into an electrical variable by means of 
adjustable core transformer coupled 
the mechanical device. It will be nc 
that these systems represent two 
sometimes more stages for the conver: 
of the measurand into an electrical v 
able. On the other hand pickups sucl 
thermocouples and resistance therm' 
eters are primary detectors which ] 
form the conversion in a single stage. : 

Figures 2(A), (B), and (C) show p 
ups whereby pressure, differential p 
sure, and flow respectively may be m 
ured and converted to proportional \ 
ages. For variables such as pressure j 
level the arrangement shown in Fiji 
2(D) is suitable if a record is desire< 
the transmitter. In the case of nonlii 
variables a secondary pickup is emplc 
as indicated in Figure 2(E). Here; 
primary detector is balanced by a mo 
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•driven slide-wire unit and retransmitted 
to the differential transformer by means 
•of the correct nonlinear function cam to 
give an output voltage which varies 
linearly with the measurand. 

Figure 3 shows several pickups of the 
resistance type. The resistance thermom¬ 
eter and the strain gauge are primary 
detectors in which the resistance varies 
directly with the measurand. The pres¬ 
sure pickup shown in Figure 3(C) func¬ 
tions by positioning the slider of a low- 
torque potentiometer. Figure 3(D) shows 
the manner in which a retransmitting 
■slide-wire is actuated by the motor of a 
self-balancing recorder. 

The Scanner 

The scanner shown in Figure 4, con¬ 
sists of a synchronous motor driving a 
linear rise cam which actuates the core 
■ of the differential transformer in such a 
way that its output voltage varies pro- 
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Figure 4. Elementary diagram of the scanner 



portionally with time. A second cam, 
mounted on the same assembly as the 
scanning cam, actuates the timing micro - 
switch. The two cams are so arranged 
that at the instant when the differential 
coil output voltage is passing through 
zero the micro-switch will close and will 
remain closed until after the scanning 
voltage has exceeded 100 per cent. The 
opening of the timing micro-switch 
initiates the switching sequence. After 
an interval required for this operation it 
closes once more to begin the timing cycle 
for the next variable. 

The secondary of the scanning trans¬ 
former is connected in turn to each of the 
pickup transformers by a sequence switch, 
or sequentially operated relays. The 
scanning transformer and the pickup 
transformer are connected to obtain a 
voltage proportional to the difference of 
their outputs. This difference or error 
voltage is applied to the electronic relay. 
This relay is so arranged that when the 
input voltage to it is of one phase its out¬ 
put relay coil is de-energized, and at the 
instant the phase changes the relay is 
energized, closing its contacts. 

The scanner operates in the following 
manner. Starting at time equal to zero, 
when the timing micro-switch closes, the 
scanning voltage is zero and begins to in¬ 
crease linearly with time. The voltage to 
the electronic relay is a maximum at 
T=0 since 

E r — Ei n — E s and at jT=0, E s = 0 

This relay voltage will decrease linearly 
with time as the scanning coil voltage in¬ 
creases and approaches the value of the 
output voltage of the pickup. At the in¬ 
fant when the two voltages are equal, 
e electromagnetic relay on the output 
the electronic unit will be de-energized 
us opening the circuit to the transmitter 
r transmission line and completing the 
period of transmission of intelligence. 

The curves of Figure 5 indicate dia- 
grammatically the operation of the elec¬ 
tronic relay. Three curves of relay input 
voltage (for three different signal inputs) 
on a time base are shown in Figure 5(B) 


showing how the magnitude of the signal 
voltage determines the pulse length. 
Figure 5(C) depicts one cycle of the 
scanning voltage. 

Over the measurement interval the ref¬ 
erence voltage of the scanner is given by 
the equation 

e s /E s =t/T 0 (1) 

where 

e s = the voltage output of scanner at time t 
£ s = the voltage output of the scanner at 
time To, the end of the scanning 
period 

The output of the linear transducer is 
given by 

?x/Ex = x/X 0 (2) 

where 

e x ~ voltage output of transducer at a value 
of variable x 

E x — voltage output of transducer at a value 
of variable X Q (100 per cent) 

The duration of the measurement con¬ 
tact is obtained by equating e s to e x and 
solving for t. Hence 

t = (E x /E s )(x/X a )T 0 (3) 


If E S = E X 

t = (x/X 0 )T 0 (4) 

Equation 4 shows that the duration of 
the measurement contact is directly pro¬ 
portional to the measurand x. 

If the output of the transducer is a non¬ 
linear function of the measurand, this 
function can be extracted by using a 
servo-actuated receiver to operate 
through a cam, the core of a retransmit¬ 
ting transformer. For example, the volt¬ 
age output of the differential transmitter 
shown by Figure 2(B) is proportional to 
the square of flow. In this case the re¬ 
transmitter cam has a square root rise 
characteristic so that the output of the 
retransmitter transformer will be pro¬ 
portional to flow. If, however, all of the 
quantities to be scanned have the same 
nonlinear characteristic, the cam of the 
scanning unit may be shaped so that its 
rise is a similar nonlinear function of 
motor shaft rotation. For example, in the 
case of the differential-type flow meter 

Cx — E x (.x/X q) 2 (5) 

If the rise of the scanning cam is made 
proportional to the square of the angle 
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Figure 7. Simplified 
block diagram of station 
arrangement 



It is entirely feasible to use a rotating 
slide-wire as the scanner device, and for 
high-speed telemetering a calibrated elec¬ 
tronic generator could be used. For the 
requirements of the pipeline telemetering 
system the cam-operated scanner has a 
number of basic advantages such as elim¬ 
ination of sliding contacts, ruggedness, 
and ease of calibration. 

Transmission Link 

One method by which the intelligence 
iay be transmitted and coded is by the 
i of the teletype system. While a de¬ 
led explanation of the various circuits 
id to achieve this is beyond the scope 
x the paper, the general operation of the 
system is given in the following para¬ 
graphs. 

An electric contact whose duration is 
proportional to the measurand is supplied 
by the electronic balance detecting relay 
of the scanner to the teletype transmit¬ 
ting equipment. The teletype apparatus 
sets up coding letters for each variable, 
such as is shown next, using the letter S 
for suction pressure, T for temperature, 
F for flow, L for level, and so forth. The 
equipment starts printing at a constant 
rate the appropriate letter at the begin- 
ling of the scanner timing cycle and stops 
vhen the scanner contacts open. Thus, 
he number of letters in any one sequence 
proportional to the variable. At the 


end of the transmission of intelligence, a 
code letter such as PQ is sent to indicate 
the location of the transmitter in cases 
where there may be similar intelligence 
channels from other remote stations. 

SSSSSSSSSSSSSSSSSSS PQ 
TTTTTTTTTTTTTTTTTTTTTTTTT 

PQ 

FFFFFFFFFFFFFFF PQ 
LLLLLLLLLLLLLLLLLLLLLLLLLL- 
LLLLLLLLLLLLL PQ 

The teletype channel has the advantage 
of using comparatively standard equip¬ 
ment and easily furnishing the intelligence 
required both at the transmitter and re¬ 
ceiver and all other teletype sets con¬ 
nected to the particular network. It has 
the disadvantage that the resolution is 
limited by the speed at which the tele¬ 
type may be operated. In most cases, 
however, this is not a serious drawback. 

Pipeline Application 

This system was first used on automatic 
p um ping stations for a pipeline of the 
Shell Oil Company. It is installed, as 
’ shown in Figure 6, on the four automatic 
booster stations located at New Douglas, 
Dennison, and Effingham, Ill., and North 
Salem, Ind. The stations make a much 
higher throughput possible than usual in 
an 8-inch line as well as increasing over¬ 
all operating efficiency. 

The operation of the four stations is 


entirely controlled from the central dis¬ 
patch center in New York City. Various- 
products are pumped in turn through the 
line and it has been found that the slight 
intermixing produces no harmful results 
in the quality of any of the products. The 
location of these products is traced on a 
control board which simulates the pipe¬ 
line. Thus, the dispatcher in New York 
has a completely integrated picture of the 
location and quantity of all products in 
the whole pipeline as well as the pressures 
and pump motor current of each of the 
four booster stations. Since the dis¬ 
patcher has control of these stations, as 
well as complete operating information, 
he is in a position to operate the line in a 
more efficient manner. 

The telemetering equipment, as pre¬ 
viously described, furnished information 
regarding suction pressure, discharge 
pressure, and pump motor current. A 
block diagram, Figure 7, indicates the 
telemetering equipment supplied for each 
of the four booster stations. The suction 
and discharge pressures are measured by 
gauges installed in the pump house and 
locally telemetered to the control room 
where they are recorded. A motor cur¬ 
rent recorder also is located on the same 
panel. Electrically isolated voltages from 
the three transmitting differential trans¬ 
formers are connected to the sequence re; 
lay and in turn to the scanner for long 
distance telemetering. A simplified sche- 
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matic diagram, Figure 8, outlines the 
circuits used in the measuring and telem¬ 
etering equipment supplied by Bailey 
Meter Company. 

To obtain readings of telemetered 
variables, the dispatcher at New \ork 
City dials a code number which operates 
the equipment at any one of the four 
chosen stations and automatically telem¬ 
eters the pressures and pump motor 
current in the form as shown previously. 
At the end of one telemetering cycle of the 
three variables, the equipment automati¬ 
cally shuts off and the dispatcher, if he de¬ 
sires, then may obtain similar informa¬ 
tion from the three other stations. The 
dispatcher also may start and stop the 
main pump motor, lubricating oil pump, 
suction and discharge valves, and other 
necessary equipment. In an emergency, 
he also may sound a siren to summon the 
operator who is located at each station 
for maintenance and emergency purposes. 

All the information received in the New 
York dispatcher’s office is simultaneously 
printed on all the other teletype machines 
which are connected to the same system. 
The teletype machine is operated at a 
speed of 368 characters per minute or 163 
milliseconds per character. The total 
cycle period for the scanning of one meas- 
urand is 15 seconds, or 45 seconds for a 
complete intelligence report from any one 
station. The 100-per-cent measurement 
time was selected at 60 characters or 
9.783 seconds. Thus, the resolution of 


the system is one part in 60 or approxi¬ 
mately 1.7 per cent. The remaining por¬ 
tion of the cycle, 5.2 seconds, is allowed 
for the various switching operations in¬ 
volved to transfer from one pickup to the 
next, change code letters, and so forth. If 
a lower resolution were permissible, the 
cycle time could be reduced proportion¬ 
ally, and similarly for a higher resolution, 
a longer cycle would be employed. 

Operating Experience 

This system has been in operation since 
September 1950. Some trouble was ex¬ 
perienced at the beginning with servo- 
amplifier failure, caused by faulty com¬ 
ponents, but this was easily remedied 
since these units are standard plug-in as¬ 
semblies. Since then the equipment has 
been found to be completely reliable due 
to its inherently simple nature, lack of 
delicate components, and over-all rugged¬ 
ness.- It serves its purpose admirably in 
integrating all information and super¬ 
visory control at one centralized point and 
thus permitting much more efficient 
utilization and operation of the pipeline 
system. 

Conclusions 

It will perhaps appear obvious to the 
electronic engineer in the field of mobile 
telemetering that the utilization of elec¬ 
tronic techniques could have been made 


IIO vr 

60 v 

in any or all of the elements of the system 
described. For example, an electronic 
scanning system could have been used in 
place of the electromechanical system, re¬ 
sulting in much greater speed of scanning. 
Acceptance of complex electronic systems 
for the control and telemetering of indus¬ 
trial process variables has been somewhat 
slow, particularly where the requirements 
do not exclude the use of simpler electro¬ 
mechanical mechanisms. Where high 
speed and sensitivity are required and a 
great number of variables are to be telem¬ 
etered, the electronic systems probably 
would be required. 
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Subbarmonics in a Series Nonlinear 
Circuit as Influenced by Initial 
Capacitor Charge 


WILLIAM J. McKUNE 
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S ERIES circuits containing a capaci¬ 
tance and nonlinear inductance have 
assumed considerable importance- in 
modern engineering design. However, 
knowledge of the behavior of such cir¬ 
cuits is extremely limited because of the 
inherent nonlinearity of the differential 
equations encountered. 

While considerable work has been done 
by the investigators of nonlinear vibra¬ 
tions and celestial mechanics on nonlinear 
differential equations in which the non¬ 
linearity is involved in the damping term 
(the equation of van der Pol) or in the 
displacement term (Buffing's equation) 
or in the forcing function, very little has 
been done on differential equations in 
which the nonlinearity is involved in the 
inertia term. The subject is further com¬ 
plicated by the fact that transients in 
such networks may be extremely long so 
that solutions by a differential analyzer 
are lengthy and subject to cumulative 
errors. 

Plence it is felt that at present the most 
reliable method of predicting with a 
reasonable degree of accuracy the per¬ 
formance of proposed circuits is from a 
wealth of experimental data. The data 
now available give considerable insight 
into the type of response that may be en¬ 
countered. From the work of Rouelle, 1 
McCrumm, 2 and others, 3-5 it is known 
that three different types of current re¬ 
sponse may be encountered in a series cir¬ 
cuit containing a resistance, a capaci¬ 
tance, and a nonlinear inductance upon 
shock excitation. The first of these is the 
normal fundamental-frequency exciting 
current of the inductance. The second is 
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a much higher value of current of funda¬ 
mental and higher frequencies, which has 
been labeled ferroresonance. The third 
type of response contains currents which 
are integer submultiples of the applied 
frequency. This phenomenon is called 
subharmonic response, and the order of 
this response indicates the ratio of the ap¬ 
plied frequency to the lowest frequency 
component of the resulting current. 

The work of these investigators gives 


information as to the ranges of circuit 
parameters for which the different types 
of response may be obtained. Angello 6 
reported the effect of the initial switching 
angle on the occurrence of subharmonics 
with constant initial capacitor charge. 
However, the effect of initial capacitor 
charge or of initial flux linkages has not 
been investigated by previous authors. 
The purpose of the investigation reported 
in this paper was to determine the effect 
of the initial capacitor charge on the oc¬ 
currence of subharmonics, with the initial 
switching angle and initial flux linkages 
held constant. The investigation was re¬ 
stricted to third-order subharmonics as 
this is the order most likely to be en¬ 
countered. 


Test Circuit 

The circuit investigated is shown in 
Figure 1. Figure 1(A) shows the non- 
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Figure 1. Test circuit 
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Effective Open Circuit Current 


Figure 2 (left). 
Saturation curves 
of inductive ele¬ 
ments at funda¬ 
mental frequency 


Figure 3 (right). 
Hysteresis loops 
at 57 volts rms 




1-T 


linear circuit containing a resistance R, a 
capacitance C, and a nonlinear inductance 
(a transformer with open secondary) 
connected through a relay contact X to a 
low-impedance variable voltage source. 
The capacitor was charged prior to clos¬ 
ing the relay contact by applying a uni¬ 
directional voltage. Inasmuch as the oc¬ 
currence of subharmonics was found to be 
influenced by initial flux linkages, the 
transformer was demagnetized by ap¬ 
plying a decreasing-amplitude sinusoidal 
voltage to the secondary of the trans¬ 
former prior to each reading. 

A high-persistence oscilloscope was 
used to set the relay contact X to close at 
the desired angle and to monitor this 
angle continuously. To observe the 
angle it was necessary to blank the oscil¬ 
loscope approximately one cycle after re¬ 
lay X operated. It was found that if the 
beam of the oscilloscope was left on prior 
to operation of the X relay the initial 
trace on the screen was objectionable. 
This was eliminated by applying a nega¬ 
tive voltage to the grid of the cathode-ray 
tube, removing this bias just prior to the 
closing of X, and putting it back approxi¬ 
mately 1 cycle later. Because of the 
high potential to ground of the grid cir¬ 
cuit, it was necessary to isolate this cir¬ 
cuit from the X relay. 

The circuit used for control of the 


switching angle is shown in Figures 1(B) 
and 1(C). The grid of the 6SJ7 tube was 
supplied through a phase-shifting net¬ 
work such that the phase position of the 
approximately square wave output of the 
6SJ7 tube could be varied with respect to 
the applied voltage. This output was ap¬ 
plied to the grid of the thyratron through 
a differentiating network with a time con¬ 
stant of about 5 microseconds. The 
phase position of the resulting pulse 
could thus be regulated by means of the 
phase shifter in the grid of the pentode. 
The magnitude of the pulse was controlled 
by adjusting the input to the phase 
shifter and was set to such a value that 
the positive pulse would raise the grid 
voltage of the thyratron slightly above 
cutoff. Thus the thyratron was trig¬ 
gered once every cycle with a pulse of 
very short duration and of variable phase 
position and would fire at this point if the 
plate circuit was closed. 

Method of Taking Data 

The data were taken as follows. The 
residual flux in the inductance was de¬ 
stroyed and the desired potential was 
placed on the capacitor. As soon as the 
initiating switch was closed, relay Y 
closed, which unblanked the oscilloscope 
and completed the plate circuit to the 



1 —,—,—,—|—| j | i r 

thyratron. The next positive triggering 
impulse from the pentode fired the thyra¬ 
tron, energizing relay X. Relay X then 
energized the nonlinear circuit and picked 
up relay Z which again blanked the oscil¬ 
loscope. The proper time interval was 
obtained by adjusting the spacing of the 
relay contacts. 

Since the work of Angello® indicated 
that the probability of subharmonic re¬ 
sponse was greatest at zero switching 
angle, it was decided to energize the cir¬ 
cuit at this point. The wave form of the 
current in the nonlinear circuit was ob¬ 
served on a low-persistence oscilloscope 
connected across the resistive elejnent. 

Two different inductances were used in 
order to determine the effect of different 
saturation curves. The saturation curves 
for the two transformers are shown in 
Figure 2. The hysteresis loops for the 
two transformers at an applied voltage of 
57 volts rms, GO cycles, are shown in 
Figure 3. 

Three tests were made with transformer 
number 2: one with constant capacitance 
and variable applied voltage, and two 
with constant applied voltage and vari¬ 
able capacitance. With transformer 
number 1, one test was made with con¬ 
stant capacitance and variable voltage 
and one with variable capacitance and 
constant applied voltage. .The values of 
voltage and capacitance which were held 
constant were selected because in a rough 
exploration of the complete regions of re¬ 
sponse they proved to be the values most 
conducive to subharmonic response. 

The results obtained are shown in 
Figures 4 through 8 and Appendix I 
shows the conditions under which they 
were obtained. 

Conclusions 

Several facts are worthy of note. In 
his experimental work Rouelle 1 found 
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Figure 4. Areas of stable response. Transformer number 1—applied voltage, 60 volts rm 


that, for random excitation, subharmonic 
response could not be obtained as easily 
for a transformer with a sharply curved 
saturation curve as for one in which the 
saturation curve had less severe curva¬ 
ture. 

The investigation herein described 
supports this finding as can be seen from 
a, comparison of the areas for the two 
transformers. 

In his work McCrumm 2 found that 
for certain regions of capacitance, sub- 
harmonics could not be obtained for cer¬ 
tain definite values of applied voltage al¬ 
though they could be obtained for values 
of applied voltage above or below this 


value. Dehors’ work 7 supported this 
conclusion. However, although a lengthy 
search was made on the two inductances 
used in the present investigation, this 
gap could be found on neither. One pos¬ 
sible explanation is that transformer 
number 2 had a high value of conductor 
resistance; McCrumm pointed out that 
high values of resistance could destroy 
this "lake.” And the curvature of the 
saturation curve for transformer number 
1 may have been so sharp that this 
phenomenon may not have existed. 

Angello 6 found that the initial flux 
linkages in the inductances had no effect 
on the wave forms of stable subharmonics. 


This was confirmed; however, the initial 
flux had a definite effect on whether or 
not subharmonics would occur. 

It should be pointed out that the re¬ 
sponse in the immediate neighborhood of 
ferroresonance was difficult to determine 
because of the magnitude of the transients 
encountered. In no case did a transient 
ferroresonant current yield a steady-state 
subharmonic current. 

Several authors have reported that a 
prolonged subharmonic transient may be 
obtained. This was observed here, with 
yet another phenomenon. In certain re¬ 
gions the current for the first few cycles 
was that of normal exciting current which 
would then grow into the higher magni¬ 
tudes of subharmonic response. 

Subharmonics of the second, fourth, 
sixth, and ninth orders were also ob¬ 
tained, although no attempt was made to 
map the regions except to note that they 
were extremely small. It is of interest to 
note that the sixth- and ninth-order sub¬ 
harmonics were always accompanied by a 
predominant third-order subharmonic. 
Similarly, the fourth order was always ac¬ 
companied by a predominant second 
order. 

It was hoped that the results would be 
of such a type that definite conclusions 
could be reached regarding a means of 
avoiding ferroresonance and subhar¬ 
monics by a proper initial charge on the 
capacitor. In the case of subharmonics 
no such conclusions could be drawn. It 
can be seen, however, that by putting- 
zero charge on the capacitor, ferrores¬ 
onance is not likely to be obtained 
whereas it may be obtained with a high 
value of initial charge. The placing of a 
resistor in parallel with the capacitor as 
suggested by Butler and Concordia 8 ac¬ 
complishes this by allowing a path over 
which the residual charge may be dis¬ 
sipated. 


Appendix I 

Frequency of applied voltage = 60 cycles per 
second 

.Series resistance A = 0.80 ohm 

Switching angle=0 degree 

Initial flux linkages =0.0 j 

Transformer number 1 =2-kva distribution 
transformer, Hypersil core, 0.22-ohm 
conductor resistance 

Transformer number 2=power supply plate 
transformer, 2.20-ohm conductor re¬ 
sistance 

Sign convention = switch closed with the 
voltage rise from a to b, see Figure 
1(A), passing through zero from nega¬ 
tive to positive 

=charge on capacitor C 
arbitrarily defined as positive when d 
was positive with respect to e 
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In Figures 4 through 8, the area of stable 
third-order subharmonic response is repre¬ 
sented by dots. The area of stable ferro- 
resonant response is shown by crisscross 
lines and the area of normal exciting current 
is clear. 
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Solution of Electrical Engineering 
Problems by Southwell’s 
Relaxation Method 


MISS E. M. GRAD 

NONMEMBER AIEE 


Synopsis: The solution of simultaneous 

linear algebraic equations by the basic re¬ 
laxation method is described in detail and 
the use of short-cut methods is explained 
and demonstrated. It then is shown that 
the same technique can be applied to the 
solution of ordinary and partial differential 
equations, and its accuracy is discussed. 
Several problems taken from various 
branches of electrical engineering are solved 
to illustrate the application of the method. 


T HE method of “relaxation” as de¬ 
vised by Southwell 1 ”* for the solution 
of framework problems has proved to be 
of very wide application. Used first in the 
numerical solution of systems of simul¬ 
taneous linear algebraic equations, it now 
lias been extended, with outstanding suc¬ 
cess, to the solution of ordinary and par¬ 
tial differential equations. The tech¬ 
nique is essentially the same for all equa¬ 
tions, and many otherwise quite intract¬ 
able problems have been solved by this 
method. Occasionally conditions are met 
in which the successive solutions con¬ 
verge only very slowly, or even not at all. 
This question of “ill-conditioning” in re¬ 
lation to the relaxation method was dis¬ 
cussed by Fox. 4 

In electrical engineering the most 
likely applications of the method are to 
the solution of sets of simultaneous equa¬ 
tions arising, for example, in the analysis 
of coupled circuits, and to the solution of 


electrical and magnetic field problems, de¬ 
termined by the partial differential equa¬ 
tions of Poisson and Laplace, or the wave 
equation. By other methods, these dif¬ 
ferential equations can be solved only in 
simple or in “idealized” cases, but relaxa¬ 
tion can give a numerical solution for most 
boundary shapes and conditions. 

Several accounts of the method and its 
application have been published 1 " 7 but 
few deal with the actual technique. The 
purpose of this paper is to present an out¬ 
line of the basic principles of the method 
and to demonstrate the technique by ex¬ 
amples taken from various branches of 
electrical engineering. 

Solution of Simulataneous Linear 
Algebraic Equations 

Basic Process of Relaxation 

Consider as an example of the method 
the pair of equations 

79-3x+y=0 W 

147 -f-x — 3y=0 

Table I. Operations Table 


These equations may be written in the 
form 

iq=E79-3.%-+y ( 3 ) 

i?o=147+x-3y (4) 

where, for the correct values of x and y, 
R x = R, = 0. For any other values of x and 
y, Ri and R* are not zero. Ri and R» are 
called the “residuals” and are a measure 
of the error of the solution for any given 
x and y. Ri and R» can he calculated for 
any pair of x and y, but one wants to find 
that unique pair which makes R\ = R-i = 0. 

First, an operations table is set up. 
From equations 3 and 4 it is seen that unit 
positive increment to x changes Ri by —3 
and Ri by -ff 1 and unit positive increment 
to y changes Ri by +1 and R« by —3. 
This is given in Table I. 

One has to make an initial guess of x 
and y to start the computation (a good 
initial guess is helpful as it will shorten 
the subsequent numerical work, but is by 
no means essential). In this particular 
example, take for convenience x = 0 and 
y = 0 as initial values. From this one cal¬ 
culates the residuals Ri and lousing equa¬ 
tions 3 and 4 and enters these values tn 
the relaxation table, Table II. 

The biggest residual then is picked out, 
which is 147 in this case, and it reduces to 



ARi 

ARi 


.-3. 

.... 1 

Ay = 1. 

_ 1. 

. . . . -3 
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Tabic II. Relaxation Table 




Ri 

Ri 

s = 0. 

... y = 0. 

. 79_ 

. ... 147 


Ay = 49. 

.128_ 

. . . . 0 

A a; = 42. . . . 


. 2.... 

. . . . 42 


1) 

h- 1 

. 16_ 

. . . . 0 

A* = 5. 


. 1 

5 


Ay = 1... 

. 2_ 

2 

A* = 1. 


. -1_ 

. . . . 3 


Ay = 1. 

. 0_ 

. . . . 0 

£ = 48. ... 

... y = 65 



Check Ri — 

79-144 + 65 = 0. 




147+48-195 = 0. 




zero. Since Ay = 1 gives AR- 2 — — 3, Ay— 
49 gives AR 2 = —147, thereby making R« 
zero, but at the same time it makes ARx — 
49, so increasing R x to 128. This opera¬ 
tion is entered in the second line of the 
Relaxation Table. Proceeding in the 
same way, take Ax =42, giving ARi = — 
126. This does not make R\ quite zero, 
but as near to zero as necessary. The 
corresponding increment to Ri — 42, and 
this operation is written as line 3. The 
working is continued in this manner until 
both residuals are zero. (In line 7 a nega¬ 
tive residual occurs, but this presents no 
difficulty.) The increments to x and y 
then are summed and the values so ob¬ 
tained substituted in the original equa¬ 
tions, to check that Ri=R 2 =Q. 

If a mistake had been made in Table 
II, Ri and R 2 would not be found to be 
zero, but some other small numbers. 
These values then would be taken as new 
starting values for the residuals and “re¬ 
laxed” to zero, the further increments re¬ 
quired to do this being added to the previ¬ 
ous ones to give the correct solution. 

Relaxation is used to get a numerical 
solution as quickly as possible. The 
method involves a large number of very 
simple operations. Since the results can 
be checked easily during the progress of 
the work, it is better to work quickly and 
risk arithmetical errors, making checks of 
the residuals only occasionally, than to 
stop and check each step by which the re¬ 
siduals have been calculated. To make 
the method even more speedy, several 
shortening devices have been developed. 

Shortening Devices 

Over relaxation. From the example it 
can be seen that whenever one residual is 
decreased the other is increased though 
by a smaller amount. If, therefore, in¬ 
stead of making just sufficient increment 
to make a residual zero, one makes a 
larger increment, that is, one overrelaxes 
to make the residual negative, then when 
it receives an increment as a result of an¬ 
other operation it still will be near zero. 


This can be seen from the alternative 
method, Table III, of working the previ¬ 
ous example. 


Table III. Alternative Method 




Ri 

Ri 

x = 0 . 

... y = 0. 

... 79.... 

.... 147 


Ay = 60.... 

...139.... 

. . .-33 

As = 45_ 


... 4_ 

.... 12 


Ay = 5. 

... 9. .. . 

. . . . -3 

As = 3. 

s = 48_ 

... y = 65 

... 0.... 

. . . . 0 


Table III, of course, gives the same 
answer but it took fewer lines of work¬ 
ing. 

Block Relaxation. By definition this is 
the application of unit increment to more 
than one of the variables at a time. That 
is a superposition of the effects Ax= 1 and 
Ay= 1, giving a new line in the operations 
table shown in Table IV. 

Table IV 


A2?i ARi 


As = 1.—3. 1 

Ay= 1. 1.—3 

As = Ay= 1.—2.—2 


When the total residual is zero, individ¬ 
ual residuals must he evenly distributed 
among positive and negative residuals, 
and such residuals are much easier to 
“liquidate” than a set which all have the 
same sign. Block relaxation is very use¬ 
ful in reducing the total residual (that is, 
the sum of all residuals) to zero. To dem¬ 
onstrate this the same example is again 
solved. 


Tabic V 




Ri 

Ri 

* = 0. 

... y = 0. 

_ 79. ... 

.... 147 

As = 56.... 

. . . Ay = 56.... 

...-33.... 

. . .. 35 


Ay = 9. 

.. .-24_ 

.... 8 

As= -8... 


_ 0.... 

.... 0 


Group Relaxation. By definition this is 
the application of unequal increments to 
more than one of the variables at a time. 

Any new line can be formed in the oper¬ 
ations table to meet a particular require¬ 
ment; an example is given in Table 
VI. 

This particular operation would be used 
if one residual were very large compared 
with the other. 

Equalization of Residuals. In any set of 
linear equations, if the initial residuals are 
Rob), R 0 (y), Ro(z) ... it is possible by ap¬ 


plying appropriate increments A+r, A n y, 
A n z ... to reduce these values to the values 
R n (x), R n (y), R n (z) . . . such that 

... =Z(0<A<1) 

Ro(x) Rob) 

It is then found that further increments 
(K/l-K)A n x, (K/l-K)A n y, (K/l-K)- 
A n z .. .will reduce all the residuals to zero. 

Application of the Method to 
Coupled Circuits 

For the general case of a £-mesh linear 
network, the steady state is described by 
the p equations 

Ei=*Znh +Znh + • • • f-ZipIp 
Ei— Z n Ti -\~Z<i 2 I 2 -f~ • • • ~\~ZopIp 


Ep — Z v \h AZ-pR’i + . . . +Z VP I P 
where 

E’s are impressed voltages (constant) 

Z’s are impedances (constant) 

Ps are mesh currents (unknown) 

The currents then can be found by 
solving this system of simultaneous linear 
algebraic equations. While the solution 
can be obtained by elimination, it is 
found that in most cases relaxation is 
simpler and quicker to carry out. 

Relaxation has the added advantages 
that 

1. It very quickly gives an approximate 
solution which may be all that is required. 

2. Starting from this approximate solution, 
a solution of any desired accuracy can be 
obtained. 

3. The progress of the solution can be 
checked easily by recalculating the residuals 
from the original equations. 

Its disadvantage is that all the currents 
must be determined simultaneously 
whether they are all of interest or not. 

Application of the Method to 
Characteristic Value Problems 

Systems of simultaneous equations 
dealing with characteristic value prob¬ 
lems, for example, vibrating engineering 
structures, can be solved by a combina¬ 
tion of the relaxation process as described, 
and of Rayleigh’s principle. This method 
has been discussed fully in a paper by Fox 4 


Table VI 




AR\ 

AR 2 


As = 1 . 

.... —3..... 

.... 1 


Ay = 1 . 

.... 1_ 

.... -3 

oi 

II 

R 

<1 

.... Ay = 1 . 

.... —5.... 

--1 
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Figure 1. Subdivision of 
range into finite intervals 


in which examples are given showing how 
certain devices can be used to shorten the 
process, and also how modes other than 
the lowest can be found. 

Solution of Ordinary Differential 
Equations 

Nature of the Problem 

Ordinary differential equations of the 
second order can be written in the form 



where z is a known function of x, y, and 
dy/dx. For a unique solution of the equa¬ 
tion two boundary conditions also must be 
specified. Sometimes the forms of the 
equation and boundary conditions are 
such that an analytical solution is impos¬ 
sible and then numerical methods must 
be used. Usually, either y and dy/dx are 
known at one end of the range of solution 
or y is known at each end. In the first of 
these cases, an initial-value problem, the 
problem is usually solved by step-by-step 
methods, while in the second case, a 
boundary-value problem, relaxation is a 
much better method. 

o—©—o 

Pa Po Pi 

Figure 2. Incomplete relaxation pattern 
(first step) 


Although differential equations require 
for their solution a continuous function in 
the given range, relaxation, like all other 
numerical methods, gives only numerical 
values at a series of points. These points 
are taken for convenience as equidistant 
points along the x-axis, and values at in¬ 
termediate points then can be found by 
graphical or algebraic interpolation. 

Relaxation cannot solve differential 
equations as such, but by replacing them 


by the corresponding finite-difference 
equations they can be converted into sys¬ 
tems of simultaneous linear algebraic 
equations which are amenable to the 
standard relaxation techniques, provided 
that the solution does not contain singular 
points within the range of the solution. 

Differential equations leading to char¬ 
acteristic values and characteristic solu¬ 
tions can be dealt with as mentioned pre¬ 
viously in the section entitled “Applica¬ 
tion of the method to characteristic 
value problems” after the differential 
equation has been replaced by a set of 
simultaneous linear equations. 


Finite Difference Approximation 
Consider, for example, the equation 




with the boundary conditions that y is 
known at each end of the range of x under 
consideration. This equation, with K 
being a constant, of course can be solved 
directly and the appropriate boundary 
conditions inserted, but it was chosen 
here as an example to demonstrate the 
method in a simple way. One divides 
this range of x into equal mesh lengths a. 
Taking any nodal point as the origin 
P 0 (x 0 ) where y has the value y Q , and ad¬ 
jacent points Pi(xo-fa) and P 3 (x 0 — a) 
where y = yi and y — y 3 as shown in Figure 
1, one can show from Taylor's expansion 
that the finite difference approximations 
are 



and 

/ d 2 y\ yx+y 3 —2y 0 
\dx 2 /o a 2 

The given equation therefore can be ap¬ 
proximated by 

yi(2-\-ak') —Fys(2 — uk ) —4yo 2a 2 zo — 0 

This is a linear algebraic equation, and 
since an equation of this form connects 


o—o— G> 

Pi 




Figure 3. Incom¬ 

jT a _ak\_ 

r\ 

plete relaxation pat¬ 

v 


tern (second step) 

Po 

Pi 



(ES> —CD—<CD 


Figure 4. Complete relaxation pattern 


©-©-© 

Figure 5. Relaxation pattern for differential 
equation y" — z(x) = 0 



Figure 6. Assumed initial values of residuals 



0 

6 

10 

20 
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Figure 7. Increments and resulting residuals 
after first relaxation step 



Figure 8. Relation pattern for block relaxation 


every set of three consecutive points, the 
given d i ff erential equation has been re¬ 
placed by a set of simultaneous linear 
algebraic equations, the number of equa¬ 
tions in the set being given by the number 
of points at which the solution is deter¬ 
mined. 

The residual R 0 at the point P 0 is given 
by 

| Ro = yi(.l+ak/2)+y a (l—ak/2)—2y a — a 2 za | 

(5) 


Basic Process of Relaxation 

Since all the equations 5 (for different 
points Po, Pi... P n ) are of the same form, 
every line of the operations table is of the 
same form.. It is best visualized not as a 
table but as a “relaxation pattern.” 

This relaxation pattern is a graphical 
representation of the effect on residuals 
at various points of unit increment at any 
given point. Patterns can be formulated 
for any residual formula, and they are the 
essential feature of the relaxation method. 
Once the pattern has been made, it can 
be applied to the numerical work in which 
all the residuals are reduced to negligible 
values. It is left to the discretion of the 
computer to decide at which points these 
patterns will be used most profitably. 
This makes relaxation a more interesting 
and more flexible computational tech¬ 
nique than the purely mechanical itera¬ 
tion methods. 

The character of the relaxation pattern 
depends on the form of the given differen¬ 
tial equation. It may be the same for all 
points in the region, or it may depend on 
the position of the origin. The method of 
deriving the pattern, however, is the same 
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Figure 12. Potential values, incre¬ 
ments, and residuals after first relaxa¬ 
tion step 


Figure 9. Mesh points surrounding point P„ 

in all cases. As an example, consider the 
formulation of the pattern corresponding 
to equation 5. 

It is seen from equation 5 that AP 0 = — 
2 for Ay 0 = + 1. This is represented 
graphically by the (incomplete) pattern, 
Figure 2, the increment Ay being written 
above the circle, the resulting change of 
residual inside the circle. However, the 


complete relaxation pattern, representing 
the change of all residuals due to a unit 
increment of y at P Q , is shown in Figure 
4. 

Consider now the simple example where 
K = 0. Here the relaxation pattern re¬ 
duces to that shown in Figure 5. Suppose 
the residuals are initially as shown in 
Figure 6. As in the case of simultaneous 
equations, one starts by picking out the 
largest residual 20. From the relaxation 



Figure 10 (left). Relaxa¬ 
tion pattern for Laplace’s 
equation 


Figure 1 3 (right). Potential 
values, increments, and re¬ 
siduals after further relaxa¬ 
tion 



increment Ay 0 does not only .change the re¬ 
sidual at the point P 0 but also the resid¬ 
uals Pi and P 3 at the neighboring points 
Pi and P 3 . The changes of. the residuals 
at these points due to the increment 
Ay 0 — 1 are obtained from equation 5 by 
changing the origin first to P x and then to 
P 3 . Taking the new origin at Pi, the 
value y 0 takes the place of y 3 in equation 5, 


pattern one sees that to make this zero 
one has to make an increment +10, so 
changing this residual by 10x(—2) and 
those on either side of it by 10x(+l). 
This operation is shown in Figure 7, the 
increments being written to the left of the 
vertical lines, the new values of the re¬ 
sidual to the right. (The ends of the 


range are given by fixed boundary condi¬ 
tions and so there can at no time be any 
residual there.) The process is continued 
until all the residuals may be deemed 
negligible, for the accuracy desired. It is 
also necessary to impose the condition 
that the total residual shall be small. 

Shortening Devices 

As before, overrelaxation is quite 
straightforward and very useful. 

Patterns for block relaxation are formed 
by superposing individual patterns and 
are seen to be of the form shown in Figure 
8 . 

Patterns for group relaxation can also 
be formed to suit any particular case; 
see the section on “Numerical Example” 
given later in the paper. 

Accuracy of the Method 

The difference equations can be solved 
to any degree of accuracy but the ultimate 
accuracy of the solution depends on the 
closeness of the approximation of these 
difference equations to the given differen¬ 
tial equation. Since terms of order a 3 
were neglected in this approximation, the 
solution can be accurate only when 
worked for values of a which are small 
enough for these terms to be truly neg¬ 
ligible. The usual practice is to work 
solutions for a— (# 2 — xi)/2n{n — 1,2,3 . . .) 
where a is the mesh length and x 2 —xi the 
complete range of solution, with n in¬ 
creasing until two subsequent solutions 
differ by a negligible, amount. Even 
though initial solutions worked for a large 
value of a may be quite far from the final, 
correct solution, they are useful and time¬ 
saving in giving a good initial guess when 
the interval is halved, and the conver¬ 
gence of the solution can be seen. 

As an example, the equation 

d 2 y 

--;+»(1 —3C)10 4 = 0 
dx- 


and hence the increment Ay 0 produces in 
the residual Pi the change AP X = (1—a/a/2) 

Ay 0 . This is shown graphically in Figure 
3 for a unit increment Ay 0 = 1. Similarly 
one finds that the change in the residual 
P 3 is AP 3 = (1+ ah/ 2)Ay 0 . The resulting 

Figure 11 (left). Initial 
potential values and re- 
60 siduals for solution of 
Laplace's equation 
SO 

45 Figure 14 (right). Sub¬ 
division of relaxation 

40 ne ts 
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Tabic VII 




Extra- 

Analyt- 



polated 

ical 

n 

12 3 

4 Value 

Solution 

y. ■ 

. .313...273...264, . 

.261...260.. 

.260.4 



Figure 15. Relaxation pattern at point A in 
Figure 14 


GKm© 


Figure 16. Relaxation 
pattern at point B in 
Figure 14 



© - © 
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Figure 17. Relaxa- 
tion pattern at point 
E in Figure 14 



Figure 18. Relaxation pattern at point C in 
Figure 14 


© - © 


Figure 19. Relaxation 
pattern at point D in 
Figure 14 



Figure 20 (left). 
Relaxation pa tternat 
point F in Figure 14 


Figure 22 (right). 
Fictitious boundary 
point to represent 
prescribed boundary 
gradient 
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with y = 0 at x = 0 and y = 0 at x — 1, which 
has the analytical solution y— 260.4 at 
x= 1/2, has been solved for various n and 
the results are shown in Table VII. 

Often, the solution approaches upon re¬ 
peated reduction of the interval size a the 
correct value in an oscillatory manner and 
extrapolation must be carried out with 
care as was pointed out by v Salvadori. s 

Solution of Poisson’s and Laplace’s 
Equations in Two Dimensions 

Nature of the Problem 

Partial differential equations may also 
be of initial-value or boundary-value 
types. As in the case of ordinary differen¬ 
tial equations it is the boundary-value 
types which are solved most conveniently 
numerically by relaxation. v Such prob¬ 
lems arise frequently in electrical engi¬ 
neering. For instance, in electric ma¬ 
chines, the shape of the boundary and the 
field values on the boundary are known, 
and the field configuration inside the re¬ 
gion is required. Here, electric and mag¬ 
netic potentials satisfy Poisson’s or La¬ 
place’s equations. 

The first step in solving these equations 
is to replace them by their finite difference 
approximations, that is, a set of simul¬ 
taneous linear algebraic equations holding 
at various points at which the solution is 
to be obtained. The most convenient 
points to take are the nodes of a symmet¬ 
ric and systematic network; in practice, 
square meshes generally are used, though 
triangular and hexagonal meshes are geo¬ 
metrically possible. 


may be written in finite difference form as 

4 

E —4Fo = a%o 

n = l 

a being the mesh size (see Figure 9) 
whence the residual R 0 at the point P 0 is 
given by 

4 

Ro = E Vn-iVo-aho (6) 

n=i 

For Laplace’s equation z~ 0, and the 
equation reduces to 

Ao = E !/ n-4V, (7) 

n = 1 

Basic Process of Relaxation 

By arguments similar to those in the 
preceding section of the same name, it is 
found that the appropriate relaxation pat¬ 
tern in this case is that shown in Figure 
10 . 

In solving Laplace’s equation, for given 
boundary values, in the region shown in 
Figure 11, V values first are guessed and 
entered to the left of the nodes, the cal¬ 
culated residuals being entered to the 
right. 

The largest residual then is relaxed as 
shown in Figure 12, the process continued 
as in Figure 13 and further continued until 

1. All the residuals are small. 

2. The total residual is small. 

3. The total residual in any area is small. 

Finally, all increments A V are summed for 
each node. 

Shortening Devices 


Finite Difference Approximation 


Poisson’s equation 
dx 2 ""^dy 5 


• =z(xy ) 



Any device which will increase the 
speed of the method is of great importance 


Figure 21 (left). Unequal network star for 
curved boundary 

Figure 23 (below). Section through instru¬ 
ment magnet 
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Figure 24. Relaxation net for magnet, Figure 23, showing computed field values and residuals 


and again overrelaxation proves very use¬ 
ful. 

Local''bumps” in the field are smoothed 
out readily by point relaxation but most 
residuals left may be of the same sign, 
giving a large value for the total residual. 
The only way to change the total re¬ 
sidual in a given area is to move residuals 
away over the boundary. In a large 
field most points axe a considerable dis¬ 


tance away from the boundary and it 
takes a great number of point relaxations 
to shift residuals as far as the boundary. 
Block relaxation in which a number of 
points are treated together and hence any 
point becomes relatively nearer the 
boundary, therefore proves very useful in 
reducing the total residual to zero. The 
patterns for this are formed by superpos¬ 
ing individual patterns, and easily can be 



Figure 25. Refinement of relaxation net for 
magnet. Figure 23 


made to suit any shape of area. It is 
found that if the area to which unit incre¬ 
ments are made is enclosed by a contour, 
then one unit of residual goes down each 
mesh line cut by the contour. If, there¬ 
fore, the contour encloses a total residual 
R, and cuts N lines, by making incre¬ 
ments of magnitude R/N, the residual in¬ 
side the contour will be reduced to zero. 

Accuracy of the Method 

In the derivation of the finite difference 
formulae, terms of order a 4 and higher 
were neglected. To overcome this limi¬ 
tation to the accuracy, two procedures are 
possible. One is to make a so small that 
these terms really are negligible, and the 
other is to include these terms by a 
slightly modified procedure. These two 
alternatives now will be discussed in 
turn. 

Graded Networks. The choice of mesh 
length in solving field equations is very 
important; too large a mesh length a 
gives a result which is inaccurate due to a 
finite mesh size error, while too small an a 
requires much unnecessary computation. 
Adopting the same criterion as before, 
that is, that the solution does not change 
appreciably on reducing the mesh size, it 
frequently is found that this condition is 
satisfied in some parts of the field but not 
in others, for example, in the region of 
sharp corners where the field varies very 
rapidly. To minimize labor, it is possible 
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Figure 26. Magnetic flux 
lines and equipotentials of 
magnetomotive force for mag¬ 
net, Figure 23, as finally com¬ 
puted 



to use a fine mesh only where necessary, 
but one then has the problem of joining 
the two types of net. 

In Figure 14 a net of mesh of size a is 
shown, which is to be halved in the right- 
hand region. The procedure is as fol¬ 
lows: 

1. Draw in the extra mesh lines required. 

2. Guess values at the new nodes so formed. 

3. Calculate the residuals at these points. 

To be able to calculate the residual at a 
point, the value of the function has to be 
known at that point and also at the sym¬ 
metrically surrounding points. Some of 
the nodes just formed do not have sym¬ 
metrically surrounding points and so pre¬ 
sent a special problem. Referring to 
Figure 14: 


B 



A 



C 

D 



Figure 27. Section through cylindrical air- 
gap choke 


Columns 1, 2, 3, and to the left—these are 
all in the coarse mesh. 

Columns 6, 7, and to the right—these are 
all in the fine mesh. 

Column 5—alternate points can be regarded 
as being in the coarse mesh, but the others 


have no point to their left. Produce the 
fine net one step further at these places to 
get the points shown in column 4. 

Column 4—these points are now surrounded 
by points distance a/y/2 away. The re¬ 
siduals are now calculated from equations 6 
or 7 but with a replaced by a/V2. (This 



Figure 28. Magnetic flux lines and equipotentials of magnetomotive force for air-gap choke. 

Figure 27, as computed by relaxation 
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is permissible since V 2 remains invariant 
under rotation of axes.) 

The residuals at all points now can be 
calculated. In the region of the join, 
however, modified relaxation patterns 
also are required. Points A, B, C, D, E, 
and F in Figure 14 are typical points in 
the region, 

Consider the point A. Changing the 
value there by -f-1 means changing its re¬ 
sidual by —4 units. Now consider which 
other points are affected by this change. 
The value at A is used in the calculation 
of the residuals at all the points sur¬ 
rounding it on the large mesh, and also 
at the points above and below it on the 
diagonal mesh. Thus, the complete pat¬ 
tern at A is that shown in Figure 15. 

Similarly, the patterns at points B, C, 
D, E, and F are found to be those of 
Figures 16 to 20. 

Now these patterns can be used to re¬ 
lax in both regions and across the join. 

It is sometimes possible to dispense 
with the diagonal mesh, but then relaxa¬ 
tion can only be carried out in each region 
separately and residuals are liable to build 
up on the boundary. 

Inclusion of Higher Order Terms. Fox 9 
has shown that by taking account of 
higher order terms in the finite difference 


Figure 21. In this case there is no means 
of knowing V\ (at point Pi), but Vb, the 
potential at the boundary point B, is 
known as part of the boundary condition. 
Also, the length pa can be calculated. 

From a Taylor expansion first with 
x — a, and then x=—a, equation 8 for 
the residual R Q at P 0 is obtained. 


known. For all points on the boundary, 
Vi would have a fictitious value, and 
hence the standard equation cannot be 
applied directly. However, one may say 



hence V\ may be found in terms of known 

( 8 ) 



In the case where two lines are cut by 
the boundary at distances pa and qa from 
P 0 , the corresponding residual equation 
9 is 



In obtaining these equations it was 
necessary to neglect terms of order a s , 
making them inherently less accurate 
than the standard equations 6 or 7. 

Boundary Conditions 

A slightly modified procedure is neces¬ 
sary if not the function, but its normal de¬ 
rivative is known at the boundary. 

Consider first a rectilinear boundary as 
shown in Figure 22, with A the direction 
normal to the boundary, &F/5.S being 


values and substituted in the standard 
equation. 

When the boundary is curved, two of 
the required values may be fictitious, and 
the normal to the boundary is not in 
general parallel to the mesh lines. This 
makes it a little more difficult to express 
the fictitious values in terms of known 
ones, but it is still quite straightforward. 
This method of attack is due to Fox. 10 An 
alternative approach due to v Southwell 


approximation greater accuracy can be 
obtained with a relatively large value of 
the mesh length a. 

Since the work involved increases so 
rapidly with a, this method of computa 
tion is often very useful. The calculation 
of the required differences, however, may 
be a little tedious. 

Curved Boundaries 

If the region considered has curved 
boundaries, it is impossible to draw a 
square mesh to fit it exactly. There is 
then the problem of "unequal stars” as in 
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Figure 29. Magnetic field inside dynamo- 
electric machine computed by relaxation 
(Motz and Worthy) 
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Numerical Examples 



Figure 30 (left). Mesh points 
surrounding point 0, used in finite 
difference approximation of 
“mixed" derivatives 


Figure 32 (right). Mesh points 
used in finite difference approxi¬ 
mation of “quasi-plane potential 
equation" 



and Vaisey 11 is based on the membrane 
analogy. This leads to different equa¬ 
tions, but Fox has shown that the agree¬ 
ment between solutions worked by either 
method is quite good. 

Numerical Example 

Consider as an example the bar magnet 
with pole pieces as shown in Figure 23. 
This is similar to the type used in measur¬ 
ing instruments, and so is of frequent oc¬ 
currence. The problem is really a 3 -di¬ 
mensional one, but if the width is reason¬ 
ably large compared with the other di¬ 
mensions, a section through a mid-plane 
can be regarded as experiencing no trans¬ 
verse end-effects, and can be treated as a 
2 -dimensional case. This means the 
governing equation has the form 

3 2 » dH' 
hx 2 3y 2 

where 11= — grad being the mag¬ 
netic potential and H the field strength. 
The shape of the boundary is specified 
accurately, but the values of II — — grad 
T on the boundary are not known ex¬ 
actly. As a first approximation it is as¬ 
sumed that there is a linear drop of mag¬ 
netic potential T r along the magnet and 
that the iron has infinite permeability. 

Working on a net of mesh length one- 
quarter the gap width, the residuals were 
found to be predominantly negative giv¬ 
ing a large value for the sum of the re¬ 
siduals. The following device was found 
particularly useful in reducing this total 
residual. The total residual in each 
column was calculated and, the number 
of points in each column being known, it 
was possible to make patterns for relaxing 
the whole column at once. In this way 
the distribution shown in Figure 24 was 
obtained. Smaller meshes were used 
then, as shown in Figure 25, and since the 
residuals on the boundaries of each mesh 
were found to change very little, diagonal 
nets for joining the meshes were unneces¬ 


sary. Equipotentials and flux lines are 
shown in Figure 26. 

Solution of 3-Dimensional Partial 
Differential Equations 

Although the finite difference equations 
can be formulated easily for 3-dimen¬ 
sional problems, the number of points in¬ 
volved and the difficulties in recording 
make the method impracticable in gen¬ 
eral. 

Several special cases' of 3-dimensional 
partial differential equations, however, 
can be treated in this way. Tranter , 12 for 
instance, has shown how to reduce the 
number of independent variables in a 
given partial differential equation by one 
by using Fourier transforms, thus reduc¬ 
ing 3 -dimensional equations to 2 -dimen¬ 
sional form. Also, the important class of 
equations having rotational symmetry 
can be written in essentially 2 -dimen¬ 
sional form. 

Axially Symmetric Fields 
Laplace’s equation in this case is 

3 2 F 3 2 F 1 57 
3 r- 3k 2 dr 


whence we have the equation for the re¬ 
sidual 



( 10 ) 


where?’ 7 ^ 0 . As 

lim /I 3 F\ dW 
r— *-0\r 3 r ) 3 r 2 

One finds for the residual at a point P a on 
the axis 

| Rq= Vi 4 -F 3 T 4 F 4 6 F 0 (11) 

where r=0. 



Figure 31 (left). Relaxation pattern for equation 

3 x 2 Oy 2 bx3y 


Figure 33 (right). Relaxation pattern for 
"quasi-plane potential equation" 



Cylindrical Air-Gap Choke. Consider 
the air-gap choke shown in Figure 27. 
Assuming as boundary conditions a linear 
drop of magnetic potential drop along 
AB and CD, and constant potential on the 
rest of the boundary, relaxation was per¬ 
formed on successively finer nets, as de¬ 
scribed in the previous section on "Graded 
Networks.” The equipotentials then were 
plotted and these with the flux lines are 
shown in Figure 28. 

Dynamo-Electric Machines. As shown 
in a paper by Motz and Worthy 13 * the 
magnetic field inside a dynamo-electric 
machine can be calculated by this method. 
Their result for a 4-pole machine with a 
slotted armature and distributed wind¬ 
ings is shown in Figure 29. 

Electron Lenses. Most electron lenses 
have axial symmetry and hence the 
governing equation can be written in its 
2 -dimensional form. Motz and Klanfer 14 
have used relaxation to calculate the 
field distribution in an electrostatic lens, 
consisting of two semi-infinite cylinders 
placed end to end and separated by a 
finite gap. Hesse 15 recently applied the 
method to the calculation of the field dis¬ 
tribution inside a magnetic electron lens, 
taking account of saturation of the pole 
pieces. 

Solution of Special Partial 

Differential Equations 

Types of Equations 

Difference equations can be formed to 
replace most differential equations oc¬ 
curring in electrical engineering, and the 
appropriate relaxation patterns found by 
similar methods to those already de¬ 
scribed. If, however, the equation is very 
complicated, so is the difference equation, 
and the relaxation pattern may involve a 
great number of points. It is still possible 
to work with such patterns as South- 
well 16 has shown, but the manipulation 
then becomes rather involved and tedi¬ 
ous. In these cases, however, solution by 
any other method is practically impossible 
and so the labor involved in relaxation is 
justified. 

Equations Containing Mixed 

Derivatives 

The finite difference approximations to 
mixed derivatives can be found quite 
easily. With the arrangement of points 
as shown in Figure 30, a Taylor expansion 

* The author wishes to thank Dr. H. Motz and the 
Institution of Electrical Engineers, London, Eng¬ 
land, for permission to reproduce Figure 29 from 
one of his papers. 
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leads to this approximation 
/ d*V\ 

4a2 \ MM “ Va-Vb + Vc-V d 

\ ox uyJ y 

The use of this can best be seen from 
this example 

527 527 527 

- 4 .-j_2 --= 0 

bx 2 by 2 bxby 

This can be written approximately 

Vi+V s -2V 0 , F 2 +F 4 —2F 0 , 

_L. ■ -J- - 

a 2 a 2 

Va-Vb + Vc-Vd _ 0 


( 12 ) 

The relaxation pattern corresponding to 
equation 12 is shown in Figure 31. 


4a 2 

whence the residual equation is 
(i ? 0 = Z F t ■—4 Fo+Vs ( Va -Vb + Vc - Vd) 


Quasi-plane Potential Equation 

This is the name given to the equa¬ 
tion 



+Z = 0 


which reduces to Poisson’s equation, (that 
is, the plane harmonic equation) when e 
is constant. It is, for example, the equa¬ 
tion for the electric field F in a region of 
nonuniform dielectric constant e. 

With the arrangement of points as in 
Figure 32, and denoting the value of € be¬ 
tween P a and Pi by e a , and so forth, the 
finite difference approximation to this 
equation is found to be 


~ [«i( Fj - Fo) - eg( Fo - F 3 ) +«( V 2 - Fo) - 

«4(7o-F«)]+Z-0 


whence the residual formula 13 is ob¬ 
tained 


4 4 

i?o — e nVn ET.2 e m -f-U 2 Z (13) 

n = i 77i sa 1 

The corresponding relaxation pattern is 
shown in Figure 33. 


Resonator and Oscillation 
Problems 

.... In most oscillation' problems (for ex¬ 
ample, the field distribution inside a res¬ 
onator) the governing equation is of the 
form 

V a F+£ 2 F=0 (14) 

The boundary conditions in wave¬ 
guide and resonator problems are such 
that this equation has nontrivial solutions 
only for certain discrete values of k 2 . 
Solving the equations therefore involves a 


determination of the allowed parameters 
k 2 (the “characteristic values”) as well as 
of the corresponding function V. They 
can be solved by a combination of relaxa¬ 
tion and iteration, making use of Ray¬ 
leigh’s principle which gives for the funda¬ 
mental mode 


Sfv- V 2 Vdxdy 
•dxdy 


(15) 


The procedure is to assume values for 
F and with these values calculate k 2 from 
equation 15. This value is substituted in 
equation 14 which then can be relaxed in 
the usual way. From this a corrected 
value of F is found which is used in equa¬ 
tion 15 to obtain a corrected value of k 2 , 
and so forth. The process is continued 
until corrections are negligible, 3 or 4 
cycles usually being sufficient. 

Allen, Fox, Motz and Southwell 17 ap¬ 
plied this treatment to a 2 -dimensional 
section through a klystron tube. Motz 18 
extended the method to finding the field 
distribution and resonance frequency of 
cavity resonators, improving the accuracy 
by taking account of the analytic be¬ 
haviour of fields in the region of sharp 
corners. 


Conclusion 

Southwell and his co-workers, who de¬ 
veloped the method of relaxation, were 
interested mainly in applying it to prob¬ 
lems in mechanical engineering. How¬ 
ever, as this account has shown, the prin¬ 
ciples and techniques established by them 
are easily applicable to many problems 
in electrical engineering. The most im¬ 
portant and most frequent use of the 
method is in the solution of the electro¬ 
magnetic field problems described by 
Laplace’s, Poisson’s, and Maxwell’s equa¬ 
tions. Its great advantages as compared 
with analytical methods are: 

1. It is simple to apply. 

2. The process is essentially the same for 
all types of equations. 

3. It can be checked easily during the 
progress of the work. 

4. It can be applied to problems for which 
no analytical solution can be obtained. 

5. It can mostly be made as accurate as 
a given problem requires it to be. 

The chief disadvantage is that the 
solution obtained is a numerical one and 
applies only to the one case for which it 
has been worked. A variation of any 
parameter necessitates a complete re¬ 
calculation. The solution for the first 
case, however, often gives some guidance 
in the solution of the next, thus reducing 


the effort required for the whole range of 
solutions. 
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Notes on the Design of Eccles-Jordan 

Flip-Flops 

MORRIS RUBINOFF 

ASSOCIATE MEMBER AIEE 


Synopsis: The design of a prototype 
Eccles-Jordan flip-flop such as used in large- 
scale digital computers is discussed, with 
particular attention to the role played by 
component variability. A graphical design 
technique is described; the technique should 
prove useful in other switching-circuit de¬ 
sign problems where a large number of 
identical direct-coupled circuits are em¬ 
ployed. 

A preliminary analysis is carried through 
for the rise time of the on-going control grid. 
The results may be used to indicate the in¬ 
fluence of transpose capacitance on speed of 
switching of flip-flops. 

I T IS common knowledge that the bi¬ 
stable vacuum-tube device referred to 
as the Eccles-Jordan flip-flop was first 
described by Eccles and Jordan in 1919. 1 
It is perhaps not so well known that the 
device described by them is not identical 
with the current flip-flop in that they 
used batteries in the place of transpose re¬ 
sistors and thus dispensed with resistor 
divider networks. 

The flip-flop now in common use is 
shown in Figure 1. It is seen to consist of 
the following; two divider networks each 
containing a plate resistor R h a transpose 
resistor R 2 usually shunted by a trans¬ 
pose capacitor C 2 , and a grid resistor i? 3 ; 
two tubes, either triodes or pentodes; 
and two supply voltages, one positive 
and the other negative with respect to the 
grounded cathodes. Variations of this 
configuration are in use; note in par¬ 
ticular that it is not necessary that the 
two halves of the flip-flop be nominally 
identical. However, this paper will re¬ 
strict consideration to the illustrated cir¬ 
cuit. 

Figure 1 illustrates one stable state in 
which the left tube is completely cut off 
while the grid voltage of the right tube is 
slightly above ground. A. T. Starr 2 
showed that stability can be attained even 
though the “off” tube is not cut off and 
the “conducting” tube has negative grid. 
In fact, it is easy to show that stability 
requires only that, in the quiescent state, 
the loop gain around the two tubes be 
less than unity. On the other hand, it is 
advisable to protect against malfunction¬ 
ing due to noise and, therefore, to reduce 
the loop gain by ensuring both cutoff for 
one tube and grid conduction for the 
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other. In any case, it should be noted 
that a variation in any parameter which 
increases the grid voltage of the cutoff 
tube or decreases the grid voltage (or 
current) of the conducting tube leads the 
flip-flop closer to instability; this cri¬ 
terion will be used later to determine the 
influence of parameter variations on flip- 
flop stability. 

Scope of Paper 

The design of a single Eccles-Jordan 
flip-flop for laboratory use is relatively 
simple and straightforward. On the con¬ 
trary, the design of a prototype flip-flop 
for use, for example, in a large-scale 
digital computer requires more care be¬ 
cause of the large number of flip-flops in¬ 
volved. First, the flip-flop must fit into 
the framework set by other (computer) 
circuits, which implies specified output 
voltage levels and switching times. 
Secondly, the flip-flop must exhibit the 
desired characteristics in spite of com¬ 
ponent variations, both initial and during 
use in the machine. In this respect, at¬ 
tention must be focussed on resistor drift, 
tube aging, and fluctuations in supply 
voltages. 

The general problem of flip-flop design 
with all its ramifications is much too com¬ 
plex to permit of its discussion here. 
Hence, the following assumptions are 
made in order to restrict and to define 
the scope of the paper: 

1 . All resistors will be so selected that 
after long service in the computer they will 
not drift outside a given fractional tolerance 
p from nominal resistance. 

2. All supply voltages can be held at all 
times to within a given fractional tolerance 
<r. 

3. Knowledge of tube characteristics for 
the complete life of every tube (or tube re¬ 
jection criteria) are such that: 1. plate 
current i v for zero grid voltage % and speci¬ 
fied plate voltage e v will always exceed a 
known value, and 2. grid bias for given e v 
need never exceed a known value to ensure 
cutoff. 

4. Satisfactory operation is ensured if: 1. 
tube cutoff is ensured, 2. the grid of the 
conducting tube reaches at least zero bias, 
3 . zero grid bias implies zero grid current. 
(The grid current condition has been in¬ 
cluded only to simplify exposition. The 


graphical design procedure described in this 
paper is readily modified to take account of 
grid current at zero grid bias.) 

5. Only the grounded-cathode flip-flop is 
to be considered, and no inductors are to be 
used. 

Notation 

The notation used in this paper is illus¬ 
trated in Figure 1 . E x and E» refer to the 
positive and negative supply voltages, 
while E and G refer to plate voltage and 
grid voltage respectively. Algebraic signs 
have been appended so that all symbols 
are positive quantities. The subscripts 
C and NC refer to the C and NC divider 
networks connected to the conducting 
and nonconducting anodes respectively. 
Finally, the transpose capacitance C 2 is 
assumed to be a variable parameter 
available to the designer, but the stray 
capacitances C] and C 3 in shunt with Ri 
and Ri respectively are fixed by tube 
type and physical layout. In this regard 
it is worth mentioning that C 3 does not in¬ 
clude Miller capacitance; this is cer¬ 
tainly justifiable for the off tube, which 
is the only case treated here, but can be 
largely justified in any case. 

Derivation of the “Worst” 
Configuration 

It is proposed now to determine the 
“worst” possible configuration of the pa¬ 
rameters, that is, that configuration most 
nearly approaching instability as de¬ 
fined above. A graphical approach is to 
be employed, partly because it has proved 
to be simpler than the analytic method, 
but mainly because the graphical method 
is most illuminating and immediately ap¬ 
plicable to design of more-sophisticated 
flip-flops, of univibrators, and of all cir¬ 
cuits using resistance divider networks. 

Consider first the NC divider illus- 
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trated in Figure 2. Note in the upper 
graph the “minimum” plate character¬ 
istic for e 0 = 0, item 3 in the section ‘ ‘Scope 
of Paper.” The divider current i NC is 
drawn as if it were negative plate cur¬ 
rent, the reason to be given later. The 
solid lines labelled R lt R 2 , and i? 3 are seen 
to be “load lines” for the respective re¬ 
sistors, that is, their slopes are the re¬ 
ciprocals of the resistances represented. 
The broken lines show the effect of more 
negative-E* supply (larger E 2 ). Clearly, 
the resultant change in G tends toward 
flip-flop instability since the originally 
zero grid goes negative. 

Similarly, the lower graph of Figure 2 
shows that an increase in R% in the NC 
divider tends toward instability. By 
graphically studying the effect of varying 
each parameter, one at a time, its in¬ 
fluence on stability is rapidly ascertained. 

The C divider can be studied graphi¬ 
cally in similar fashion, as shown in Figure 
3. Here again the solid lines are load 
lines. Note in particular that the load 
line for Ri fits properly onto this graph 
only because the divider current i c has 
been treated as negative i v . The graph 
shows that decreasing Ri tends to insta¬ 
bility because it tends to raise the grid 
■voltage of the nonconducting tube. Ac¬ 
tually, it is possible to derive more in¬ 
formation from this graph by noting that 
the change AG is due to increased i c which 
is in turn due to increase in E c . It fol¬ 
lows that decreasing Ri will tend to in¬ 
stability whenever the plate character¬ 
istic has a positive finite slope at the point 
(E c , i p ), a condition generally realized. 

The effect of variation of each param¬ 
eter of the C divider can be studied 
graphically to yield information regarding 
tendency toward instability as well as 
associated conditions on the minimum 
plate characteristic. Although this is not 


to be carried through in this paper, the 
results can be summarized as follows: 

Lemma I: 

(a) . For the NC divider, the most nearly 
unstable configuration is obtained with 
Ei(min), i?i(max), i? 2 (max), Rz( min), and 
jE 2 (max): 

(b) . For the C divider, the most nearly 
unstable configuration is obtained with 
.Ei(max), 7?i(min), i? 2 (min), i? 3 (max), and 
E 2 (min). 

Study of the configurations of lemma 
I will show that they do not result in ex¬ 
treme values of E; nevertheless, the 
notation E(max), E(min) will be used 
for this case. 

One might object that 1(a) and (b) are 
inconsistent since Ei and £ 2 are the same 


for both dividers and cannot be both 
maximum and minimum simultaneously. 
However, the purpose is to design a flip- 
flop stable over long periods of time. 
Hence, it is necessary to apply both parts 
of the lemma simultaneously in the initial 
design. Each part was independently 
derived for the same reason. 

Some Design Lemmas 

Lemma I can be applied to derive an¬ 
alytic expressions for the network re¬ 
sistances. Thus from 1(a) it follows that 

,_ Fi(min) _ E 2 (max) 

tN ° i?i(max)+I? 2 (max) E 3 (min) 

while from 1(b) is obtained 



Figure 1. Grounded-cathode Eccles-Jordan Hip-flop. (Super- 
cript and subscript bars denote maximum and minimum respectively 


Figure 3. Influence of C-divider parameters on stability 
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Figure 4 (left). Determination 
of optimum Ec(max). (Super¬ 
script and subscript bars denote 
maximum and minimum re¬ 
spectively) 


Figure 6 (right). Initial con¬ 
ditions for rising grid voltage 



when a larger value of Ei is chosen. 
Lemma IV; 


tc = 


£ c (max)+G c (min) 


jRs(ttiin) 


£ 2 (min) — G c (min) 


ip + ic 


Rs( max) 

£x(max)—£ c ( max) 
£x(tnin) 


( 2 ) 


(3) 


These equations can be transferred by 
straightforward algebraic manipulation 
to 


_ 

E+G 


--(1-g) 


= _Ei(max)_ 
£i(max)-£ 

(1 —p)( 1 — <r) 
= (l+p)(l + <r) 

£=£ e (max) 


52 = 1 _^2( 1 _g) 

G 

^ £ 2 (min) 

G=£? c (min) 


(£i(tnax) —£)(1 —g) 

e z 


( 6 ) 


ij,£i(min)==£x(rnax)5 z +G 


(4) 


. „ , . . E+G 

tp £2 (min) = —-X 

1 ~£ 

1 


/3 2 


E+G 


-a -g) 


(£i(max)- £)(1 — g) 
^£ 3 (max) = £j(min) X 


(5) 


The following lemmas may be seen to 
follow almost immediately from equations 
4, 5, and 6 : 

Lemma II: 

Ri, Rz, and £ 3 are all inversely proportional 
to i v ; that is, all other things equal, the 
resistances are minimized by choosing a 
tube with maximum perveance. (This 
lemma is modified slightly when 0 for 

«j=0.) 



If £ 2 be chosen appreciably larger than G so 
that g is small, then £1 and £2 are almost 
independent of £ 2 . £3 is approximately 

proportional to £ 2 . 

Lemma IV(a ): 

A reasonable choice for £2 is £ 2 = 10 G c (min). 
Lemma V: 

All three resistors increased as either resist¬ 
ance tolerance p increases or supply toler¬ 
ance <r increases or both, for then 0 and ($ 
decrease and 0 increases. 


Since neither £ 2 nor £3 can be negative, 
it is necessary for the denominator of the 
first term inside the bracket of equations 
5 and 6 to be positive. That is 


E+G 


(£x(max)-£)(l-g) 


>0 


(7) 


and therefore 
Lemma VI: 

The minimum permissible value of £1 is 
limited by the inequality 


£x(max)> 


E+G 1 

1 ~g ' 0* 


Finally it is proposed to derive a cri¬ 
terion for choice of £ c (max). The cri¬ 
terion will be based on the assumption, 
believed to be approximately true, that 
flip-flop speed is inversely proportional to 
£ 1 , that is, that Ri should be minimized. 
By lemma III, this occurs when Ei is 
chosen as small as possible. Applying 
lemma VI to equation 4, one obtains 


fp£i(min) > 


(E+G)S Z 
(1 -gW 


+ G 



( 8 ) 


It is helpful to introduce the resistance 
R\* given by 



(9) 


Figure 5. Graphical flip-flop design. (Superscript and subscript bars denote maximum and 

minimum respectively) 


Thus, £ 1 * is a convenient measure of 

£1. 
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Consider first a flip-flop designed 
around a pair of pentodes so that G is 
fixed by the screen voltage. The upper 
graph of Figure 4 illustrates a construc¬ 
tion for obtaining Ri* as defined by equa¬ 
tion 9. The procedure consists of draw¬ 
ing a line such as AC from the point 
( —G/(5 2 ),0) to intersect the minimum 
plate characteristic at C; if E be chosen 
to equal the plate voltage at C, the slope 
of the line A C is inversely proportional to 
Ri*. Clearly, Ri* is minimized when the 
line is drawn tangent to the character¬ 
istic, as indicated by the solid line AB 
(labelled i?t*(min)). The point of tan- 
gency determines the optimum value of E. 
Since 0 < S < 1 , Ri* can never be decreased 
below the value corresponding to the 
dashed line passing through (—G c (min), 
0) and tangent at F. Clearly, E should be 
chosen somewhere on the “knee” of the 
curve. 

In the case of a triode flip-flop, the same 
reasoning leads to the result that E be 
chosen as large as possible. The upper 
limit is set by the rated anode power 
dissipation W. Hence, E should be 
chosen at the “notch point” N at which 
the minimum plate characteristic inter¬ 
sects the “rated power” curve PF—con¬ 
stant as shown in the lower graph of 
Figure 4. The broken line IN indicates 
the corresponding i?!*(min). 

It should be recalled, however, that in 
the case of a triode, G is roughly propor¬ 
tional to Ei which is in turn roughly pro¬ 
portional to E. Hence does not in¬ 
crease with decreasing E as rapidly as 
might otherwise be anticipated. The 
solid line JK shows the effect of reducing 
all voltages by a factor of two. By draw¬ 
ing a line from the origin through the new 


design point K to intersect LN at M and 
by applying the lemmas of similar tri¬ 
angles, it can be shown that the new de¬ 
sign point results in an increase in Ri* in 
the ratio of LN to LM. In practice, the 
design point must be chosen well below 
the notch point because of other con¬ 
siderations, such as resistor power dis¬ 
sipation and magnitude of supply volt¬ 
ages. Consequently, this construction is 
a useful guide to the resultant increase in 
Ri. 

The preceding comments may be sum¬ 
marized as follows: 

Lemma VII: 

(a) For a pentode flip-flop, EJ max) should 
be chosen on the knee of the minimum plate 
characteristic. 

(b) For a triode flip-flop, £c(max) should 
be chosen as large as other considerations 
permit. 

Graphical Design Procedure 

The analytic expressions for the net¬ 
work resistances given by equations 4, 5, 
and 6 of the section entitled “Some De¬ 
sign Lemmas” contain the following seven 
parameters, E h En, E, G, i v , p, cr. The 
magnitudes of these parameters may be 
assigned arbitrarily by the designer, sub¬ 
ject only to the condition on Ei stated in 
lemma VI. The assignment of values 
for p and <r is dependent upon knowledge 
of the characteristics of the resistors and 
of the power supplies to be used, while 
the selection of G and i P depends upon a 
knowledge of the characteristics of exist¬ 
ing vacuum tube types; clearly, engi¬ 
neering judgment must enter into the 
selection of values for these. At the same 
time, lemmas IV, V and VII serve as use¬ 


ful guides which restrict the ranges of 
p, a, Ei and E. 

The remarks of the preceding para¬ 
graph suggest that all seven parameters 
may be chosen in advance and that the 
analytic expressions of equations 4, 5, and 
6 may be used directly for determining 
the values of R\, Ri, and Where flip- 
flop application is not overly restricted 
by side conditions, this is probably true. 
However, it has been the author’s ex¬ 
perience in computer design that side 
conditions force the engineer to try many 
possible flip-flop designs. In this event 
the computation of the analytic expres¬ 
sions is time consuming and the graphical 
approach now to be described is much 
faster. Moreover, as we have seen, the 
graphical method immediately illustrates 
the influence of varying each parameter. 

The design procedure is based directly 
upon lemma I and the assumptions speci¬ 
fied in the section entitled “Scope of the 
Paper.” The procedure is illustrated in 
Figure 5. 

1. Set up ip, e p axes; the voltage scale is 
known but the current scale is unknown 
initially. 

2. Draw a line parallel to and below the e p 
axis to represent i c , the C-divider current. 
(It is convenient to make i c —10 divisions.) 

3. Locate —E 2 (min), — (?c(min), Ec(m ax). 
Optimum choices for Ei and E are based on 
lemmas IV(a) and VII; — Gc(min) is the 
guaranteed grid cutoff bias. Draw load 
lines for l? 2 (min) and I? 3 (max). 

4. Draw the load line for As(min) in the 
NC divider. This line passes through 
( —E 2 (max), 0) and (— Gc(min) — [£ 2 (max) 
— E 2 { min)], [1+p/l —p ]). Produce this line 
to cut the current axis; the intersection is at 
i N c, the current through all three NC- 
divider resistors. 

5. Draw the load line for i? 2 (max) in the 
NC divider; it passes through (0,0) and 
[JS c (max)+G c (min), i c (l-p/l+p)]. Pro¬ 
duce the line to i NC ; the corresponding 
voltage is Ena • 

6 . Draw the load line for i?i(max) in the 
NC divider; it passes through {E NC { min), 
0 ) and (Ei(min), «#<?)• 

7. Draw the load line for Efimm) in the C 
divider; it passes through (.Eflmax), i c ) and 
[EArc(min)-b £i(max) — £i(min), i c — i NC 
( 1 +p/l —p)]. Produce to intersect E c 
(max) at i v . Then ip is the minimum tube 
current at e v =E 0 (max) for the most nearly 
unstable configuration, that is it is the cur¬ 
rent obtained from the assumed minimum 
plate characteristic. Hence this deter¬ 
mines the current scale for i NC and i c . 
Consequently all resistance values can now 
be obtained from this graph. 

Since this procedure has been based on 
lemma I, an analysis must show that any 
variation of one or more parameters 
within the allowed limits results in a 
stable flip-flop configuration. 



Figure 7. The influence of r 2 on switching time 
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No matter how the resistance values 
are calculated, it will in general happen 
that non-Radio Manufacturers Associa¬ 
tion values will result. Lemmas II, III, 
and IV may be applied to rectify this 
situation. Thus, since i v has been deter¬ 
mined using engineering judgment, it is 
quite reasonable to alter the preassumed 
value by some small percentage; this 
may be used to effect a small proportional 
rise or fall in all resistances. Again, a 
small percentage increase or decrease in 
E\ may be used to vary the resistances in 
a nonproportional manner. Finally, R s 
may be adjusted with negligible effect on 
At and R-> by varying E 2 . (Where these 
lemmas can not be used, as for example 
when Ei and E 2 cannot be altered, the 
non-Radio Manufacturers Association 
difficulty may be alleviated by designing 
with larger tolerances, p and <r, than engi¬ 
neering judgment demands.) 

The power of the graphical method is 
illustrated by two remarks. First, the 
broken lines on the right-hand side of 
Figure 5 instantly depict the influence of 
increased E x on the divider resistors. 
Since i P is given and does not change, 
larger E\ results in larger i NC and i e and 
hence in smaller R 2 and Rs ; at the same 
time, Rx increases as is evident from in¬ 
spection of the intersection of the two 
load lines for Ah (min) with the voltage 
axis. Second, the graphical approach is 
easily extended to take account of grid 
current at zero grid volts; the method 
is indicated by the insert top left of Figure 
5. The author has used the graphical 
method for a case when the analytic ex¬ 
pressions are particularly complex; this 
was the design of a flip-flop with a guaran¬ 
teed minimum grid current and also the 
ability to drive an NE-17, connected to 
the anode, with a prescribed minimum 
current. The saving in design time was 
deemed to be substantial. 

General Switching Considerations 

There are at least two switching times 
of importance which must be distin¬ 
guished. The first is the repetition rate 
that is, the rate at which the flip-flop can 
switch repeatedly from one state to the 
other and back. Repetition rate is re¬ 
lated to the settling time of the flip-flop, 
that is, the time required for the flip-flop 
to reach its new quiescent state. 

The second switching time of impor¬ 
tance is the switching speed, that is, they 
speed at which an initially quiescent flip- 
flop changes state. The switching speed 
is intimately related to the minimum du¬ 
ration of triggerpulseforsuccessful switch¬ 
ing. It will be shown that in general 


switching speed and repetition rate can¬ 
not be maximized simultaneously. 

In both cases, the time of arrival at the 
new state must be defined a priori, say, 
through specification of output voltage. 

Switching time is necessarily dependent 
on the kind of triggering used and on the 
amount of energy provided by the trigger 
pulse. For example, the ultimate in 
speed is achieved by a trigger which, by 
direct application to all appropriate 
points, instantaneously switches all volt¬ 
ages to their new values. 

In practice, however, the trigger is ap¬ 
plied to only one electrode for a time suf¬ 
ficient to insure change of state. It is 
proposed to consider the case when the 
conducting grid is instantaneously cut 
off. The problem is then reduced to de¬ 
termining the rate of rise of grid voltage 
of the initially nonconducting tube in re¬ 
sponse to a step of current, ii, injected into 
the C divider between Rx and A 2 , as in¬ 
dicated in Figure 6 . It is easily shown 
that 


e g (t) = Gnc — (Gwc + Gc) X 

[(l+5)e-“ ( -6e“^] ( 10 ) 


where 


6 = 6 (r 2 ) = 


a(l — 7 r 2 ) 
(7 —a) 


7 .«=r^— m+Rj)r k ±S] 

2RiTjT]c 

S*= KRi+RM-^RiXRiTjn) 

Tk = RkCic 
Ri = R\ +i ? 2 +Rz 

RiTjTjc = R1T2T5 + RtrsTi+Rmn 

(Ri+ Rj)ric = (Ax-)- Ri)Tz-b(Ri-]rRs)Ti^r 

(Rs-\-Ri)r 

2 

Also it can be shown that S 2 ^ 0 and there¬ 
fore 7 ^ ce. The equality signs hold if and 
only if ri = r 2 =T 3 . 

In what follows, it is assumed that r 2 
is the only explicit variable available to 
the designer and that all other terms are 
parameters adjustable within the limits 
set by stability requirements. An in¬ 
spection of equation 10 shows that the 
transient term consists of two exponential 
terms, of which the first requires longer 
to decay since a < 7 . Moreover, the 
variation of 7 and a with t 2 can be ob¬ 
tained immediately from the expressions 
of equation 10. Thus it can be shown 


Lemma VIII: 

(a) 7 >H; 7 -*-H as nr* but limit y>H 

(b) a<H; limit a = 0 as r 2 —► k 


where 


(Rx+Rz+RsXRxts+Rstx) 

(RiTs 4 “ A 3 T 1 ) 2 + RiRzTs 2 + R 2 .RsTi 

When ri=r 3 , then H= 1/n and the 


equality signs in lemma VIII are obeyed. 
Moreover, note that the largest value of a. 
is less than the smallest value of 7 , re¬ 
gardless of the value of r 2 , for any specific 
design. 

The coefficient b has been defined in 
equation 10 ; in particular, it can be seen 
that 

6(0) =—-— >0 (11) 

7 — a 

It follows that both coefficients are posi¬ 
tive for r 2 =0. It can be shown that: 

Lemma IX: 

As T 2 increases, both coefficients decrease 
through zero and become negative. In 
particular, 6 = 0 when r 2 equals the smaller 
of ri and t 3 ; ( 1 + 6 ) equals zero when r 2 
equals the larger of n and r 3 . 

Since a can never exceed 7 , minimum 
settling time can be achieved only by 
eliminating the term exponential { — at), 
that is, by making (1+6) = 0. Thus 

Lemma X: 

To achieve minimum settling time, r 2 should 
be chosen equal to the larger of n and r ;f . 
Settling time is then dependent only upon 7 . 

Two observations are worthy of note. 
First, this value of 7 does not differ too 
markedly from 1 /ti. It follows that the 
magnitude of r x constitutes a good first 
approximation to flip-flop speed. This 
circumstance was utilized in the section 
entitled “Graphical Design Procedure” 
in choosing optimum E c (max) on the 
basis of minimum Rx. Second, when 
r 2 equals the larger of Tx and t 3 , the equa¬ 
tion for 7 simplifies to 


i?i+i? 2 +A 3 
(i?i +Ra )r 3 + Rstx 


Ti = T 2 >T 3 


i?l + f ? 2 + I ?3 . 

-y =---; Tl<.r 2 —t 3 

RiTS + ( Rz + As )ti 


( 12 ) 


As pointed out, equation 12 was chosen 
to satisfy minimum settling time. It is 
common experience, however, that rep¬ 
etition rate can be improved by intro¬ 
ducing a small amount of overshoot at the 
expense of settling time. It is not dif¬ 
ficult to show the following: 


Lemma XI: 

Overshoot occurs if and only if r 2 is greater 
than both n and r 3 . 

It has been found experimentally for 
one flip-flop design that maximum rep¬ 
etition rate was obtained when r 2 ex¬ 
ceeded the larger of Tx and r 3 by about 
10 per cent. 

To obtain maximum switching speed 
requires that the grid voltage e g rise at its 
maximum rate. By differentiating e„ 
with respect to time t in equation 10 , it 
can be shown that: 
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Lemma XII: 


For maximum switching speed, r 2 should be 
so chosen that 


7 " 2 > > 


R2T3T1 
i?ir 3 —l-i? 3 Ti 


or C 2 > > 


C\Cj 

C\-\- Cs 


In words, C<> should be chosen much 
larger than the capacitance of the series 
connection of C\ and C 3 . 

The influence of r 2 on switching speed 
is illustrated by the curves of Figure 7, 
which show e g versus time for various 
values of t->. Note that these curves are 
pictorial verification of lemmas X, XI, 
and XII. 

The use of large r 2 to obtain maximum 
switching speed implies the existence of 
appreciable overshoot. In the case of the 
rising grid, this is of no consequence since 
the grid is clamped at cathode voltage. 
On the other hand, it is the same network 
which is involved in fall time considera¬ 
tions when the flip-flop returns to its 


original state. Consequently, there will 
be a similar overshoot in the case of the 
falling grid, probably of increased mag¬ 
nitude because of the initial grid clamp¬ 
ing. The overshoot reduces the ultimate 
repetition rate since it causes the grid, 
when retriggered soon after having been 
triggered, to start put from a voltage be¬ 
low — G c and hence to take longer to rise. 
Thus it appears qualitatively that the 
two switching times, repetition rate and 
switching speed, cannot be optimized 
simultaneously. 

Conclusion 

This paper has been concerned with the 
design of Eccles-Jordan flip-flops. After 
restricting the scope of the discussion to a 
particular design and to specific assump¬ 
tions, analytic design expressions were ob¬ 
tained and a number of lemmas de¬ 
veloped with regard to optimizing param¬ 


eters while guaranteeing stability within 
stated parameter tolerances. Finally, a 
general graphical design procedure was 
developed. 

The use of a transpose capacitor to im¬ 
prove switching times was discussed and 
a number of lemmas developed to serve 
as a guide not only in the choice of trans¬ 
pose condenser but also in choice of other 
parameters (particularly Ri and i? 3 ) with¬ 
in the stability limits previously obtained. 

The paper makes no claim to complete¬ 
ness but is meant only to serve as a col¬ 
lection of rules to guide the design of flip- 
flops. 
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No Discussion 



T HE continuing need for magnetic 
core materials which can be used ef¬ 
fectively at high frequencies has stimu¬ 
lated the search for materials in which 
eddy-current effects produced by rapid 
changes in magnetic flux are at a mini- 
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Figure 1. Coil of cold rolled 
48-per-cent nickel-iron alloy 
1 /4-mil thick and ceramic bob¬ 
bins 


in powdered metals or thin laminations of 
metal strip. 

The ferrites, which have resistivities 
about 10 8 to 10 12 higher than conven¬ 
tional magnetic alloys, have been used 
extensively in television sweep trans¬ 
formers where they may be operated at 
high frequency with low power loss . 1 
However, such characteristics as rel¬ 
atively low magnetic saturation, low 
Curie point and certain dielectric proper¬ 
ties, particularly for larger masses, re¬ 
strict the scope of application of the fer¬ 
rites. On the other hand, the fine sub¬ 
division of particles in powdered metal 


Ultrathin Tapes of Magnetic Alloys 
with Rectangular Hysteresis Loops 


MARTIN F. LITTMANN 

NONMEMBER AIEE 


Synopsis: The principal soft magnetic 
alloys have been rolled to thicknesses be¬ 
tween 1/8 and 1 mil. Magnetic properties 
of toroidal cores wound from ultrathin 
tapes were studied using d-c magnetization. 
Even though the coercive force tends to in¬ 
crease with reduced thickness, the coercivity 
of ultrathin metal tapes is usually lower 
than that of ferrites or powdered metals. 
Relatively rectangular d-c hysteresis loops 
were obtained in three' of four alloys in¬ 
vestigated. 


mum. Generally speaking there are two 
ways in which eddy-current effects can be 
reduced: increasing the electrical resis¬ 
tivity as in the ferrites, or restricting the 
paths in which eddy-currents can flow, as 
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4-79 PERMALLO 


H - OERSTEDS 


Figure 2 (above). D-c magnetization curves of 1 /4-mil 
ultrathin metal tapes compared with other materials used at 

high frequencies 


l6 TB r /B* = .945 


H ( _ a 5 OERSTEDS 


-4 -2 0 


B r /B m = .940 


1/2 MIL 


1/4 MIL 


2 4-4-2 

H-OERSTEDS 


Figure 3 (right). D-c hysteresis loops of grain-oriented 3- 
per-cent silicon-iron 


cores results in very low d-c permeability 
even for the best alloys. Where high 
permeability is desirable, very thin metal 
tapes are useful and these are the basis 
of the present paper. 

The commercial use of magnetic alloys 
rolled to thin gauges grew very rapidly 
during World War II. For power trans¬ 
formers, laminated or tape-wound cores 
of strip 4- to 6 -mils thick were employed, 
and for radar pulse transfonners cores of 
1 - to 3 -mil-thick materials were used. 
This large-scale use of thin materials was, 
to a large degree, the result of develop¬ 
ment of economical methods of rolling 1 - 
to 6 -mil strip, primarily silicon-iron, in 
quantities of thousands of pounds. Al¬ 
ready during this period there were some 
applications involving rates of magnetiza¬ 
tion so high that materials even thinner 


The alloys were rolled as single strands 
in a width of 2 inches employing work 
rolls of very small diameter. After roll¬ 
ing, the strip was slit into ribbons 1/8- 
inch or 1/4-inch wide for winding into 


toroidal cores. Because of the fragile 
nature of small toroidal cores, often con¬ 
sisting of only a few wraps of tape, it was 
advantageous to wind the cores on ce¬ 
ramic bobbins as shown in the foreground 


Table I. Basic Properties of Alloys Investigated 

79-Per-Cent Ni, 79-Per-Cent Ni, 
3-Per-Cent Si, 48-Per-Ceat Ni, 4-Per-Cent Mo, 5-Per-Cent Mo, 
Balance Fe Balance Fe Balance Fe Balance Fe 

Magnetic Sat uration, Gausses.20,000.16,000.8,700.8,000 + 

Curie Temperature, Degrees Centi- 

grade. 700 . 440 . 420 . 400 

Resistivity, Microhm-Centimeters. 47. 45. 55..... . DU 

Density, Grams per Cubic Centi- 

meter.. 7.65. 8.25. 8.72. 8.77_ 


than 1 mil (0.025 millimeter) were re¬ 
quired. Small quantities of 4-79 Permal¬ 
loy tapes 1/4-mil thick were rolled on a 
small mill at the Bell Telephone Labora¬ 
tories for such applications . 2 The pro¬ 
duction of magnetic alloys in tapes as thin 
as 1/8 mil has now been undertaken by 
commercial organizations. 

Since 1950, four different alloys have 
been rolled to gauges less than 1 mil and 
fabricated into toroidal cores. The de¬ 
scriptive term ultrathin has been applied 
to these materials. The alloys chosen for 
study are commonly used soft magnetic 
materials having high permeability. The 
group was restricted to those alloys hav¬ 
ing relatively high resistivity in order to 
keep eddy-current effects to a minimum. 
The values of magnetic saturation, Curie 
temperature, and volume resistivity, as 
shown in Table I, are characteristic of the 
alloys themselves and probably do not 
change with thickness in the range of 1 
mil to 1/8 mil. 



Figure 4. D-c hysteresis loops of grain-oriented 48-per-cent nickel-iron 
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Figure 5. D-c hysteresis loops of 4-79 Permalloy showing pronounced rectangularity 



Figure 6. Comparison of hysteresis loops for ultraihin magnetic alloys 1 /4-mil thick. Data 

plotted to same scale 


of Figure 1. The bobbins furnished sup¬ 
port for the tape during annealing and 
also provided a frame for the copper wind¬ 
ings used in testing. The tape was in¬ 
sulated with a magnesia coating prior to 
annealing. Annealing was conducted in 
very dry, pure hydrogen at temperatures 
from 875 to 1,150 degrees centigrade. 
The particular magnetic properties ob¬ 
tainable depend not only on the com¬ 
position of the alloy but on the processing 
steps of the metal supplier, the fo rmi ng 
and annealing procedure used by the 
fabricator, and the stressed and tempera¬ 


ture to which the core may be subjected 
after it is completed. These factors, in¬ 
dividually or collectively as the case may 
be, all contribute to the resistivity, crys¬ 
tal and domain orientation, and stresses 
present in the lattice structure of the com¬ 
pleted core and thus govern the charac¬ 
teristics of the core material when a mag¬ 
netomotive force is applied. Naturally 
it would be most difficult to present a de¬ 
tailed picture of the magnetic properties 
of ultrathin magnetic alloys at this early 
stage of development. Our purpose is to 
present a general picture of the type of 


magnetic characteristics we may expect 
under d-c magnetization as determined by 
the ballistic galvanometer. The measure¬ 
ments were made at room temperature 
with the flux parallel to the rolling direc¬ 
tion of the strip. 

The d-c permeability, particularly at 
high inductions, of ultrathin metal tapes 
is very much higher than for other ma¬ 
terials used for high frequencies as shown 
in Figure 2. It may be noted that the d-c 
magnetization curves for the 1/4-mil 
tapes are very similar to those for the 
same alloys in the more familiar thick¬ 
nesses of 1 to 14 mils. For comparison, the 
magnetization curve is shown for a speci¬ 
men of 1-mil, grain-oriented 48-per-cent 
nickel-iron having a very rectangular 
hysteresis loop. Rectangular hysteresis 
loops are of considerable interest for ap¬ 
plications such as magnetic amplifiers 
and memory cores and, as will be shown, 
ultrathin metal tapes may be made to 
have hysteresis loops which are quite 
rectangular. 

For example, the loops shown in Figure 
3 for 1/2-mil and 1/4-mil grain-oriented* 
3-per-cent silicon-iron show a high degree 
of rectangularity. Comparing the width 
of the two loops it may be noted that the 
coercive force is almost doubled by the re¬ 
duction in thickness from 1/2 to 1/4 mil. 
At first thought this might be considered 
to be the result of an increased volume ef¬ 
fect of surface contamination during final 
annealing. Of course, these very thin 
tapes are particularly sensitive to damage 
from faulty atmosphere during annealing. 
However, the increase in coercive force 
with reduced thickness has been ob¬ 
served even under very carefully con¬ 
trolled annealing conditions. The sug¬ 
gestion is made that as thickness is re¬ 
duced there are changes in domain size 
and shape which cause an increase in co - 
ercive force. A contributing factor may be 
the type of crystal orientation present, 
particularly for material such as the 3-per¬ 
cent silicon-iron, which has a high coef¬ 
ficient of anisotropy. 

As shown in Figure 4, the coercive force 
for grain-oriented f 48-per-cent nickel-iron 
also increases with reduced thickness but 
not as rapidly as for the silicon-iron. 
Compared to the 3-per-cent silicon-iron 
the coercive force for 48-per-cent nickel- 
iron is roughly 10 times lower for thick¬ 
nesses below 1 mil. 

4-79 Permalloy is generally regarded as 
a material having high permeability par¬ 
ticularly at low inductions. Hence, the 

^Substantially that of a single crystal with a cube 
edge in the rolling direction and a face diagonal in 
the transverse direction. 

fSubstantially that of a single crystal with a cube 
edge both in the rolling and transverse directions. 
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Figure 7. Variation of ratio of residual to maximum induction with magnetizing force fori /4-mil 
cores with rectangular hysteresis loops. Data for 1-mil grain-oriented 48-per-cent nickel-iron 

included for comparison 


high degree of rectangularity of the hys¬ 
teresis loop which it was possible to ob¬ 
tain in ultrathin gauges of 4-79 Permalloy 
was hardly expected, see Figure 5. 

The hysteresis loops for 1/4-mil ex¬ 
perimental specimens of the various al¬ 
loys including Supermalloy are compared 
with each other on the same scale in Fig¬ 
ure 6 which illustrates how the ma¬ 
terials compare in a qualitative way. The 
comparison is intended to be on a qualita¬ 
tive basis since it will be understood that 
the properties for each alloy can be varied 
considerably and the particular specimens, 
data for which are presented here, were 
chosen because their hysteresis loops were 
quite rectangular. Again referring to 
Figure 6, the various shapes of loops illus¬ 
trate that the terms “square loop” or rec¬ 
tangular hysteresis loop are, at best, only 


descriptive and that important differences 
in degree exist depending on the criteria 
used. Such criteria may be the slopes of 
the sides of the loops, or the slope of the 
top of the loop as saturationisapproached. 

One useful way of measuring rectangu¬ 
larity is to note how the ratio of residual 
to peak or maximum induction B r /B m 
varies with peak or maximum magnetiz¬ 
ing force H m . The curves for H m versus 
Br/B m shown in Figure 7 represent the 
same cores used for Figure 6 and Figure 
2. Lower ratios of B r /B m may be 
achieved by variationof processing and an¬ 
nealing treatments. Very low ratios may 
be obtained by the introduction of a small 
air gap into the magnetic path. Superm¬ 
alloy is the only material among those 
studied here in 1/4-mil thickness which 
has not been found thus far to display a 


maximum ratio of B r /B m greater than 
0.9. 

The space factor or packing factor pres¬ 
ent in the experimental cores used for the 
tests described above is as low as 50 pei 
cent for 1/4-mil strip. For small cores 
containing only a few turns of strip, space 
factor is not a limitation. Where space 
and weight are limiting factors, the 
thickness of magnesia insulation would 
probably be a handicap. However, there 
is considerable promise that improved 
techniques may permit much higher 
space factors without impairing the value 
of the insulation. 

It is too early to state in detail what the 
potential capabilities of these ultrathin 
materials are under pulse excitation or 
continuous high-frequency operation. 
The48-per-cent nickel-iron and 4-79 Perm¬ 
alloy tapes show considerable promise 
for high-speed memory cores in com¬ 
puters. For magnetic amplifiers and 
radar pulse transformers, the low coercive 
force and high magnetic saturation of 
ultrathin alloy tapes are of interest, al¬ 
though it may not be desirable necessarily 
to have a rectangular d-c hysteresis loop. 
Ultra-thin 4-79 Permalloy appears to be 
particularly attractive since it has rela¬ 
tively high volume resistivity and can be 
made to have either a rectangular or 
rounded d-c hysteresis loop together with 
low coercive force. 
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The Current Status of Dynamic 
Stability Theory 

FRANK E. BOTHWELL 
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A LTHOUGPI there exists today no 
theory of dynamic stability of suf¬ 
ficient generality to satisfy the demands 
of engineering, there are, nevertheless, 
several tools which may be employed ad¬ 
vantageously by the engineer in many 


stability investigations. A short de¬ 
scription of these tools illustrated by cer¬ 
tain practical examples will be presented. 

The concepts of absolute and relative 
stabilities of dynamical systems are mat¬ 
ters of much concern in nearly all branches 


of the physical sciences and in many 
branches of the social sciences. Although 
a great amount of research effort has been 
expended in this field, it can be said safely 
that only the simplest stability problems 
have been solved satisfactorily. This 
statement applies not only to nonlinear 
systems, but also to many systems which 
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can be represented in their ranges of 
operation by entirely linear equations, as, 
for example, the linear differential equa¬ 
tion with periodic coefficients. 

The fundamental problem is the defi¬ 
nition of stability. It is difficult to find 
one definition which will suit all purposes. 
Asymptotic stability, as formulated by 
Liapounoff, 1 is a satisfactory tool for the 
investigation of the action of a regulator 
or a feedback amplifier but is unsatisfac¬ 
tory in the treatment of conservative sys¬ 
tems, unless all periodic behavior is 
termed unstable. Liapounoff’s criterion 
of stability of so-called steady motion 
rectifies this difficulty for cases (such as 
linear conservative systems) in which the 
period of oscillation is independent of 
amplitude. However, such a criterion ap¬ 
plied to the motion of celestial bodies 
would produce the embarassing conclusion 
that the motion of each planet is one of 
unstable equilibrium. Thus it is neces¬ 
sary to introduce a third concept; namely 
orbital stability. 2 

The three concepts, asymptotic sta¬ 
bility, orbital stability, and stability in 
the sense of Liapounoff (hereafter to be 
referred to as temporal stability) may be 
termed absolute stabilities. They do not 
provide, in themselves, a solution to all 
problems of definition. A servomech¬ 
anism, although absolutely and asymp¬ 
totically stable, may exhibit such a per¬ 
sistent oscillatory response as to render 
its usefulness no greater than if it had been 
unstable. Another concept must be in¬ 
troduced, therefore, and one is concerned 
with degree of damping or relative sta¬ 
bility. 

Although too little damping is objec¬ 
tionable in many systems, in a few cases 
zero damping or even a slight negative 
damping is permissible. Consider the 
case of the guided missile, the flight path 
of which is often characterized by oscilla¬ 
tions occasionally increasing in amplitude. 
Inasmuch as the total flight time of the 
missile may include no more than three or 
four oscillations, the magnitude of the 


negative damping is the proper criterion 
of merit of the control system and the 
concept of absolute stability is not im¬ 
portant. 

Methods 

Three types of motion will be discussed 
in the order of increasing complexity: 

1. Stationary points, or points of equi¬ 
librium, such as operating points of regula¬ 
tory devices. 

2. Periodic motions such as occur in the 
operation of the electric oscillator. 

3. Any bounded or unbounded, periodic or 
nonperiodic motion. 

Obviously item 3 includes 2, and item 2 
includes 1. However, items 1 and 2 are 
more commonly encountered than is item 
3, and there is some advantage in a sepa¬ 
rate treatment. 

The particular motion, the stability of 
which is under investigation, whether it 
be an equilibrium point or a periodic or 
nonperiodic wave form, is generally de¬ 
noted as the undisturbed motion, whereas 
any other motion of the dynamical system 
is called a disturbed motion. Broadly 
speaking, a system is said to be asymptot¬ 
ically stable if all disturbed motions 
sooner or later degenerate into the undis¬ 
turbed motion itself. Application of the 
phrase “all disturbed motions” must be 
made with reservation however. Con¬ 
sider the case of the periodic motion of a 
wrist watch. Any light shock will change 
its period momentarily, but eventually the 
motion reverts to its former value, and 
consequently the timepiece is labeled 
stable. However, if dropped from a sec¬ 
ond story window, the watch quickly 
reaches a condition of instability. 

The situation is exemplified by the elec¬ 
tric oscillator which will not oscillate un¬ 
less first dropped on the floor, or perhaps 
by the more practical case of the amplifier 
which “sings” when overloaded and con¬ 
tinues in this state after removal of the 
excitation. These two examples refer to a 


type of behavior called hard excitation, in 
which the system becomes unstable only 
when subjected to disturbances of suffi¬ 
ciently large value. The general problem 
of dynamic stability under hard excitation 
is very difficult. Any one problem usually 
can be solved only by obtaining an ap¬ 
proximate solution for the motion of the 
system subject to all possible boundary 
conditions. 

Liapounoff presents two methods for 
the investigation of stability under soft 
excitation, that is, stability of dynamical 
systems subject only to infinitesimal dis¬ 
turbances. These are called the first and 
second methods. The former involves 
the investigation of the behavior of all 
solutions which originate in an infinites¬ 
imal region surrounding the undisturbed 
motion. The second method involves the 
investigation of the existence of a function 
which, in connection with the given 
dynamical equations, possesses a certain 
prescribed behavior. The difficulty in 
finding such restricted functions has re¬ 
sulted in very little use of the second 
method, and consequently stability in¬ 
vestigations are effected by means of the 
first method. 

Definition 

Consider a dynamical system described 
by the set of differential equations 

dXi 

—7 =fi(t,x h x 2 ,..,,x n ), i — (1) 

at 

in which the functions /* are, in general, 
nonlinear. Let x iQ (t) denote a particular 
undisturbed motion, and let xft) denote 
any disturbed motion. The undisturbed 
motion is temporally stable if, and only if, 
for a given positive number e, no matter 
how small, it is possible to find another 
number 5 such that for all disturbed mo¬ 
tions satisfying the initial inequality 

|.Vio(0)-^(0)|<5 (2) 

the subsequent motion satisfies the in¬ 
equality 



Figure 1. Voltage regulated 
n-stage d-c generator 


224 


Bothwell—Current Status of Dynamics Stability Theory 


July 1952 




[xj'o(i) &j(/)j<Ce (3) 

If, in addition 

lim !#«>(/)— Xi(t)\ =0 (4) 

t —>- co 

the undisturbed motion is said to be 
asymptotically stable. If the undis¬ 
turbed motion x i0 (0 is periodic, equation 
3 requires the approach of the disturbed 
motion to a motion of exactly the same 
period as the period of the undisturbed 
motion. This last requirement may be 
relaxed by means of the concept of orbital 
stability. However, this particular type 
of stability will not be considered in this 
paper. 

The First Approximation 

The tool which so far has been found 
most useful (and which falls in the cate¬ 
gory of the first method) in the investiga¬ 
tion of dynamic stability, is the so-called 
method of first approximation. The pro¬ 
cedure is contained in the following six 
steps. 

Representation of the Physical 

System by Means of a Set of 
Dynamical Equations 

This first step is often the most dif¬ 
ficult. There is no general method of ap¬ 
proach, although the resultant dynamical 
representation is invariably in the form 
of differential, difference, or integral 
equations, or, in some cases, a mixed sys¬ 
tem involving two or three types. 

Search; for the Undisturbed 
Motions 

If the undisturbed motions are station¬ 
ary points, they are easily found by solv¬ 
ing the static equations obtained from the 
dynamic equations by deleting all de¬ 
rivatives and differences. However, if the 
undisturbed motions are periodically 
varying solutions of the dynamical rep¬ 
resentation, their determination is often 
most difficult. An approximate means 
of solution will be provided in the section 
on the method of equivalent lineariza¬ 
tion. When the undisturbed motion is 
nonperiodic and is prescribed by means of 
certain boundary or other characterizing 
conditions, there is no general method of 
approach for nonlinear systems. 

Perturbations 

Once the undisturbed motion has been 
found, all disturbed motions in the neigh¬ 
borhood of the undisturbed motion may 
be found by means of perturbations. The 
usual method involves the expansion of 
the original set of dynamical equations in 



Taylor series about the undisturbed mo¬ 
tion itself. Stability of the set of linear 
perturbation equations resulting from 
neglect of the higher order Taylor series 
terms generally, but not always, guaran¬ 
tees stability of the original dynamical 
equations. The cases in which the rule 
breaks down are usually those in which 
certain linear perturbation terms disap¬ 
pear. Inasmuch as one is usually con¬ 
cerned with the effect on absolute and rel¬ 
ative stabilities of system parameter 
variations, and since the rule breakdown 
generally occurs only when the system 
parameters are such that the system it¬ 
self is on the very edge of instability, such 
an occurrence requires no further atten¬ 
tion. However, there are cases in which 
the breakdown of the rule occurs for all 
parameter values, the method of first 
approximation fails, and one is forced to 
consider higher order Taylor series terms. 

Characteristic Equation 

If the original dynamical equations are 
autonomous, that is, if the dependence on 
time is not explicit, and if the undisturbed 
motion is a stationary point, the equa¬ 
tions of first approximation are linear and 
possess constant coefficients. Stability is 
consequently determined by the roots of 
the characteristic equation of the first ap¬ 
proximation. If the original system is 
described entirely by differential equations 
or entirely by difference equations, the 
characteristic equation is either an alge¬ 
braic equation or can be made so by a 
logarithmic transformation. However, if 
the system is represented by a mixed sys¬ 
tem of differential-difference equations, 
the characteristic equation is transcen¬ 


dental and presents certain difficulties of 
solution. 

Stability Criteria 

The occurrence of one or more roots of 
the characteristic equation with positive 
real part is the necessary and sufficient 
condition for instability of the original 
system. Whether or not this situation 
occurs can be determined by means of 
either the numerical method of Nyquist 3 
or the analytical method of Routh 4 and 
Hurwitz. 5 The choice of method depends 
on the complexity of the characteristic 
equation and on whether one is concerned 
with parameter variations. 6 Both the 
Routh method and the Nyquist method 
were originally derived for the treatment 
of absolute stability. However each may 
be modified to yield certain information 
concerning the actual location of the 
roots. 6 " 9 

Stability Boundaries 

The final step in any stability analysis is 
the presentation of the absolute and rel¬ 
ative stability boundaries as functions 
of basic system parameters. A few ex¬ 
amples will now be given. 

Voltage Regulation 10 

Let us follow through the six steps for 
the n-stage d-c generator, shown in Figure 
1, which employs a 3-terminal nonlinear 
regulator C, such as the step resistor, the 
resistance value of which varies in steps 
with the control voltage V . 

1 Denoting the voltage function of the 
regulator Cby/Oi, V ), the dynamical equa¬ 
tions are 
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E = Ei+f 


(A*) 


(5) 


dEs/dl = (Ri/LiXKiEi-Eo) 
dEi/dt = {R»JL,)(KiE«-E 5 ) 

dEn/dt = (Rn-l/Ln-i)(Kn-iE n __i — E n ) 
dV/dt — ( Rn / 'Ln){K<>K n E n — F) 

2. The undisturbed motion is a stationary- 
point, the operating point of the system. 
It is obtained by deleting the derivative 
terms in equation 5 and solving the resulting 
equations. Denoting the equilibrium volt¬ 
ages by subscript zero, it may readily be 
shown that 


E‘ia — KiEio Em — K\KiEiq 

E n o = KiK-i.. . K n _iEiQ 


( 6 ) 


V^K.KJU. .. K n E w 


in which Ew is given by the nonlinear equa¬ 
tion 


E = Eiq ~\-f(Eu,/R uKqKi .. ,K n Eia) 


C7) 


3. Setting Ei = E ia -\-ei and V—Vq+v, one 
finds the set of linear perturbation equations 

ei +(fn/KiRi)ez+f v v=0 

de<i/dt +(i?i/ L{)e o —KfRi / Li)e i = 0 

( 8 ) 

de,i/dt-\-(R n — 1 / Ln—l)^n b\-n _I X 

( Rn—l/En—l)^n—\ — 0 

dv/dt+(R n /L n )v-K a K n (R n /L n )e n - 0 

in which //, and f v denote the partial de¬ 
rivatives of/with respect to h and V, evalu¬ 
ated at the equilibrium point. 

4. The characteristic equation 8 is ob¬ 
tained by equating to zero the determinant 
of the coefficients of the unknown voltages 

and v, after replacing the operation of 
differentiation by the symbol p. It is 

FS )1 


Ri 

L t+ ZV 


r 

> -0 (9) 

-^1-^2. . .-Ln 

in which K=K,K U . .K n 

5. _ The Hurwitz determinant is useful in 
this case. For a 3-stage generator, equation 
9 is the cubic 


aop 3 + ai p^ -f aap +a 3 = 0 

in. which 

<Zq=1 


( 10 ) 


* \£i W 


a 2 




R S R, 



' L,L,l,\ 1+f x+ K fr) 


( 11 ) 


The stability boundary is obtained by 
il c teri n in ant Zer ° ^ penultimate Hurwitz 


l a i ao 
l a s a o 
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Using equation 11, one finds an expression 
for the stability boundary 

(l/W-F a 3 ')( 1+—■;+—t') —1-f-a/ (12) 

\ «2 aa / 

shown in Figure 2, in which 


9 ( l-f-ai' / —l)-j-a2" 


X 


l+|(\/iw-l) x 


ai'**K- 


fv 


(I4) 


1+ 


t R% Li 

jiu «2 = — 


Ri 


Ri R\ ~ffii 

Rz Li 


in which 


«3'=— 


Lz Ri -+-//, 


(13) 


„ , E-c E 

*-*&=»**'■ "**7. 


The region of stability lies below the surface 
considered as a contour map. 

6. Except for the case of a linear control 
element C, the partial derivatives//, and fv 
are functions of the location of the operating 
point, which itself is determined by the 
circuit parameter values. Therefore equa¬ 
tion 12 and Figure 2 do not express, in 
general, the stability boundary explicitly in 
terms of basic system parameters, and, in 
fact, the problem cannot be completely 
solved until the form of the control function 
f(R, V) is known. A common control ele¬ 
ment for high-voltage generators is the step 
resistor. The essential features of its sta¬ 
bility in connection with Figure 1 may be 
observed by assuming that the resistance 
varies linearly and continuously with con¬ 
trol voltage when the latter exceeds a critical 
value, V c - Thus 


a.% 


n _ 


_ R’tLi 

K c -1 URt 


az r, = 


1 RzLx 
K c ~ 1 L 3 Ri 


Rc = 0 


if V< V e 


Ri 


=E C —(V-V C ) if V>V C 
so that 


(15) 

In Figure 3, the region of instability lies in¬ 
side the mound enclosed by the coordinate 
plane a/' — G and the contour surface. 

Shunt Feedback Amplifier 6 

An excellent example of the concept of 
relative stability is afforded by the treat¬ 
ment of the circuit of Figure 4. For the 
case in which the admittance Y is com¬ 
posed of a parallel combination of resist¬ 
ance and capacitance R and C, one sys¬ 
tem time constant is always equal to RC. 
Denoting by RC/b, where — oo<&< 
1, the value of the maximum time con¬ 
stant of the system, the region of relative 
stability is given by the inequality 

1 - 2 

tt3<- for a 2 <- Oil 

Oil — 012 3 


f(h,V) = 0 


if V<V C 


Ri 


< 


=E c -f(V-V c )l\ if V> V c 

y c 

and the stability boundary takes the form of 
Figure 3, or 


2(2aicx-z) 2 
2oci-\-oi‘z 


for cti 


(16) 


in which 


^ KiRSz R*SiS 2 
«i = Ao, az — —-—, a 3 = - 


1-5 


(1-5)= 



If equality signs are inserted in equation 
16, one obtains the locus of parameter 
values for which the maximum time con¬ 
stant in RC/b. The results are shown in 
Figure 5. 

Periodic Solutions 

So far in the examples only the undis¬ 
turbed motions which were stationary 
points have been considered. The opera¬ 
tion of a nonlinear oscillator is often de¬ 
scribed by Van der Pol’s equation 
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d 2 x , dx 


(17) 



For fi small, one undisturbed motion is 
given approximately by 

.v 0 (/) = 2 sin / (18) 

The equations of first approximation are 
obtained by the substitution of 

x(/)=xo(J)+5(0 

together with equation 18, into equation 
17 and omitting all but linear terms in 5. 
The result is the linear differential equa¬ 
tion with periodic coefficients 

d*8 , . , d5 

~+m( 4 sm*1-1)-+ 

(1+8 sin t cos t)8 — 0 (19) 

Here the methods of Nyquist and Routh 
are of no use, and the general method, as 
outlined by Stoker, 11 for example, stops 
far short of the actual solution for sta¬ 
bility boundaries. In fact, the bound¬ 
aries can be obtained only after one has 
evaluated a certain infinite order deter¬ 
minant called a Hill determinant. 

For the very simple case of the Mathieu 
equation 


-—b(5 + « cos t)x — Q (20) 

at 2 

which represents a variety of physical 
systems, 12 the stability boundary has been 
completely determined and is shown in 
Figure 6. Even in this simple case the 
boundary is a very complicated function 
of the parameters 8 and e. 

Equivalent Linearization 

A very powerful tool for the investiga¬ 
tion of stability of periodic motion is 
based on the so-called method of equiva¬ 
lent linearization as formulated by Kryloff 
and Bogoliuboff, 13 and apparently first 
linked with the Nyquist criterion by 
Kochenburger 14 in his treatment of the re¬ 
lay servomechanism. The method in¬ 
volves the systematic replacement of all 
nonlinear elements in the system by 
quasi-linear elements which have the 
property of preservation of waveform but 
nonlinear magnification. For example, in 
the Van der Pol equation 17, the non¬ 



linear term 


(* 2 -l)^ (21) 

dt 

is replaced by a linear term which is the 
result of considering only the fundamental 
term in the response of equation 21 to the 
excitation x~a sin cat. But the funda¬ 
mental term is (a 2 /4— l)acocos cat, so that 
the equivalent linear term is (a 2 /4—1) 
(, dx/di ), and equation 17 is replaced by 


d 2 x ( a 2 



( 22 ) 


Thus oscillations which have amplitude 
greater than a —2 are obviously damped, 
while oscillations of magnitude less than 
a —2 obviously are characterized by 
negative damping. Therefore, the solu¬ 
tion x—2 sin t is stable. 

Inasmuch as the disturbed motion (and 
for that fact even the undisturbed mo¬ 
tion) is not exactly sinusoidal, there is 
some doubt as to what is meant by the 
amplitude of motion. In fact, the whole 
procedure rests on rather flimsy mathe¬ 
matical ground. However, in systems 
which are essentially low-pass filters, such 
as servomechanisms, wherein harmonic 
distortion is low, the procedure has been 
very successful, and this is its great re¬ 
deeming feature. 

A remarkable example of the use of 
equivalent linearization occurs in the 
field of economics. A model of the busi¬ 
ness cycle has been hypothesized by 
Goodwin 16 in the form of the autonomous 
nonlinear mixed differential-difference 
equation 

73>(2+0)+(l — oc ) y ( t + d ) = (23) 

in which y(t) denotes present national in¬ 
come, and the induced investment func¬ 
tion <j>(x) is the unsymmetrical saturation 
curve of Figure 7. If c denotes the am¬ 
plitude of the quasi-linear oscillation of 
y(/), the characteristic equation of the 
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Figure 7 (left), 
induced invest¬ 
ment function 



corresponding linearized equation is found 
to be of the form 

h(z)=g(c) (24) 

The Nyquist plot of each side of equation 
24 is shown in Figure 8. Here c is plotted 
as a function of g(c) so that the intersec¬ 
tion of the spiral h(iu) (co=2 ir/p, where p 
is the period of the assumed oscillation) 
with the g{c) axis determines an infinite 
number of admissible values of period p. 
That is, it is possible that the system os¬ 
cillates in any one of an infinite number of 


Figure 8 (right). 
Nyquist plot of 
quasi-linear rep¬ 
resentation of 
Goodwin’s mod¬ 
el of the busi¬ 
ness cycle 



mess 



/i=(l "a) 6/y 


X-"V m - '/3 


r= Vy * 4.0 


modes. The admissible values of p are Undisturbed motions which are sta- volume li, January 1932, puges 126. 47. 


approximately 8 years, 1 year, 1/2 year, 
1/3 year, and so forth. 

The question of stability is decided by 
the following reasoning. For values of c 
for which the spiral encloses the point 
g(c), the motion is stable, c will decrease, 
md g(c) will increase. For values of c for 
vhich the spiral does not enclose the 
point g(c), the motion is unstable, c will 
increase, and g{c) will decrease. Obvi¬ 
ously each of the infinite number of modes 
is stable. It is also apparent that a rela¬ 
tively small disturbance of a high order 
mode is sufficient to cause oscillation in a 
different mode, a result which is not 
available in the analysis of infinitesimal 
disturbances. That these modes do exist 
and are stable, has been demonstrated 
experimentally by means of analogue 
equipment. 

Conclusions 

Of the three types of stability, orbital, 
temporal, and asymptotic, the last has re¬ 
ceived the most attention. Almost ail the 
successful stability analyses have em- 
ployed the first of Liapounoff’s two 
methods, and have been concerned mainly 
with systems subject only to infinitesimal 
disturbances. 


tionary points almost always are treated 
successfully by a sequence of first order 
perturbation and a stability criterion of 
Nyquist or Routh type. Relative sta¬ 
bility as well as absolute stability may be 
treated by this method. 

Undisturbed motions which are periodic 
are currently treated by two methods. 
The perturbation method invariably 
produces a linear equation with periodic 
coefficients, the stability of which may be 
investigated by Hill’s method. However 
the general method stops far short of a 
solution, and only the simplest cases have 
been solved. Equivalent linearization has 
been currently successful in the analysis 
of low-pass systems. However, it is a 
very approximate method, refinements 
are involved, it fails completely on high- 
pass systems, and mathematical justifica¬ 
tion is lacking. The general problem of 
the stability of periodic motion remains 
unsolved. 
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Types of Magnetic Amplifiers 

Survey 


JAMES G. MILES 
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Synopsis: Many forms of magnetic ampli¬ 
fiers have been described in the past 50 years. 
The patents and writings relating to these 
devices have been examined to determine 
what types of magnetic amplifiers exist and 
what properties these devices have in com¬ 
mon. The objects of this paper are: (1) 
to describe the various types of magnetic 
amplifiers, and (2) to present a method by 
which the basic principles relating to these 
devices may be codified according to an 
orderly system. 

Magnetic Amplifier Chart 

A chart which is basic to a scheme for 
classifying the various types of mag¬ 
netic amplifiers is given below. Any 
known magnetic amplifier may be repre¬ 
sented by a combination of the key letters 
and numbers listed in the chart. Each of 
the letters and numbers represents infor¬ 
mation relating to the magnetic amplifier 
class, circuit type, and special means for 
affecting the characteristics of the par¬ 
ticular amplifier to be described. 


In defining classes A and B of the chart 
an arbitrary distinction has been made 
between inductive-impedance devices (A) 
in which the a-c power is applied to the 
core via the load or reactance winding, 
and (B) in which the a-c power is applied 
to the core via a separate primary wind¬ 
ing. According to this distinction, there¬ 
fore, key A refers to saturable reactors 
controlled by signal magnetomotive force 
and key B refers to saturable transformers 
controlled by signal magnetomotive force. 

Two classes of signal-magnetomotive- 
force-controlled saturable transformers 
appear. The first, key B', relates to a 
class of saturable transformers in which 
control is effected by varying the incre¬ 
mental permeability or by otherwise dis¬ 
torting the flux-time wave form in the 
core. The second class, key B", relates to 
the class of transformers in which control 
is effected by varying the reluctance of a 
magnetic path which shunts the trans¬ 
former leg linking the secondary winding 


and thus controlling the mutual coupling 
between the primary and secondary wind¬ 
ings. 

In defining the specific circuit types in¬ 
dicated on the chart, push-pull (symmet¬ 
rical) circuits, key 1, are to be distin¬ 
guished from single-sided (unsvmmet- 
rical) circuits, key 2. 

It is to be understood that “direct cur¬ 
rent” is a generic term which represents 
the signal, bias, or control effects that 
actually may be transitory or otherwise 
variable with respect to time. A general 
characteristic of magnetic amplifiers is 
that the power-supply frequency is con¬ 
siderably higher than any component of 
input-signal frequency. 

To simplify the classification, only 
single-stage magnetic amplifiers are con¬ 
sidered. Cascaded arrangements of stages 
are not included, as such, in this classi¬ 
fication, but considerations are stated 
which must be taken into account in the 
design of single magnetic-amplifier stages 
which are to be used in cascade with other 
stages. 


Paper 52-216, recommended by the AIEE Magnetic 
Amplifiers Committee and approved by the AIEE 
Technical Program Committee for presentation at 
the AIEE Summer General Meeting, Minneapolis, 
Minn., June 23-27, 1952. Manuscript submitted 
December 4, 1951; made available for printing 
May 14, 1952. 

James G. Miles is with Engineering Research 
Associates, Inc., St. Paul, Minn. 


Classification Chart for Magnetic Amplifiers 

General Classes 

A. Saturable reactors controlled by signal magnetomotive force 

B. Saturable transformers controlled by signal magnetomotive force 

B'. Control by varying incremental permeability or by distorting wave form 
Bh Control of mutual inductance by varying the reluctance of a saturable shunt 

C. Magnetic modulators with negative resistance characteristics 

General classes A and B contain specific circuit types such as polarized and nonpolarized magnetic amplifiers comprising: 

1. Single-sided circuits 

2. Push-pull circuits 

3. Wheatstone-bridge-type circuits 

4. Special magnetic amplifier circuits (such as cross-connected pairs of magnetic amplifiers, magnetic amplifiers in com¬ 
bination with transformers, and so forth) 

Amplifier sensitivity, stability, linearity, and speed of response may be affected by incorporating any one or more of the follow¬ 
ing means: 

a. External feedback (amplitude-related and/or phase-related external feedback means including impedance coupled 
means terminating in a single-control winding or special feedback winding) 

b. Internal feedback (by rectifiers seriesed with load windings) 

c. Unidirectional bias magnetomotive force of fixed magnitude 

d. A-c bias magnetomotive force at fundamental or harmonic of power supply (carrier) frequency 

e. Off-resonant and ferroresonant means 

f. Frequency-discriminating means 

g. Additional co-operative wave-form distorting means 

h. Ultra-carrier frequency magnetomotive force means for reduction of hysteresis effects 
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Fisure 3. Variable mutual-inductance type d-c saturable transformer 


Figure 1. Single-sided d-c saturable reactor without feedback 


Examples of Types of Magnetic 
Amplifiers 

Figures 1 through 20 show typical mag¬ 
netic amplifiers which may be classified 
according to the chart. At least one ex¬ 
ample of a circuit relating to each key is 
shown in the figures. The examples have 
been chosen not on the basis of dates of 
invention, but because they illustrate par¬ 
ticularly well the circuit types described. 

Figure 1 (class A, 1, c) is a schematic of 
the first known d-c controlled saturable 
reactor for which gain was claimed. E. F. 
W. Alexanderson 1 described this device in 
United States Patent 1,206,643 (1916), 
and as a basic device it has been used 
widely. Alexanderson’s device comprises 
two variable-reactance windings as load- 
current control means, the reactance 


windings being placed in series between 
the a-c load and the a-c power supply. 
The total impedance of the circuit is con¬ 
trolled according to the reactance of the 
reactance windings. Control is accom¬ 
plished in a well-known manner based on 
the change of incremental permeability 
of the core according to the amount of d-c 
saturating magnetomotive force applied 
to the core. The saturating windings are 
connected so that opposing a-c electromo¬ 
tive forces are induced in them by the re¬ 
actance windings; the result is that zero 
alternating voltage of fundamental power- 
supply frequency appears across the ter¬ 
minals of the d-c control circuit. Alex¬ 
anderson provided separate d-c bias wind¬ 
ings for the purpose of introducing a d-c 
bias magnetomotive force to locate an op¬ 
timum operating point on the magnetic 
characteristics. 


Figure 2 is an example of a d-c con¬ 
trolled saturable transformer of the vari¬ 
able incremental-permeability type (class 
B', 1). This circuit is due to A. E. Bowen, 2 
United States Patent 2,230,558 (1941). 
Here the primary and secondary windings 
are wound on the outside legs of a 3-leg- 
ged core in such directions as to be mu¬ 
tually coupled and so as to pass no flux 
of power-supply fundamental frequency 
through the center leg of the core. Con¬ 
trol is effected by imposing, via the satu¬ 
rating winding, a condition of d-c pre¬ 
magnetization. If a core is premagnetized 
by a d-c component of magnetomotive 
force, a given a-c magnetomotive force in¬ 
troduced by the primary winding will ef¬ 
fect less change in flux than if the core 
were not premagnetized. Secondary 
voltage therefore may be controlled in a 
core by introducing a controlled amount 


"Back-Contact" "Front-Contact" 



Figure 2. Variable-permeability type d-c saturable transformer Figure 4. Variable mutual-inductance type d-c saturable transformer 
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External Feedback 
Signal and d-c Bias 




Figure 5 (right). Single-sided d-c saturable reactor with combined 

external feedback 

Figure 7 (above). Flux-normalizing method of control for internal- 
feedback magnetic amplifier 


of d-c premagnetization so as to effect a 
control of the a-c flux variation in the 
magnetic circuit linking the secondary 
winding with the primary winding. The 
amplitude of the a-c flux which mutually 
couples the primary and secondary is con¬ 
trolled by the degree of d-c premagnetiza¬ 
tion. 

A d-c saturable-shunt type of variable- 
mutual-inductance transformer is shown 
in Figure 3. This circuit, due to E. R. 
Stoekle, 3 United States Patent 1,376,978 
(1921), (class 3", 1 ), shows a 4-legged 
transformer in which the flux that couples 
the secondary with the primary may be 
shunted, to a controlled degree, by a vari¬ 
able reluctance shunt. The mutual flux 
that links the secondary with the primary 
winding thus is decreased and controlled 
effectively. Control is effected in this 
typical circuit by means of saturating 
windings arranged to cause “circulating” 
d-c flux in the two saturable shunt legs. 
In this manner, the reluctance of the 
shunt legs to the passage of a-c flux is 
controlled effectively. In Figure 4, the 
transformer primary winding is inserted 
as a controlled reactance in series between 


the load and a-c power source. The re¬ 
actance of the primary winding in this cir¬ 
cuit is controlled according to the com¬ 
bined reflected impedances from the two 
secondary windings with which it is imme¬ 
diately coupled, the reflected impedances 
being controlled according to the relative 
mutual coupling between the primary and 
secondary windings. 

Another example of the d-c saturable 
shunt-type variable-mutual-inductance 
transformer is that shown in Figure 4. 
C. M. Hines 4 in United States Patent 
2,215,820 (1940) described a saturable 
transformer (class B u , 1 ) similar to that 
shown in Figure 3 but having another set 
of secondary windings placed on the d-c 
saturable legs as shown. These “back 
contact” secondary windings are arranged 
to be coupled effectively with the primary 
windings when the control legs are de- 
saturated. On the other hand, if the d-c 
control legs are saturated, a minimum of 
primary flux will link the back contact 
secondary windings and a maximum of 
primary flux will link the “front contact” 
secondary winding. This circuit, with the 
primary winding excited by the a-c power 


supply and with the secondary windings 
exciting separate loads, is, in effect, a 
single-pole double-throw a-c power 
switch. 

Means are abundant in the magnetic- 
amplifier art by which either positive or 
negative feedback may be introduced to 
affect gain, stability, linearity, and speed 
of response of magnetic amplifiers. Feed¬ 
back means may be divided into two 
principal classes: external feedback and 
internal feedback. These classes are de¬ 
scribed in the following paragraphs. 

In external feedback circuits, a measure 
of output voltage and/or current is re¬ 
introduced as a feedback magnetomotive 
force via an external circuit which may 
terminate either in the control winding or 
in a special feedback winding. The 
amount of external feedback may be 
modified by adjusting either the ampli¬ 
tude or the phase of the feed-back signal 
with respect to the input or control signal. 

One of the earliest examples of the ap¬ 
plication of amplitude-related external 
feedback to a d-c saturable-reactor mag¬ 
netic amplifier is shown in Figure 5 (class 
A,l,a,cJ). Here E. F. W. Alexanderson 5 


Figure 6 (below). Impedance-coupled external feedback means 


Figure 8 (right). Internal-feedback magnetic amplifier with saturating 
winding control for a-c load 




Impedance 
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D-C Control 



Figure 9. internal- 
feedback magnetic 
amplifier with satu¬ 
rating winding con¬ 
trol for d-c load 


in United States Patent 1,328,797 (1920) 
showed a method for obtaining load-cur¬ 
rent control by placing a controlled re¬ 
actance winding in parallel with the load 
and across the secondary terminals of an 
a-c supply transformer. A measure of 
magnetic-amplifier output potential (also 
load potential) is applied via a second 
transformer and rectifier circuit to an ex¬ 
ternal feedback winding in such manner 
as to aid the effect of the control magneto¬ 
motive force. D-c bias is introduced via 
a parallel circuit using the same winding 
as used for external feedback, as shown. 

An alternative means for introducing 
amplitude-related external feedback is 
shown in Figure 6 which illustrates, in 
block diagram form, an impedance- 
coupled external feedback means (class 
-, l,a) described by A. S. FitzGerald 6 in 
United States Patent 2,464,639 (1949). A 
measure of load current creates a potential 
drop across the impedance Z. The im¬ 
pedance is so connected that the potential 
drop across it will aid (or oppose) the con¬ 


trol signal applied at the magnetic-am¬ 
plifier input. This action amounts to 
feedback. 

A phase-related external feedback cir¬ 
cuit is described by J. Slepian 7 in United 
States Patent 1,645,302 (1927), not illus¬ 
trated here. In this patent a series-res¬ 
onant inductance - capacitance phase 
shifter in the control circuit fixed the 
phase of the control signal with respect 
to the feed-back signal. 

The so-called “internal feedback’’ 
means is the second important means for 
effecting feedback. In order to achieve 
internal feedback, rectifiers are inserted in 
series between the load windings and the 
load circuit. Since only unidirectional 
current impulses can flow in the load 
windings of such circuits, unidirectional 
pulsating flux components are created in 
the reactor cores in magnitude according 
to the magnitude of the load current. 
Control of such devices is effected in one 
of several ways to be discussed in the fol¬ 
lowing paragraphs. 


One of the earliest described means for 
effecting control in an internal feedback 
circuit is shown in Figure 7 (class A,l,b,d). 
This method is due to F. G. Logan, 8 
United States Patent 1,997,179 (1935). A 
“normalizing shunt” is placed in shunt 
with the rectifier circuit so that a con¬ 
trolled amount of reverse magnetomotive, 
force may be applied to the reactor core 
during its quiescent half-cycle. By means 
of this shunt circuit, a current is caused to 
flow at such times as the reactor otherwise 
would be inactive, the direction of this 
current being opposite to that in which 
the current will flow when the reactor 
again becomes active in carrying load 
current. The flux and energy conditions 
thus tend to establish, to a degree deter¬ 
mined by the adjustment of the normaliz¬ 
ing impedance, the conditions which 
would have existed if the reactor had not 
been subjected to an inactive condition. 
According to the value of the normalizing 
impedance, therefore, a controlled amount 
of energy and current is delivered by the 
load circuit. The changes in nornmlizing- 
impedance current required to effect con¬ 
trol are claimed to be very slight com¬ 
pared to the changes in load current 
which may be obtained. 

Another and more widely used means of 
control for internal-feedback magnetic 
amplifiers is shown in Figures S and 9. 
This means of magnetic-amplifier control 
employs a d-c saturating winding to pre¬ 
saturate the core by an amount adjusted 
by the control current. The control wind¬ 
ing links the magnetic core arrangement 
in such manner that the d-c control mag- 



Figure 10. Wheatstone-bridge type magnetic amplifier 


Figure 11. Wheatstone-bridge type magnetic amplifier with internal 

and external feedback 
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Figure 12. Push- 
pull magnetic 
amplifier with in¬ 
ternal feedback 
and with fre¬ 
quency-discrimi¬ 
nating means 


Flux vs. MMF 



Characteristics 


Sec. E vs. Time 



Three -Element Polarized Telegraph 
Signals to be Amplified 


Note: This is twice the power of the 
output of a single secondary and is 
Double Frequency with respect to 
a-c power supply frequency. 


Figure 13. “Sec¬ 
ond - harmonic 
type" push-pull 
magnetic ampli¬ 
fier with internal 
feedback and 
with frequency - 
discriminating 
means 


netomotive force either will be aided or 
opposed by the unidirectional internal- 
feedback magnetomotive force set up by 
the load windings. Figure 8 shows an in¬ 
ternal-feedback magnetic amplifier (class 
A,l,b,d) for the control of an a-c load. 
This circuit is due to F. G. Logan, 9 United 


States Patent 2,126,790 (1938). Figure 9 
illustrates a similar d-c signal-winding- 
controlled internal-feedback magnetic am¬ 
plifier (class A,l,b,d ) for the control of a 
d-c load circuit. This circuit was de¬ 
scribed by S. M. Hanely, 10 United States 
Patent 2,144,290 (1939). 


Magnetic amplifiers may be arranged in 
the form of a Wheatstone bridge. Such 
circuits are advantageous in that they 
may be set up readily to have polarized 
characteristics, that is, they may be made 
to deliver positive output corresponding 
to a positive control signal, zero output 
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D-C Control 



Load Windings 


Note: When the core is under the influence of 
control mmf, A~C Power Supply is suf¬ 
ficient to saturate appreciably the core 
once each cycle. 


Symmetrical non-linear impedance elements 
without unidirectional properties per se. 
Examples: Thyrite or Parallel Back-to- 
Front Non-Symmetrical Recti¬ 
fiers. 


Characteristic: 


Voltage 



Time 


Figure 14. Un¬ 
balanced har¬ 
monic device 
employing 

1. Symmetrical 
current-distorting 

impedances 

2. Asymmetrical 
current-distorting 
impedances with 
polarizing means 
(d-c saturable re¬ 
actors) 



Figure 15. Regu¬ 
lating apparatus 
with 

1. Inductance- 
capacitance off- 
resonant vari¬ 
able - impedance 

circuit 

2. Resistance-in¬ 
ductance - capa¬ 
citance ferroreso- 
nant control 

circuit 


corresponding to zero control signal, or 
negative output current corresponding to 
a negative control signal. These bridges 
have the advantage that the control, out¬ 
put, and power-supply circuits are prac¬ 
tically isolated one from the other. 

Many Wheatstone-bridge circuits have 
been described in which one or more arms 
of the bridge comprise the reactance 
windings of magnetic amplifiers. A 
specific example is the class A,3 circuit 
shown in Figure 10, this circuit being due 
to L. W. Thompson, 11 United States 


Patent 1,844,704 (1932). All four arms 
of the bridge are the reactance windings 
of d-c magnetomotive-force-controlled 
saturable reactors. Internal and external 
feedback means are shown in Figure 11 
as applied to the Wheatstone-bridge type 
magnetic amplifier due to Uno Lamm, 12 
United States Patent 2,403,891 (1946) 
(class A,3,a,b). 

A class B',2,b,f,g device is disclosed by 
E. M. Boardman 13 in United States Pat¬ 
ent 2,108,642 (1938). This device, sche¬ 
matically shown in Figure 12(A), is for 


repeating and amplifying 2-element polar¬ 
ized signals such as those used in teleg¬ 
raphy. The magnetic amplifier com¬ 
prises a single core in the form of a d-c 
saturable transformer having a control or 
saturating winding, a primary winding,, 
and two secondary windings connected in 
a circuit which is the equivalent of a push- 
pull system. The amount of a-c power 
magnetomotive force introduced via the 
primary winding normally is not sufficient 
to saturate the core unless control magne¬ 
tomotive force also exists. 
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Note; ‘Gain is 
effected by caus¬ 
ing negative re¬ 
sistance to appearj 
across points A 
and B. 



Note; 

1. Passes signal fre“ 
quency both ways. 

2. Blocks internvxiu- 
lation products of 
Magnetic Modulator 


Note; 

Resonate* difference 
frequency. 

(This helps to effect 
power transfer* 


o a 

A-C Power 


Figure 16. 


Magnetic amplifier with ultracarrier frequency 
means for hysteresis-effect reduction 


‘agitating’ 


If the control signal is zero, equal and 
opposite potentials of the same frequency 
as the power-supply frequency are in¬ 
duced in the secondary windings. Equal 
and opposite currents thus are caused to 
flow through the rectifiers and through the 
two halves of the output device; the net 
output signal is consequently zero. 

Assume that a positive control current 
is made to flow in the saturating winding. 
The wave shapes of the voltages induced 
in the secondary windings will now assume 
the forms shown in Figure 12(C) because 
of the displaced excitation conditions illus¬ 
trated in Figure 12(B). The rectifier as¬ 
sociated with secondary S A , will pass pos¬ 
itive half-cycles of current, and the recti¬ 
fier associated with secondary S B will pass 
negative half-cycles. Because these rec¬ 
tifiers exhibit a considerably reduced im¬ 
pedance when the applied forward voltage 
is increased, a notably greater current 
flows through the rectifier which is con duct¬ 
ing than through the rectifier which is not 
conducting. Due both to the distortion in 
secondary voltage wave form and to the 
nonlinear resistance characteristic of the 
rectifiers, therefore, a considerable un¬ 
balance of voltage drop is made to occur 
between the two halves of the output 
device. 

The unidirectional currents flowing 
through the output devices also flow in 
secondary windings S A and S B in un¬ 
balanced magnitudes such that a net uni¬ 
directional magnetomotive force exists in 
a direction aiding the original control 
magnetomotive force. This amounts to 
positive internal feedback which in¬ 
creases the amplifier output. 


Negative control-current signals have 
similar, though opposite, effects. 

A magnetic flip-flop based on the prin¬ 
ciple of the circuit of Figure 12, but with 
external feedback added is described by 
E. T. Burton 14 in United States Patent 
2,147,688 (1939). 

Figure 13 is the schematic of a unique 
magnetic amplifier claimed by E. M. 
Boardman 13 in United States Patent 
2,108,642 (1938) and belonging to class B', 
2,b,f,g. Before demodulation the amplified 
signal is an amplitude-modulated wave 
form whose carrier frequency is double 
that of the a-c power-supply frequency. 
This type has been called a “second har¬ 
monic device’’ by some writers. 

The magnetic amplifier of Figure 13 is 
claimed to be particularly adaptable for 
repeating and amplifying 3-element signal 
impulses, that is, positive, zero, and nega¬ 
tive impulses. Two cores in the form of 
d-c saturable transformers are employed. 
Each core has a control or saturating 
winding, a primary winding, and two 
secondary windings connected in the 


Figure 17. Series negative-resistance type magnetic amplifier 


equivalent of a push-pull system. A-c 
power is introduced via the two primary 
windings, and the power-supply frequency 
is considerably higher than any compo¬ 
nent of the input signal frequency. An 
essential feature of the primary excitation 
of this circuit is that it must be sufficient 
to saturate the cores partially during a 
portion of each half-cycle of the power- 
supply frequency. The saturating wind¬ 
ings are connected so that, with zero con¬ 
trol signal, equal and opposing electromo¬ 
tive forces of power-supply frequency are 
induced in them. Decoupling of the con¬ 
trol circuit from the a-c power supply is 
provided in this manner. 

Parameters are chosen so that, with 
zero input signal, potentials of equal and 
opposite wave forms containing a mini¬ 
mum of even harmonic distortion are in¬ 
duced in the seriesed secondaries. The re¬ 
sulting potential across the seriesed 
secondaries is therefore zero under zero 
signal conditions, and no current flows in 
either of the output circuits. 

When a positive current impulse, as 
shown in Figure 13(B), is applied to the 


saturating windings, however, the net 
operating points on the magnetic charac¬ 
teristics of the two cores are shifted so 
that potentials induced in the secondaries 



Figure 18. Buck/ 
boost transformer 
controlled by d-c 
saturable reactor 
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of one core are no longer exactly opposed 
to the potentials induced in the second¬ 
aries on the other core. Instead, a wave 
form with strong second-harmonic con¬ 
tent is induced in both secondary windings 
on each core, and the series combination 
of one secondary on one core and one 
secondary on the other core yields a net 
electromotive force having essentially 
second-harmonic wave form. The dis¬ 
placed excitation conditions are illustrated 
in Figure 13(B), and the individual and 
net effects are illustrated in Figure 13(C). 
The second-harmonic wave form of volt¬ 
age is applied to the demodulators in the 
respective halves of the output circuit. 
One rectifier is connected in the output 
circuit comprising windings and S B in 
such a manner that it will conduct under 
the influence of high-potential positive 
peaks induced in this circuit. Another 
rectifier is oppositely connected in the 
output circuit comprising windings $/ 
and S B ' in a direction so that it will con¬ 
duct under the influence of the negative 
peaks, but not the positive peaks. The 
rectifier elements have nonlinear charac¬ 
teristics (copper-oxide rectifiers, for ex¬ 
ample) in which the forward impedance 
decreases rapidly as forward potential 
increases. The forward current is thus 
much greater through the conducting 
rectifier than is the reverse current 
through the nonconducting rectifier. 
Therefore, because of the combined char¬ 
acteristics of the circuit, the upper half 
of the output device is predominantly 
energized. 


The unidirectional demodulated cur¬ 
rents flowing through the output devices 
also flow through the secondary windings 
in the direction shown by the arrows in 
Figure 13(A). The unidirectional current 
flow in the secondaries produces magneto¬ 
motive forces which in S A and S B aid the 
magnetomotive force due to the positive 
control current applied to the saturating 
windings, but in S A and S B ' the unidirec¬ 
tional current produces a magnetomotive 
force which opposes the control magneto¬ 
motive force. Since the current in wind¬ 
ings S A and S B greatly exceeds that in 
windings Sfi and S B ', however, the net 
effect is that the control magnetomotive 
force is aided by the feedback magneto¬ 
motive force which results from the flow 
of demodulated output current in the 
secondary windings. This action amounts 
to positive internal feedback and operates 
to increase the output of the amplifier. 

Negative current impulses applied to 
the control circuit have similar though 
opposite effects. 

A group of devices similar to that 
shown in Figure 13 are the subjects of E. 
T. Burton’s United States Patents 2,147,- 
688 (1939) 14 and 2,164,383 (1939) 1B which 
include magnetic flip-flops. These de¬ 
vices are not illustrated in this report. 

An important feature of the devices 
shown in Figures 12 and 13 is the inclusion 
of the frequency-discriminating networks 
labelled “net” in the control, power sup¬ 
ply, and output circuits. The purpose of 
the network shown in the a-c power-sup- 
ply circuit is presumably to resonate the 


primary circuit to the frequency of the a-c 
power supply (a series capacitor for the 
same purpose is included in the reactor 
circuit of Figure 6). The purpose of the 
network in the control circuit is to pre¬ 
vent the power-supply frequency har¬ 
monic voltages induced in the saturating 
windings from affecting the control cir¬ 
cuit, and to prevent changes in the im¬ 
pedance of the control line from affecting 
the operation of the amplifier. The. net¬ 
work in the output circuit, Figure 13, is to 
improve the wave forms of the amplified 
signals which are transmitted to the out¬ 
put device. 

A class A,l,g magnetic amplifier, which 
is possibly unique in the art, is shown in 
Figure 14. Described by II. D. Middel 10 
in United States Patent 2,388,070 (1945), 
this device comprises a d-c saturable reac¬ 
tor in which the saturating flux produced 
by the control winding causes a second 
harmonic electromotive force of power- 
supply frequency to be induced in the 
load windings. Under the influence of 
control saturation, an asymmetrical volt¬ 
age is caused to act across a nonlinear re¬ 
sistor seriesed with the load winding of 
each stage. This departure from sym¬ 
metry causes the nonlinear resistor to act 
as a rectifier which introduces a d-c com¬ 
ponent in the load circuit. A minute 
change of d-c signal applied to the satu¬ 
rating winding of the reactor will thus 
cause a much larger change in d-c load 
current. 

The symmetrical nonlinear impedance 
elements in this device have no unidirec- 


Figure 19 (below). Magnetic flip-flop, unbalanced d-c saturable-re¬ 
actor type with cross-external feedback 


Signal 
W inding 


Bias 

Winding 
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tional property or rectifying tendency per 
se. For example, they may be either 
Thyrite or parallel back-to-back rectifiers 
such as copper-oxide rectifiers. The basic 
elements of the amplifier are thus a sym¬ 
metrical current-distorting impedance, an 
asymmetrical current-distorting imped¬ 
ance with polarizing means, an a-c power 
supply, and a d-c load circuit. It may be 
that this device should be considered as a 
general case within which the conven¬ 
tional internal-feedback magnetic am¬ 
plifier, as shown in Figure 10, is a special 
case. The reader is referred to the patent 
for a detailed description of this device. 

Important off-resonant and ferrores- 
onant means for increasing the gain of 
magnetic amplifiers have been described 
from an early date. Figure 15 shows an 
example of the use of such sensitivity in¬ 
creasing means in a (class A,l,a,d,e) volt¬ 
age-regulating device due to L.R. Runal- 
due, 17 United States Patent 2,278,151 
(1942). A parallel inductance-capacitance 
off-resonant circuit is used for load-cur¬ 
rent control. The range of impedance 
obtainable across the terminals of the 
parallel inductance-capacitance circuit as 
it approaches resonance is considerably 
greater than the range which might be ob¬ 
tained by use of the saturable reactor 
alone. The impedance of this network is 
changed by varying the inductance of the 
d-c saturable reactor. If the inductive 
reactance of the reactor approaches the 
capacitive reactance of the paralleled 
capacitor, the network will approach a res¬ 
onant condition and the impedance will 
approach a very high value. If the in¬ 
ductive reactance of the reactor is either 
increased or decreased from the resonant 
value, the impedance of the network will 
decrease. 

The sensitivity of the control circuit in 
Figure 15 may be increased by the inclu¬ 
sion of a ferroresonant combination¬ 
comprising a series-connected capacitor, 
resistor, and self-saturating inductance. 
By use of the ferroresonant combination, 
the current in the control circuit under the 
influence of a varying voltage excitation 
may be made to increase at a much greater 
rate than the excitation voltage. In 
other words, the ferroresonant combina¬ 
tion amplifies the current or voltage varia¬ 
tions occurring in the circuit from which 
it is excited. 

Figure 16 illustrates an unusual and 
little-exploited means for increasing the 
gain of a magnetic amplifier by altering 
the hysteresis effects of the core material. 
In this class B',l,c,h, device, due to R.V. 
L. Hartley, 18 United States Patent 1,287- 
982 (1918), the inventor describes means 
of altering the hysteresis effects in the 
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core material by the application of a 
special agitating field. The core is sub¬ 
jected to a weak high-frequency cross¬ 
magnetization at a frequency several 
times that of the a-c power supply (car¬ 
rier). This patent appears to be based 
first on the principle that a distribution of 
hysteresis loss can be made between two 
frequencies (the agitating frequency and 
the carrier frequency), and second on the 
corollary principle that at low flux densi¬ 
ties the dynamic incremental permeability 
of the material at the low frequency (car¬ 
rier) is considerably increased by the pres¬ 
ence of the higher frequency excitation. 
For a pertinent theoretical treatise on this 
phenomenon, see the article by R. M. 
Kalb and W. R. Bennett. 19 See also G.W. 
Elmen and C. F. Ort, 20 United States 
Patent 1,544,381, 30 June, 1925, which 
applies the phenomenon. 

A significant and novel type of mag¬ 
netic amplifier is the negative resistance 
type shown in Figure 17 and designated as 
class C,l,f. This device was described by 
E. Peterson 21 in United States Patent 
1,884,844 (1932). The operation of this 
circuit depends upon the intermodulation 
of power-supply and signal frequencies in 
such manner that energy is transferred 
from the power source to the signal cir¬ 
cuit. The energy transfer results when a 
difference side band of the two frequencies 
(for example, twice power frequency 
minus signal frequency) is permitted to 
flow in the circuit. See K. Heegner, 22 
United. States Patent 1,656,195 (1928). 
Other intermodulation products must not 
be allowed to flow and thereby to consume 
power. An impedance network allows 
signals of the difference side-band fre¬ 
quencies to flow internally in the device. 
Filters suppress the flow of signals of the 
other generated frequencies except the 
signal frequencies. A general limitation 
of the device is that the bandwidth of the 
signal must be less than the power-supply 
frequency; for wider bandwidths it ap¬ 
pears that the suppression of unwanted 
frequencies will be increasingly difficult. 

Classification 4 of the chart designates 
“special magnetic amplifier circuits.” 
An example of a circuit of this class is 
shown in Figure 18 due to H.E. Young, 23 
United States Patent 2,154,020 (1939). 
This is a class A,1,4, circuit comprising a 
d-c controlled saturable reactor which in 
turn controls a buck/boost transformer 
the secondary winding of which is inserted 
in series with the load. The characteris¬ 
tics of this circuit are such that zero load 
current may be obtained for zero control 
signal provided that the voltage induced 
in the secondary of the transformer is 
equal and opposite to the voltage of the 
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a-c power source under zero-signal condi¬ 
tions. 

Another special circuit is that shown in 
Figure 19 due to M. L. Wood, 24 United 
States Patent 2,524,154 (1950). This is a 
magnetic-amplifier counterpart of the 
Eccles-Jordan trigger pair in which flip- 
flop action is obtained by using part of the 
output signal of one amplifier as the input 
signal to the other amplifier and vice 
versa. In this class A,4,a circuit, the re¬ 
actance winding of saturable reactor 1 
controls (via its associated rectifier) the 
d-c magnetomotive force applied via the 
saturating winding to core 2. Corre¬ 
spondingly, the reactance winding of satu¬ 
rable reactor 2 controls the d-c magneto¬ 
motive force applied via the saturating 
winding to core number 1. This flip-flop 
may be said to be a “cross external-feed¬ 
back type” of magnetic amplifier. The 
state of this flip-flop may be changed by 
applying unidirectional magnetomotive 
forces by means such as the control wind¬ 
ings shown. 

Another type of magnetic flip-flop is 
shown in Figure 20. This is a single¬ 
sided magnetic amplifier similar to that 
shown in Figure 5, but with a sufficient 
amount of external feedback to cause in¬ 
stability. This class A,l,a,c device is de¬ 
scribed by FitzGerald 25 in United States 
Patent 2,027,312 (1936). It has been sug¬ 
gested that the signal-sided flip-flop cir¬ 
cuit of Figure 20 may be used to advan¬ 
tage in that it does not require equal loads 
as does the balanced Eccles-Jordan type 
flip-flop of Figure 19. 
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